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ABSTRACT 


Three lower confidence interval estimation procedures for system ;reliability of co¬ 
herent systems with c> die components are developed and their accuracy measured using 
Monte Carlo techniques. The procedures use either the Poisson approximation to the 
Binomial distribution, the lower Binomial confidence limit procedure, or a modified 
procedure using the Poisson approximation to the Binomial distribution to obtain an 
equation for the lower confidence limit. The accuracy of the interval estimators were 
evaluated using standard computer simulation methods for series, parallel,_ series- 
parallel, and \\'heatstone Bridge systems. The method determined to be most accurate 
can be combined with similar procedures for components that have continuous failure 
times and applied to yield a lower confidence interval procedure for the reliability of 
coherent systems with cyclic and continuously operating components. 
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THESIS DISCLAIMER 

The reader is:cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er¬ 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


Coherent systems are those systems for which the system performs if all components 
function, the system fails if all the components fail, replacing a failed component with a 
working component does not cause the system to fail and similarly replacing a working 
component with a failed component does not cause the system to work [Ref. 1: p. 343]. 
The reliability of a coherent system is not reduced when the reliability of its components 
is increased. Cyclic components are those whose function is measured on a pass-fail 
basis and whose reliability is computed from a discrete probability distribution, usually 
the Geometric distribution. 

The problem of obtaining confidence limits on the reliability of a coherent system 
based on data gathered on its individual components has attracted considerable interest. 
Confidence bounds for the reliability of series systems have been obtained 
asymptotically, based on methods such as Likelihood Ratio, Maximum Likelihood, or 
Modified Maximurn Likelihood. Asymptotic methods are inaccurate at higher 
percentiles unless the component sample sizes are large; Bayesian methods have been 
developed for this problem, but they are extremely sensitive to the selection of the prior 
distribution. [Ref. 2: p. 21] 

Exact confidence limits have been obtained for simple systems. This solution re- 
a.uires that the reliability of at least one of the components has to be greater than that 
of the system. However, identifying such a component becomes quite complex for all 
but simple systems composed of no more than two or three components [Ref. 3; p. 220]. 
Exact methods have also been developed for series systems using as\mptotic approaches 
and the unbiased minimum \ariance estimators of the probability of success, p, based 
on Binomial data [Ref. 4: p. 782]. 

Unfortunately, none of the above referenced interval estimation procedures based 
on discrete data can be readily used in conjunction with the data of components that 
have continuous failure times. Therefore it is difficult to obtain interval estimates for 
the reliability of complex systems that have mixtures of cyclic components and compo¬ 
nents that operate continuously. The methods developed in this thesis can be 
combined with similiar methods that use continuous data, namely those developed by 
Lee [Ref. 5]. The combination of these methods may provide inten’al estimation for the 
reliability of systems with cyclic and continuously operating components. 
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In tins thesis, three procedures that provide lower confidence limits for the reliability 
of coherent systems with cyclic components were analyzed. These procedure<: use only 
discrete data. There is a problem with using component data to establish system reli¬ 
ability, especially for a system that has quite a bit of redundancy. liven though the 
tested components fail individually and their estimated reliabilities arc moderate, had 
these components been assembled into a system, the system could very well have 
worked. In such a case, the system has a very high degree of reliability and methods that 
work well in estimating the more moderate component reliability will not woik well in 
estimating the system reliability. Thus, each procedure that we study has modilications 
to accommodate component test data which when assembed into systems would exhibit 
zero system failures, one systemifailure, or more than one system failure, i.c. 

• no component failures occur or only component.': fail that arc redundant in the 
system, so that no system failure could occur if all the components were combined 
to form systems, 

• exactly one component fails that would result in a system failure or redundant 
components in the system fail in a quantity, so that no more than one system fail¬ 
ure could occur if the components were combined to form systems (this modifica¬ 
tion is explained in ntorc detail in the following chapters), 

• for any component, /, more than one component fails out of /;, tested which would 
lead to more than one system failure. 

The systems analyzed in this thesis arc as follows: 

• all components arranged in scries (.sec Figure 1)_ 



Figure 1. Series System 

• all components arranged in series with at least one component consisting of two 
parallel subcomponents (sec Figure 2) 










Figure 2. Series-Parallel System 


• all components arranged in series with at least one component consisting of two 
parallel subcomponents and at least one component in a two out of three sub¬ 
component parallel arrangement (sec Figure 3) 
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Figure 4. Parallel System 














• a five component redundant system commonly referred to as a Wheatstone Ihidge 
(see Figure 5)^ 



Figure 5. Wheatstone Bridge 
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II. PREFERRED LOWER CONFIDENCE LIMIT FOR SYSTEM 

RELIABILITY 

A. METHODOLOGY 

The procedure selected to determine the lower confidence limit of relatively simple 
complex systems uses a method employed by Lomnicki (Ref. 6: p. 109] and extended by 
Myhre and others [Ref 2: p. 213). It uses the Poisson approximation to the Binomial 
distribution when numerous component failures occur that would result in multiple sys¬ 
tem failures if the tested components were assembled into systems. In cases where zero 
system failures occur or only one system failure could result from the failed components, 
the Poisson approximation appears to be conservative and the lower confidence limit is 
computed using the Binomial distribution directly. 

Suppose a coherent system has k different npes of components that are statistically 
independent. Test data is available on each type of component. The n, tests for com¬ 
ponent i are assumed to be independent Bernoulli trials with probability of failure, q„ for 
each test. Therefore the number of failures, F„ in these n, tests has a Binomial distrib¬ 
ution. 

Any of the k components, say component m, can be chosen as a base component 
and q, expressed as a fraction or multiple, a„ of q„. That is 

qi = aiq„. ( 2 . 1 ) 

Suppose it is appropriate to assume that the probability distribution of the number of 
failures of each component, F„ can be approximated by the Poisson distribution with 
mean n,q,. If ZE, then from Equation (2.1), the distribution of F is approximately 

/kI a 

Poisson with parameter q„Y.n,a,. Consequently the mean of F, E(F), is given by 

k 

E{F) = qn^n-fl; (2.2) 

i=l 

The system reliability, R„ can be defined as a function of the unreliability, q„ and the 
q, 

associated values a, = That is, 

Hm 

= (2-3) 
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By definition, h is non-increasing in q„, because the system is coherent. If a„ 
are known, an approximate lower confidence limit, niay be obtained from an up¬ 
per confidence limit, for by the equation 

A ^ 

• ^ 2 >—. ( 2 - 4 ) 


When appropriate, the upper confidence limit, may be obtained from the well 
known upper confidence limit for the mean, x, of a Poisson distribution, namely, ifFis 
distributed P01SS0N(/.) then 


/ 




X 


2 

a, 2(\+F) 


2 


(2.5) 


where is the (I - a) quantile of a Ghi-square probability distribution with degrees 

of freedom equivalent to 2(1 + F), where F is the number of system failures. From 

a’ 

Equation (2.2), substitution of q^T,n,a, for 2 into Equation (2.5) gives 

imi 


?m,t;(a) 


■v2 

■^s(,2(l+f) 


k 



1=1 


( 2 . 6 ) 


If the a, are unknown then an approximate upper confidence limit, is given by 


A 


X 


2 

3, 2<l+f) 



(2.7) 


where a, is an estimator for a„ i= 1,2,...,k. 

The Poisson approximation to the Binomial distribution is conservative when all F, 
equal zero or redundant components of the system fail in such a way that results in zero 
system failures. In this case, let n* represent the number of system tests equivalent to 

testing each component n, times, / = I k. Then the distribution of s, the number of 

system tests that would have been successful, is approximately Binomial(rt*,/J,). If s is 
distributed BINOMIAL(«*,i?,) then we can use the binomial lower confidence limit 


- 



( 2 . 8 ) 


6 






10 compute the lower confidence limit for system reliability. The following two methods 
for the calculation of the number of equivalent system tests, w’ere selected and each 
applied separately using Equation (2.8). 

«,* = min(H,,H2,-. «a) (2-9) 

k 

Z- 

« 2 *=-^ ( 2 . 10 ) 

w'here n, is the number of tests of component /. 

Some instances of component failures could only result in one equivalent system 
failure if all of the tested components were combined into complete systems. In this 
case, we again define and treat system trials with one failure. The resulting lower 
confidence limit, for system reliability is the solution for p = {\- q) in the equation 

»= Z (”■)(>-wV”-' (2.I1) 

In the Wheatstone Bridge case, two or more component failures among the tested 
components could result in one sjstem failure if the components are assembled into 
Wheatstone systems. In this case we set n* = n* for one interval procedure and 
n* = n* for the second interval procedure. 

In a series system, tv' is equal to the number of tests performed on the failed single 
component, because the reliability of the system is determined largelj by the reliability 
of the least reliable component. Since it is difficult to solve for p in Equation (2.11), an 
equivalent equation using the Snedecor F distribution is used. Thus, 

^s, L{x) ~ T+TT+Tif (2-12) 

where s is the number of system successes, / is the number of system failures, and 
F(|..,.jj^.„j, is the a quantile of the Snedecor F distribution with 2{f + 1) and 2s degrees 
of freedom [Ref. 7; p. 43). 











B. RESULTS 

The accuracy of.this procedure was evaluated using computer simulations for each 
oTthe following systems described in Section A; 

• series systems 

• series sytems where the second component is composed of two parallel subcompo¬ 
nents 

• series systems where the second component is composed of two parallel subcom¬ 
ponents and the fourth component consists of a two of three subcomponent par¬ 
allel arrangement 

• parallel systems 

• Wheatstone Bridge 

Groups of test data were generated where the parameters, q, and n, were chosen to 

A 

control the expected number of failures, £[/l = Confidence levels of 0.20 and 0.05 

i=l 

were used in each case. A total of 1000 replications were generated for each set of pa- 

A A 

rameter values. Each replication produced one value of R,These 1000 values, R, 

A 

were ordered and used to get the simulated probability distribution of The simi’- 
lation procedures are described in Chapter IV. The 80''’ and 95''' percentile point of the 

A 

simulated probability distribution of was compared to R, for determining the ac- 

A 

curacy of the procedure. This comparison is made because R,,i.(„ is the lower 
• * 

100(1 — a) percentile confidence limit for R, if 1 - a = R(R,.l(.) ^ R/)- This equation states 

A 

that R, is the 100(1 — a) percentile point of the probability distribution of R,,lw The 
"true confidence level" is the percentile point of the simulated distribution corresponding 
to-the true value of R,. 

The parameter values determine a case number and are labeled as 

such in the tables that describe the simulation results. A summarx table that provides 
the parameter values, q, and n,, is given in Appendix B. 

All tables report the SO"" and 95"' percentile points of the simulated distribution of 

A A / 

R/.i(a) and appear in the tables under the column labeled R,x(.)- 
1. Series System 

By definition, the reliability, R„ for a series system of k independent components 
is 

k k 

(=i /=i 
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The corresponding lower confidence liniii is given by 

~ ^ “ ^i^m,U{3)) ( 2 -' 

i=l 

These formulae are used to calculate the reliability when at least two components fail. 
If zero components fail Equation (2.8) is used and if one component fails Equation 
(-2.12) is used. The results are presented in Table I. In Table 1, column 1 of is 
calculated using Equation (2.10) and column 2 is calculated using Equation (2.9) for 
when the component failures equate to zero system failures. 


Table I. SERIES SYSTEM 


Case 

U 

Compts 

E[_F] 


a 

Level 


True Confidence 
Level 

1 

5 

1.03 

.95572 

.20 

,91447 

.90345 

75 

75 

.05 


.85140 : 

BIBB 

100 

■ 

10 

1.2 

.93206 

.20 

.88491 


100 

100 

.05 


.85140 . 

100 

100 

3 

5 

4.9 

.95034 

.20 

.95361 

.95361 

76 

76 

.05 

.95969 

.95969 

90 

90 

■ 

10 

B 

.85951 

.20 

.85369 

.85369 

85.2 

85.2 

.05 

.86718 

.86089 

94.5 

KSl 

5 

5 

5.74 

.95084 

.20 

.95036 

.95036 

76.3 

76.3 

.05 

.95779 


91.1 

91.1 

6 

To 

5.75 

85951 

.20 

.85910 

.859i0 

81' 1 

80.1 

.05 

.88114 

.88072 

92.5 

92.6 

■ 

5 

10.5 

.85828 

.20 

.87174 

.87174 

72.4 

72.4 

.05 

.88366 

.88366 

86.6 

8- ^ 

8 

10 

10.5 

.85951 

.20 

.87137 

.87137 

73.3 

73.3 

.05 


.89049 

86.4 

86.4 


In some cases the true confidence level equals 100. This is not uncommon for confi¬ 
dence intervals based on discrete data and is further explained in Appendix A. 
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2. Series-Parallel Systems 

A modification was made to the series systems to form the series-parallel system. 
The second component in the series was modified to consist of two parallel subcompo¬ 
nents of equal unrelibility, . For the second component to fail both subcomponents 
must fail. The unreliability of component two is Since component two is in series 
with,the other components the corresponding lower confidence limh is calculated using 
Equation (2.8) when the equivalent number of system failures is ?,ero, Equation (2.12) 
when the number of system failures is one, and Equation (2.14) is used in all other cases. 

A / 

The results are presented in Table 2. In Table 2, column I of is calculated using 
Equation (2.10) and column 2 is calculated using Equation (2.9) for n’* when the com¬ 
ponent failures equate to zero system failures. 


Table 2. SERIES-PARALLEL SYSTEM 


Case 

n 

Compts 


D 

a 

Level 


True Confidence 
Level 

9 

5 

1.03 

.95572 

.20 

.91447 

.90345 

75 

75 

.05 

.85140 

.8^140 

iuO 

100 

10 

10 

2.72 

.93206 

.20 

.'•5770 

.95770 

64.5 

64.5 

.05 

.95344 

.95344 

74.4 

64.5 

11 

5 

B 

.85828 

.20 

.85369 

.85369 

83.5 

83.5 

.05 

.86685 

.86089 

94.2 

94.3 

12 

10 

4.9 

■ 

.20 

.85546 

.85546 

83.1 

83 1 

.05 

.86685 

.86635 

94.1 

94.1 

13 

5 

9.79 

.95084 

.20 



76.S 

76.8 

.05 

.95779 

.96029 

88.6 • 

88.6 

14 

10 



.20 

.85943 

.86943 

75.2 

75.2 


.05 

.88606 

.88606 

87.6 

87.6 


Note the accuracy of the lower confidence interval is somewhat reasonable when the 
expected number of failures exceeds 4.9 in both the series and series-parallel .systems. 

3. Series-Parallel Systems with a 2 of 3 Parallel Component 

Another modification was made to the five and ten component series systems. 
The third component in the series was modified and now consists of three parallel sub¬ 
components of equal unreliability, g^^. For this component to fail, two or three parallel 
components must fail. The unreliability, gj, of component three is (1)(1 + 
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The fourth component of each sjstem ;s composed: o^ two parallel'subcomponents as 
defined in the series-parallel system. Since these components are in series with the other 
components, the corresponding lo\\er=confidence limit is calculated using the same series 
of equations as the series-parallel system. T. ••esults are presented in Table 3. The 
term /i* is computed from Equation (2.10) wir » equivalent number of system-fail¬ 
ures is zero. 


Table 3. SERIES-PARALLEL WITH A 2 u jT OF 3 COMPONENT SYSTEM 


Case 

# 

Compts 


R. 

II9 


True Confidence 

I evel 

H 

5 

.88 

.96525 

.2C 

.91447 

100 

.05 

.85140 

100 

16 = 

10 

2.6 

.94136 

.2\> 

.95770 

63.7 

.05 


85.7 

17 

5 

3.21 

.96035 

.20 

.95937 

80.4 

.05 

.96877 

88.5 

18 

10 

mi 

.88945 

.20 

.88582 

80.9 

.05 

.91735 

89.5 

■ 

5 

6.17 

.96035 

.20 

.96126 

78.8 

.05 



20 

10 

8.37 

■ 

.20 

.89734 

76.7 

.05 

.91817 

87.8 


Another series of simulations were conducted on these systems. For each case 
previously run, new cases were defined by decreasing the unreliability of selected com¬ 
ponents without changing the number of component tests. These changes result in a 
decrease in the number of expected failures and an increase in the reliability of the sys¬ 
tem. The results are presented in Table 4, The term n* is computed from Equation 
(2.10) when the equivalent number of system failures is zero. 
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Table 4. SERIES-^^MRALLEL WITH A 2 OUT OF 3 COMPONENT SYSTEM 
(CONSTANT NUMBER OF MISSION TESTS, DECREASING RELI¬ 
ABILITY’) _ 


m 

W 

Compts 

Em 

’ 

R. 

Level 

s / 

Frue. Confidence 
Level 

21 

10 

.51 

.99136 

.20 

.97661 


.05 

.95690 

bhe^hbi 

22 

10 

■ 

.96836 

.20 

.97661 

72.3 

.05 

.95690 

96.3 

23 

10 

6.57 

.95430 

.20 

:95497 ‘ 77.2 

.05 

.96169 ' 91.3 

24 

10 

9.90 

.90249 

.20 

.90515 

78.9 

.05 

^91486 

91.6 


Tables 3 and 4 indicate that for £(/0 > 5 the lower confidence limits are reasonably ac¬ 
curate. 

4. Parallel System 

The accuracy of the lower confidence procedure was evaluated when it was ap¬ 
plied to a five component systc ' with each component, i, in parallel with the others. 
The unreliabilities, q, of each component are equal and the unreliability v-.f the system is 
defined as qK For the system to fail all five components must fail. The ress.!.;^ of the 
computer simulations arc presented in Table 5. The term is computed from Equation 
(2.10) when the equivalent number of system failures is zero. If n of these components 
are tested and four or fewer failures occur, then no system failures would ha\e occured 
had hese components been assembled into systems. The lower confidence limit is cal- 
cula '.d using Equation (2.8) when the equivalent number of system failures is zero, 
Equation (2.12) when the number of system failures is one, and Equation (2.14), '.vith 
/: = 1, in all other cases. 







































Table 5. PARALLEL SYSTEM 


; Case 


£[F] 

K 

a 

Level 

/ 

True Confidence 
Level 

; 25 

5 

■ 

.99000 

.20 

.98403 

100 

■■ 

.97049 

100 

; 26 ; 

5 

4 

.96000 

.20 

.95770 

91.1 

.05 

.95344 

98.3 

; 27 : 

5 

: 

6 

.94000 

.20 

.93370 

83.9 

mm 

.93838 

98.2 

; 28 

5 

9 

.94000 

.20 

.94012 

77.6 

.05 

.94002 

93.6 


5. Wheatstone Bridge 

By definition, the reliability, P,, for the Wheatstone Bridge, in Figure 5, with five 
independent components is 


= P\Pi + P 2 P 5 + P\P2l>5 + P 7 P 3 PC - P\PiP2Pi p j 

-P\P2PiP5~P\P3PiP5-P2P3PiPs-i-'^PiP2p3PiPs ^ 

where /?, = (I — q). In terms of q, Equation (2.16) becomes 

«i.-. 05 ) 

= 1 - -b 

4 \ I*" * ' / 

- 2<7^(aiC203Cias) 

[Ref. 2: p. 215). By substituting a,q^.,m„ for q„ the corresponding lower confidence limit 
U obtained. 


^s. Ufl) - K^m, 0'(a)' %"•< “s) 


( 2 . 18 ) 
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The reliability of the Wheatstone Bridge system is normally high due to the re¬ 
dundancy of the system. This system experiences zero .ystem failures in the following 
five difTerent failure patterns. 

• F, = Fj = F3 = f, = = 0 

• F, = F, = F, = 0 

• F, = F, = 0 

• F, = Fs = 0 

. F3 = F3 = F, = 0 

w'here F, is the number of failures of component i, among its tests. Any other failure 
patterns will produce one or more system failures. The lower confidence limit is calcu¬ 
lated: using Equation (2;8) when the equivalent number of system failures is zero, 
Equation (2.12) when the number of system failures is one, and Equation (2.18) in all 

other cases. The results of the computer simulations are presented in Table 6. In Table 
^ / 

6, column 1 of F,x(,j is calculated using Equation (2.10) and column 2 is calculated using 
Equation (2.9) for when the component failures equate to zero system failures. 


Table 6. WHEATSTONE BRIDGE SYSTEM 


Case 

n 

Compts 


R, 

a 

Level 

» / 

Ri.ut) 

True Confidence 
Level 

29 

5 

5.15 

.99776 

.20 

.99658 


88.6 

88.6 

.05 

.99766 

.99166 

95.4 

95.4 

30 

5 

6 

.99976 

.20 

.99960 

.99960 

91.4 

91 

.05 

.99966 

.99966 

98.3 

98.3 

31 

5 

.5 

.mil 

.20 


.72478 

100 

100 

.05 

.98678 

.98678 

100 

100 

32 

5 

.5 

.99977 

.20 

.92622 

.72478 

100 

100 

.05 

.99339 

.99339 

100 

100 


The system reliability values in cases 31 and 32 are too large for these interval estimation 
procedures to be accurate using the sample sizes given in Table 15 in Appendix B. Also, 
more than 1000 replications would be needed to assess the accuracy of any lower confi¬ 
dence limit procedure for system reliability when the true system reliability is as large 
as it is in these cases. 
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in. ALTERNATE PROCEDURE A FOR THE LOWER CONFIDENCE 
LIMIT FOR SYSTEM RELIABILITY 

A. METHODOLOGY 

This procedure determines the lower confidence limit for simple and complex sys¬ 
tems using only the Poisson approximation to the Binomial distribuiton. It is used ex¬ 
tensively by Bellini [Ref. 8: p.4-6]. 

This procedure uses the principles outlined in Chapter II of this thesis except in the 

cases where there are zero or one system failure. In the case where there are zero system 

failures (no components fail) the estimated unreliabilty, q„ becomes zero because 
F. 

q, = where F, is the number of failures in n, mission tests of component /. Therefore, 
the value of a„ the estimated value of a„ becomes zero and Equation (2.7) becomes un¬ 
defined. When this occurs the estimated lower confidence limit of the sjstem reliability, 

A 

R,is defined as 






(3.1) 


where n* is defined as in Equation (2.10). 

If only component, m, fails the value of is equal to 1. The a, for all other com- 

A 

poncnis are zero, and the lower confidence limit, is defined as 




^3.2(i+n 


2/r 


'm 


(3.2) 


w'here ti„ is the number of mission tests of the failed component. 

B. RESULTS 

The accuracy of this procedure was evaluated only for series and Wheatstone Bridge 
systems. Testing was limited because a comparision of results of the three procedures, 
discussed in this thesis, indicated this procedure to be less accurate in determining the 
lower confidence limit for the reliability of simple systems. This observation is noted in 
those systems experiencing zero system failures or one system failure. Selected results 
illustrating the accuracy of this procedure compared to the accuracy of the '’preferred 







procedure" (Chapter II) are indicated in Tables 7 and 8. The term n* is computed from 
Equation (2.10) when the equivalent number of system failures is zero. 


Table 7. COMPARISON OF PREFERRED PROCEDURE AND ALTERNATE 
PROCEDURE A (SERIES SYSTEM) 



Preferred Procedure 

Alternate Procedure 

Case 


a 

Level 

D * 

True Confi¬ 
dence Level 

R ^ 

True Con¬ 
fidence 
Level 

1 

.95572 

.20 

.91447 

75 

.83140 

100 

.05 

.85140 

100 

.92356 

100 

2 

.93206 

.20 

^ .88491 

100 

.83140 

100 

.05 

.85140 

100 

.92356 

100 

3 


.20 

.95361 

76 

.95361 

76 

.05 

.95969 

90 

.95969 

90 

■ 

.85951 

.20 

.85369 

85.2 

mam 

85.5 

..05 

.86718 

94.5 

.85083 

95.2 

5 

.95084 

.20 


76.3 

.95036 

80.3 

.05 

.95779 

91.1 

.96359 

89.6 

6 

.85951 

.20 

.85910 

80.1 

.86910 

SO 

.05 

.88114 

92.5 

.88072 

92.9 

m 

.85828 


mism 

72.4 

.87174 

72.4 



86.6 .88366 

86.6 

8 

.85951 

.20 

.87137 

73.3 

.87137 

73.3 

.05 

.89049 

86.4 

.89049 

86.4 


I 


16 










































































Table 8. COMPARISION OF PREFERRED PROCEDURE AND ALTERNATE 
PROCEDURE A (WHEATSTONE BRIDGE) 



Preferred Procedure 

.Alternate Procedure 

Case 

R, 

a 

Level 


True Confi¬ 
dence Level 

* / 

R>,lm 

True Con¬ 
fidence 
Level 

29 

.99776 

.20 

.99658 

88.6 

1.000 

63.6 

.05 

.99766 

95.4 

1.000 

70.4 

30 

.99976 

.20 

.99960 

91.4 

1.000 

66.2 

.05 

.99966 

98.3 

1.000 

73.1 

31 


.20 

.91447 

100 

.72478 

91.8 

.05 

.98678 

100 

l.OOO 

91.8 

32 

.99977 

.20 .92622 

100 

.99549 

95.7 

.05 .99339 

100 

1.000 

93.7 


This alicrnaic procedure produced lower confidence limits that were more conservative 
than those produced by the “preferred procedure”. 
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IV. ALTERNATE PROCEDURE B FOR THE LOWER CONFIDENCE 
LIMIT FOR SYSTEM RELIABILITY 

A. METHODOLOGY 

The purpose of this alternate procedure was to construct a procedure that accom¬ 
modates zero component failures in a different manner than that employed in other 
procedures. This procedure redefines the estimated value, a„ of and in so doing uses 
the Poisson approximation to the Binomial distribution. 

Suppose a component, i, undergoes //, tests where there is a probability of success, 
on each component test. Let F, be the number of failures in these n, tests. We define 

ia\ 

The 50 percent lower binomial confidence limit for the reliability of component /, 
p,^i, can be determined by using its number of failures, F„ and number of tests, //,. The 
upper Biniomial confidence limit, is equal to 1 - In this alternate procedure we 
define a, by 


^m,U 

where = max(^,,„, •••> ^kv)- Note that the index will be determined by the data 

and a, is well defined even if no components fail. The probability distribution of F, is 
approximated by the Poisson distribution. The estimated upper confidence limit, 
can be calculated using Equation {2.1) and the system reliability lower confidence limit, 
/.(.)» can be obtained from Equation (2.4) where a, is defined by Equation (4.1). 

This procedure is used without exception, regardless of the number of system fail¬ 
ures. 

B. RESULTS 

This procedure was evaluated only on series and Wheatstone Bridge systems. 
Evaluations were limited because a comparision of the results with other procedures 
discussed in this thesis, indicated this procedure to be less accurate. Selected results il¬ 
lustrating the accuracy of this procedure compared to the accuracy of the "preferred 
procedure" (Chapter II) are indicated in Tables 9 and 10. The term m* is computed for 
Equation (2.10) when the equivalent number of system failures is zero. 
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Table 9. GOMPAR, .‘N OF PREFERRED PROCEDURE AND ALTERNATE 
PROCEDU L B (SERIES SYSTEM) 





Preferred Procedure 

Alternate Procedure 


H 


a 

Level 

p ^ 

True Confi¬ 
dence Level 


True Con¬ 
fidence 
Level 

i 

■■ 

.95572 

.20 

.91447 

75 

.97m 

37.1 

i 

■ 


.05 

.85140 

100 

.99489 

7.4 

j 

2 

.93206 

.20 

.88491 

100 

.97197 

70.3 

i 

1 



.05 

.85140 

100 

.99350 

33.2 


: 3 

.95084 

.20 

.9536 L 

76 

.97294 

25.7 

1 



.05 

.95969 

90 


9.7 


4 

.85951 

.20 

785369 

85.2 

.93224 

26.7 




..05 


94.5 

.97896 

8.3 


5 

.95084 

.20 

.95036 

76.3 

.97497 

30.1 




.05 

.95779 

91.1 

.99155 

9.2 


6 

.85951 

.20 

.85910 

80.1 

.91470 

39 




iBSI 


92.5 

.97405 

7.1 


7 

.85828 

.20 

.87174 

72.4 

.89848 

44.5 




.05 

.88366 

86.6 

.94607 

19.4 


8 

.85951 

.20 

.87137 

73.3 

.88412 

56.3 




.05 

.89049 

00 

.93886 

24.9 



Cases 4, 5, and 6 in Table 9 clearly illustrate a more accurate "preferred procedure" for 
a series system. 
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Table 10. COMPARISION OF PREFERRED PROCEDURE AND ALTERNATE 
PROCEDURE B (WHEATSTONE BRIDGE) 


- - 

Preferred Procedure 

Alternate Procedure 

Case 

R. 

a 

Level 

if* 

■“>. 1 ( 1 ) 

True Confi¬ 
dence Level 

Ri,LM 

True Con¬ 
fidence 
Level 

29 

.99776 

.20 

.99658 

88.6; 

98123 

90.1 

.05 

.99766 

95.4 

98235 

97.6 

30 

.99976 

.20 

.99960 

91.4= 

99810 

92.0 

.05 

.99966 

98.3; 

99856 

99.1 

31 

.99977 

.20 

.91447 

: 100 _ 

^90478 

91.8 

.05 

.98678 

100 

.99900 

91.8. 

32 

.99977 

.20 

.92622 

100 

90932 

100 

.05 

.99339 

100 

99221 

100 


In the case of the Wheatstone Bridge, Table -10 illustrates that the "preferred procedure" 
is more accurate than the "alternate procedure". 
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V. SIMULATION 

Standard simulation -techniques are used to determine the accuracy of the lower 
confidence^limit procedures. Some of the basic simulation programs were developed by 
Bellini [Ref. 8: Appendix A]. Each program was modified to incorporate the necessary 
mathematical formulae that are needed to define the lower confidence limit for a partic¬ 
ular procedure. 

Input parameters needed to;run the computer programs are 

• k - number of components in the system 

• n - vector of component tests (/ii./ii...., wO 

• q - vector of component unreliabilites 

• 0. - level of confidence 

Note that the value of system reliability, is determined by the vector q. 

The NON-IMSL random number generator, SRND, was used to simulate the suc¬ 
cess or failure of each test of the k components. From this data, the values of q„ a„ 
q„,vifp and are calculated. Each scenario was replicated 1000 times to generate 
1000 ordered values of The 1000(1 — a)'* ordered value, ^,.i(,),iooo(i-.) front smallest 
to largest denotes the 100(1 - a)''’ percentile point of the probability distribution of 

A A 

R,,n,y If the lower confidence limit procedure is exact, lowo-.) should equal R,. The 

A 

percentiles in all of the tables are the values of R,,u,).\mi.\->Y A true confidence limit is 
then determined by finding the element of the vector of replications which is closest to 
R, and noting its index number, j . The true confidence level is then calculated to be 

1000 ^ 

. A system report is generated and reports to the analyst the following information. 

• qi - unreliability of each component, / 

• a, - fraction of unreliability of base component, m 

• n, - number of mission tests for each component, / 

• R,- true system reliability 

A 

• R,,L(t) - estimated lower confidence limit for the 100(1 — a) percentile 

• difference between R, and R,,l(,) 

• true confidence level 
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VI. CONCLUSIONS AND RECOMMENDATIONS 

The accuracies of three approximate interval estimation procedures, based on dis¬ 
crete component data, for the reliability of coherent systems were analyzed in this thesis. 
Computer simulations were used to perform this analysis. Each interval: estimation 
method emi-». , o ratios of component failure rate estimates when two or more different 
component types have at least one-failure. This specific ratio feature is needed'to extend 
this work to more complex systems with mixtures of cyclic components and 
componenets whose failure times have a continuous probability distribution. 

The simulations reveal that the methodilabeled the "preferred method" in this thesis 
appears to be reasonably accurate if four or more failures are expected to occur among 
all components tested. However, any general interpretation of this t\pe is not warranted 
at this time. The potential for error can be significant if several components have rela¬ 
tively small samples (less than 15) and zero failures. Zero failures joined with small 
sample sizes will always be the bane of classical interval estimation procedures. 

The ratio procedure does allow the possible use of information extraneous to the 
data. Previous test programs on similar hardware operating under similar environments, 
as that present for the current test data, might be used to modify the component failure 
rate estimates or prehaps only the ratios of the failure rate estimates. So called "off the 
shelf' hardware purchased in accordance with existing DOD specifications would be 
prime candidates for this type of failure rate modification. Supplementing current test 
data with other existing "similar" data has become more common as resources for reli¬ 
ability demonstration testing has been reduced. 
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APPENDIX A. DISCRETE CONFIDENCE LIMIT PROPERTIES 

Equations for confidence limits on parameters of discrete probability distributions 
are not exact. If Aw is the lower 100(1 — a) percent confidence limit for the parameter 
p in the Binomial distribution, then is defined so that 

1-a (^1) 

If the parameter p, is the probability of success on each trial in a sample of size n 
and 5 is the obsetA'ed number of successes then Aw is the solution for p in the equation 

n 

-Pf~^ = « (^2) 

y=-f 

if s>0, and Pu,) = 0 if s = 0. Specifically, suppose s = n then in Equation (A2), 
Aw = • ^i^is is largest value of Aw- Consequently if the true value of p is 

greater than then P{pu,) ^p)= !• This has important implications when analyzing 
computer simulations of confidence limit procedures based on discrete data to assess 
their accuracy. If the value of p used to generate the data on the computer is greater 
than then all of the p^,^ values vyill be smaller than p and the analysis will show the 
procedure has confidence level 100 . This is to be expected when evaluating these con¬ 
fidence interval procedures for some choices of sample sizes and parameter values. 

■ The exact value of P{pu,y ^ p) depends on the sample size, n, and the true value of 
p. For fixed n, the possible values of s are 0,l,2,...,rt. Each value of s yields a specific 
value of Awi say p{s). Consequently, 

PiPu,) = P{s)\p) = P{S = s\p) (A3) 


and 


= (A4) 

If the true value of p equals AwW ^ ^hen, the probability in Equation (A4) 

has the value 1 — a because this value of p satisfies Equation (A2). Consequently if the 
true value of p equals any of the values p{ti), p(n — 1),..., p(\), then P{pu,) <p) = 1 — a. 
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For all other values of;?, Pipn,-^ </))>!— a. Figure 6 is a sketch of the behavior of this 
phenomena for a = .10 and /; = 5. 


1 

P(Puo<PlP) 

.90 


0 .1 .3 .5 .7 -9 1 

P(2) p.(3) p{4) p(5), 


Figure (i. Behavior of F(/4(,) < p) 
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APPENDIX B. INPUT PAPji.METERS 


Table 11. SERIES SYSTEM INPUT PARAMETERS 


Case 

Components 



I 

2 

3 

BB 

5 ^ 

6 

BB 

8 

9 

10 

B 


.02 

.01 

.005 

wsm 







n, 

30; 

25 

20 

10 

; 5 


B 




fl 

(Ji 

.02 

.01 

.005 

.005 

.005: 

.005 

.005 

.005 

.005 

.005 

n, ^ 

30= 

25’ 

20 

wm 

■I 

5 

5 

5 

5 

5 

: 3 5 


.01 



'.008 

.025 







200 

400 


265 

mm 

S 





B 

mm 

.005 

.0_1 

.015 

.02 

mm 

.005 

.01 

.015 

.02 

.025 

BB 

100 

ma 

30 


20 

100 

50 

30 

25 

’ 20 

B 


.01. 

.005 

.003 

.008 


IBl 





n, 

150 

80 

240 

265 

50 






6 

mm 

WiMM 

.01 

.015 

.02 

.025 

.005 

.01 

.015 

.02 

,025 

n, 

150 

20 

20 

. 50 

wem 

150 

20 

20 

50 

25 

B 


•01 

■B 

.03 


.05 ^ 



- 



ih 

300 

20 

20 

100 

50 






8 

mm 

.005 

.01 

.015 

.02 

.025 

.005 

.01 

.015 

.02 

.025 

tl; 

eb 

20 

20 

100 

50 

300 

20 

20 

100 

50 


25 









































































































































































































































Table 13. SERIES-PARALLEL WITH A 2 OUT OF 3 
INPUT PARAMETERS 


Components 


COMPONENT SYSTEM 





7E-5 2E-5 


005 .005 .005 


10 10 10 10 






2E-5 6E-5 .025 


240 265 50 


7E-5 4E-5 .025 


20 50 25 


2E-5 6E-5 .025 



.01 .015 .02 .025 


20 20 50 25 








.005 

.01 

7E-0 

4E-5 

.025 

.005 

.01 

.015 

.02 

.025 

«,■ 

300 

20 

20 

100 

50 

300 

20 

20 

100 

50 

IHI 

lE-4 

mm 

7E-5 

BH 

.001 

.001 

.001 

.001 


4E-4 

n,- 

200 

100 


50 

40 

40 

40 

40 

50 

60 

m 

.001 

.002 

7E-5 

.000 

.005 

.005 

.005 

.005 

.006 

.003 


IBII 

100 

OrT" 

50 

40 

40 

40 

40 

50 

60 

<7/ 

.025 


.001 

.001 

.001 

.001 

.002 

-.005 

n, 

200 

100 

60 

50 

40 

40 

40 

40 

50 

60 

(J; 

.025 

.01 

7E-5 

4E-4 

.005 

.005 

.005 

.005 

.006 

.035 

«/ 

200 

loo"’ 

60 

50 

40 

40 

40 

40 

50 

60 





























































































































































































































































APPENDIX G. FORTRAN CODE FOR THE PREFERRED LOWER 
CONFIDENCE LIMIT FOR SYSTEM RELIABILITY (SERIES SYSTEM 
AND WHEATSTONE BRIDGE SYSTEM) 

PROGRAM ZFYSCN 


* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE * 

* ZERO FAILURES ALLOWED; NO SCALING * 

* AUTHOR: E. F. BELLINI, LT, USN * 

* MODIFIED BY: LT VALERIE A. COVINGTON,USN (MAR 90) 

* DATE: NOV 89 * 

it it 

* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 

* RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 

* OF THEIR COMPONENTS * 

it it 

* IN ITS PRESENT CONTIGUEATION THIS PROGRAM IS SET UP TO RU-N 12 * 

* TIMES EACH TIME PRODUCING 1000 REPLICATION’S USING A DIFFERENT * 

* SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC'.* 

* THE REXX EXEC PROGRAM * 

* 'Bl' CALLS THE INPUT FILES TO BE READ ANT) NAMES THE 12 OUTPUT * 

* FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE * 

* INT)EX COUNTERS I, J, K OF THE 'Bl' EXEC ON'E CAN RUN ANT USER- * 

* SPECIFIC RUN FROM JUST ONT RUN TO ALL 12. * 


* VARIABLES USED * 

* * 


it 

AHATI 

it 

AI 

•it 

ALFA 

•it 

BIGF 


CHISQ 


C1C15 


DEGFR 

•it 

DELBRG 

•it 

DELSTR 

•it 

DELTAR 

* 

DIFF 

* 

EPS 

* 

ERROR 

•it 

FAILS 

. 2 . 

FI 


FLAG 

* 

INC 

* 

KEYl 

* 

KEY2 

* 

KEYS 

* 

KEY4 

* 

KK 

* 

LOOP 

it 

MAXALF 


WEIGHT ESTIMATES FOR EACH COMPONENT * 
INPUT WEIGHTS FOR EACH COMPONENT * 
LEVELS OF SIGNIFICANCE * 
TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHI-SQUARE RANDOM VARIABLE VALUE * 
FORMAT LABEL * 
DEGREES OF FREEDOM * 
DIFFERENCE FOR BRIDGE SYSTEM * 
DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DIFFERENCE FOR SERIES SYSTEM * 
DIFFERENCE (TRUE REL. - ESTIMATED REL.) * 
SMALL QUANTITY(CONSTANT) * 
PARAMETER FOR IMSL ROUTINE * 
COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INCREME.NT STEP SIZE FOR ROUTIN'E USMNTIX * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTIN’E SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARR.4Y SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
COUNTS NO. OF REPLICATION PERFORMED * 
MAX NO. OF SIGNIFICANT LEVELS DESIRED(ARRAY SIZING)* 
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* 

MAXREP 


MAX NO. OF REPLICATIONS 

k 

* 

MAXHUN 


MAX NO. OF PROGRAM ITERATIONS ALLOWED 

k 

* 

MSTRQ 


MASTER U.NRELIABIEITY(USED WITH Al'S TO CALC. Ql’S) 

k 

it- 

MULT 


MULTIPLIER FOR RANDOM NO. GENERATOR SRND 

k 

* 

N 


NO. OF MISSION TEST FOR EACH COMPONENT ’ 

r 

•k 

NIMAX 

• 

MAX NO. OF MISSION TESTS 

k 

•k 

NIMIM 


MIN NO. OF MISSION TESTS 

k 

* 

NINDX 


INDEX NO. OF MAX NO. OF MISSION TESTS 

k 

IT 

NIREAL 


NO. OF MISSION TESTS TRANSFOR.MED TO REAL 

k 

•fe 

NMAX 


MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL- 

k 

k 

NPRNT 


FLAG FOR DETAILED REPORT OUTPUT 

k 

k 

PRNT 


SAME AS ABOVE(PARAMETER) 

k 

k 

QHATI 


UNTiELIABILITY ESTIMATES FOR EACH COMPON’ENT 

k 

k 

QHTMAX 


LARGEST QHATI 

k 

k 

QHTUPR 


UPPER LIMIT ON SYSTEM UNTIELIABILITY 

k 

k 

QI 


INPUT UTiRELIABILIY FOR EACH COMPONENT 

k 

* 

QINDX 


INDEX 

k 

* 

QUANTL 


QUANTILE 

k 

* 

REPSHD 


REPLICATIONS HEADING FORMAT NTJMBER 

k 

* 

RHTSTR 


SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 

k 

k 

RS 


TRIjE SERIES SYSTEM RELIABILITY 

k 

k 

RSBRDG 

* 

TRUE BRIDGE SYSTEM RELIABILITY 

k 

k 

RSHAT 


SERIES SYSTEM RELIABILITY ESTIMATE 

k 

k 

RSHTBR 


BRIDGE SYSTEM RELIABILITY ESTIMATE 

k 

k 

SEED 


PARAMETER 

k 

k 

SELCTA 


SIGNIFICANCE LEVEL SELECTION 

k 

k 

SELCTB 


SIGNIFICANCE LEVEL SELECTION 

k 

k 

SORT 


PARAMETER FOR ROUTINE SRND 

k 

k 

SUMNAI 


SUM OF THE PRODUCT OF NT’S ANT) Al’S 

k 

k 

TEMP 


TEMPORARY ARRAY 

k 

k 

TOTREP 


TOTAL NUMBER OF PROGRAM ITERATIONS 

k 

k 

TRANBR 


TEMPORARY ARRAY 

k 

k 

TRANSQ 


TEMPORARY ARRAY 

k 

k 

TRANSR 


TEMPORARY ARRAY 

k 

k 

TRIALS 


BERNOULLI TRIALS ARRAY (2-DIM) 

k 

k 

TRNSTR 


TEMPORARY ARRAY 

k 

k 

TRUQNT 


TRUE QUANTILE 

k 

k 

UNIRV 


UNIFORM RANDOM DEVIATES (2-DIM) 

k 

k 

ZFAILS 


TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES 

k 

k 

TfTfTwt 

ZFPREP 

rkkkkkkkkkkkkk^ 


NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION 

k 

k 

kk 


PARAMETER (KK=10,MAXALF=2,NPRNI^O) 

PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=. 000001) 

REAL*4 UNIRVC15,1000),TEMP(1000),QI(KK),AI(KK),AHATI(KK) 

REAL*4 QHATI(KK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFACMAXALF) 
REAL*4 DF(5),AALFA(5),SUMNAI,RSHAT(MAXALF,MAXREP),RS 
REAL*4 KEYICMAXREP),KEY2(MAXREP),KEY3(MAXREP),TRNSTR(MAXREP) 
REAL*4 DEGFRCMAXREP), QHTUPR(MAXALF,MAXREP),CHISQ(MAXALF,MAXREP) 
REAL*4 QUPAl(MAXREP), QUPA2CMAXREP),RHTSTR(MAXALF,MAXREP) 

REAL*4 DELTAR(MAXALF), TRANSQ(HAXREP),TRANSR(MAXREP),DIFF(MAXREP) 
REAL*4 DELSTRCMAXALF),NIMIN,NIMAX,NIREAL(KK) 

REAL*4 RSHTBRCMAXALF,MAXPJEP),DELBRG(MAXALF),KEY4(MAXREP) 

REAL*4 TRANBR(MAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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REAL*4 AVGN,sue,STUD 


INTEGER SEED, MULT, SORT, TRIALS(15,1000), BIGF, FlCKiC), NCKiC) 
INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTRE? 
INTEGER C1C15, REPSHD, SELCTB, ALF, FLAG, LOOP,PRNT 
INTEGER QUANTL(MAXALF), TRUQNT(MAXALF),ZFAILS,FAILS,INC 
INTEGER NTEST,FCT,BFLAG 

CHARACTER*8 LOOPSO(MAXREP) 

DATA SEED/123457/, MULT/l/, INC/i/ 

DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 

DATA ALFA/. 20,. 050/ 

DATA- SORT/0/ 

ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 


* CALL CO.MPRS 
PR-NT = NPR.NT 


DO 12 1=1,KK 

AI(I) = 9999. 
NCI) = 99999999 
12 CONTINUE 

READ(03,*)K,MSTRQ 

DO 11 1=1,K 
READC03,*) AI(I),N(I) 
11 CONTINUE 


IF(K.NE.5) THEN 

WRITECl,’("WARNING: BRIDGE STRUCTURE ", 

+"ONLy USES THE FIRST 5 COMPOiNENTS" )' ) 

ELSE 
ENT) IF 

***INTTIALIZE THE QHTUPR ARRAY OF UNHELIABILITY REPLICATIONS, *** 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 

* A.NT) RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 

***ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NTJMBER *** 

DO 172 ALF=1,HAXALF 

DO 173 REPS=1,MAXREP 

QHTUPRCALF,REPS) = 0. 

RSHATCALF,REPS) = 0. 

RHTSTRC ALF, PEPS) = 0. 

RSHTBRCALF,REPS) = 0. 

LOOPSO(REPS)=' **** ** ' 

173 CONTINUE 

172 CONTINUE 

***SET FLAG TO 1 IF ALL COMPONENTS HA\E SAME NO. OF MISSION TESTS**** 
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FLAG=1 

DO 50 1=1,K -1 

IF((N(I) T N(I+l)).NE.O) THEN 
FLAG=0 

ELSE 
END TF 
50 CONTINUE 

PRINT 'flag IS:', FLAG 

***MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)*** 

ZFPREP = 0. 

ZFAILS = 0 
FAILS =0 
TOTREP = 0 
LOOP = 0 

10 IFCLOOP. LT. MAXREP) THEN 
LOOP = LOOP + 1 
IF(TOTREP.LT.MAXRUN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT’'^** 
***0F THE SUBROUTINE SHSORT iVs'c* 

DO 95 REPS=1, MAXREP 
KEYl(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
95 CONTINUE 

***CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE*** 
***CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO.*** 

CALL IMAX(N,K,NMAX,NINDX) 


***CALCULATE THE Ql'S FROM THE GIVEN MASTER Q AND THE AI'S*** 

DO 115 1=1, K 

QI(I) = MSTRQ * AI(I) 

115 CONTINUE 

DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINUE 

*** DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS*** 
DO 130 1=1, K 
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CALL SRND(SEED, TEMP, N(r), MULT, SORT) 
DO 135 J=l, N(I:) 

UNIRV(I,J) = TEMP(J) 

IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 
TRIALS(I,JO “ 0 

ELSE 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
130 CONTINTJE 


***CALCULATE THE NO. OF FAILURES FOR EACH COMPONENT*** 

DO 150 1=1, K 
FI(I) = 0 
150 CONTINUE 
lONECT = 0 

***CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES*** 
BIGF = 0 
DO 155 1=1, K 

DO 160 J=l, N(I) 

FI(I) = FI(I) + TRIALS(I,J) 

160 CONTINUE 

IF(FI(I).EQ. 0) THEN 

ZFPREP = ZFPREP + 1 

ELSE 
END IF 

***CALCULATE THE QHAT SUB I'S: F SUB I'S DIVIDED BY N SUB I'S*** 
QHATI(I) = REALCFKD) / N(I) 

BIGF = BIGF + FI(I) 

155 CONTINUE 

***COUNTS NUMBER OF COMPONENTS THAT HAVE FAILED*** 

DO 156 1=1,K 

IF (FI(I) .NE. 0) I0NECT=I0NECT+1 

156 CONTINUE 

CALL CPARE(FI,K,BFLAG) 


***CASE WHERE NO COMPONENTS HAVE ANY FAILURES*** 

AVGN=0. 0 
DO 200 1=1,K 

AVGN=AVGN+REAL(N(I)) 

200 COOTINUE 

AVGN=AVGN/REAL(K) 

IF(BIGF.EQ. 0) THEN 

LOOPSO(LOOP)=' *ZERO* ’ 

ZFAILS = ZFAILS + 1 
DO 205 ALF=1, MAXALF 

RSHAT(ALF,LOOP)= ALFA(ALF)**(1. /AVGN) 
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lECFLAG.EQ. 1) THEN 

RHTSTR(ALF,L00P)=ALFA(ALF)**( l../N( 1)) 

ELSE 
END IF 

IF (BFLAG .EQ. 0) THEN 

RSHTBRC ALF i LOOP )=ALFA( ALF) **( 1.-/AVGN) 

ENDIF 

205 CONTINUE 

DEGFR(LOOP) = 2. 

GO TO 10 

ELSE 

' FAILS = FAILS + 1 
END IF 

***COUNTS NUMBER OF COMPONENTS THAT FAIL RECORDS NO. COMPT TESTS*** 
/FCT=0 

DO 202 1=1,K 

IF (FI(I) .NT:. 0) THEN 
FCT=FCT+1 
NTEST=N(I) 

ENDIF 

202 CONTINUE 

***FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES*** 

CALL RMAXCQHATI, K, QHTMAX, QINDX) 

***CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES)*** 

***IF COMPON'ENT HAS NO FAILURES AHAT SUB I IS ZERO*** 

SUMNAI = 0. 

DO 165 1=1, K 

AHATI(I) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + N(I) * AHATI(I) 

165 CONTINUE 

***1 COMPONENT FAILURE SERIES SYSTEM*** 

IF (FCT .EQ. 1) THEN 

LOOPSO(LOOP)=' *ONECF* ' 

DO 305 ALF=1,MAXALF 
SUC=REAL(NTEST-BIGF) 

STUD=FIN(1.-ALFA(ALF),2.*(REAL(BIGF)+1. ),2.*SUC) 

RSHATC ALF,LOOP)=SUC/(SUC+(REAL(BIGF)+1. )*STUD) 

IF (FLAG .EQ.1) THEN 

RHTSTR(ALF,LOOP)=RSHAT(ALF,LOOP) 

ELSE 

ENDIF 

305 CONTINUE 
ENDIF 

***CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY*** 


DEGFR(LOOP) = 2 * (1 + BIGF) 
DO 170 ALF=1, MAXALF 







CALL MDCHId - ALFA(ALF) ,DEGER(LOOP) ,CHISQ(ALF,LOOP) ERROR) 
QHTUPR(ALF-,LOOP) = CHISQ(ALF,LOOP;) / (-2 * SUMNAI) 

IF(FLAG. EQ. 1) THEN 

RHTSTR(ALF,LOOP) = 1 -(CHISQ(ALF,LOOP) / REAL(2*N(1))) 

ELSE- 
END IF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES*** 
IF (FCT .NE. 1) THEN 

CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF,LOOP)) 
ENDIF 

^CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE*** 

***IF NO SYSTEM FAILURE AND BRIDGE SYSTEM*** 

IF (BFLAG .EQ. 0) THEN 

RSHTBR(ALF,L00P)=ALFACALF)**C1./AVGN) 

ENDIF 

***IF MORE THAN 1 SYSTEM FAILURE AND BRIDGE SYSTEM*** 

IF OBFLAG .EQ. 2) THEN 

CALL RHTBRG(QHTUPR(ALF,LOOP),AHATI,K,RSHTBRCALF,LOOP)) 
ENDIF 

**’'fEXACTLY I SYSTEM FAILURE AND BRIDGE SYSTEM*** 

IF (BFLAG .EQ. 1) THEN 
SUC=REAL(AVGN-1) 

STUD=FIN(1, -ALFA(ALF),2-*2. ,2-*SUC) 
RSHTBR(ALF,L00P)=SUC/(SUe+2.*STUD) 

ENDIF 

170 CONTINUE 

***THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN*** 

ELSE 

WRITEC1,'( " " ,/ " PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 

+ '' ALLOWED OF: ” ,I6)‘) TOTREP 
GOTO 9999 

END IF 
GOTO 10 
END IF 

C WRITE(2,'( "UNSORTED RSHAT 1 IS: " ,/10(F8. 5)) ' ) 

C +(RSHAT(1,L00P), L00P=1, MAXREP) 

C WRITE(2, ("UNSORTED RSHAT 2 IS: " ,/10(F8. 5)) ' ) 

C +(RSHAT(2,L00P), L00P=1, MAXREP) 

C IF(FLAG.EQ. 1) THEN 

C WRITE(2,’( "UNSORTED RHTSTR 1 IS: " ,/10(F8. 5))') 

C +(RHTSTR(1,L00P). L00P=1, MAXREP) 

C WRITE(2,'(‘'UNSORTED RHTSTR 2 IS: " ,/10(F8.5))') 

C +(RHTSTR(2,L00P), L00P=1, MAXREP) 

C ELSE 
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END IF 

IFCK.-EQ. 5) THEN 

WRITE(2,'( "UNSORTED RSHTBR 1 IS: " ,/10(F8. 5)) ') 
+(RSHTBR(1,L00P). L00P=:i, MAXREP) 

WRITE(2,'(''UNSORTED RSHTBR 2 IS:" ,/10CF8.5))') 
+(RSHTBR(2,L00P), LOOP=l, MAXREP) 

ELSE 
END IF 

WRITE (2,'( " ZERO AND ONE FAILURE REPS:'',/10(A8))') 

+ (LOOPSO(LOOP),LOOP=l,MAXREP): 

*a*soRT the ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL)*** 

DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPR(ALF.REPS) 

TRANSR(REPS) = RSHAT(ALF,REPS) 

TRNSTR(REPS) = RHTSTR(ALF,REPS) 

TRANBR(REPS) = RSHTBR(ALF,REPS) 

800 CONTINUE 

CALL SHS0RT(TRANSQ,KEY1,MAXREP) 

CALL SHSORT(TRANSR,KEY2,MAXREP) 

CALL SHSORT(TRNSTR,KEYS,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 

RSHAT(ALF,REPS) = TRANSR(REPS) 

RHTSTR(ALF,REPS) ='TRNSTR(REPS) 

RSHTBR(ALF,REPS) = TRANBR(REPS) 

900 CONTINUE 

700 CONTINUE 

***PRINT OUTPUT REPORT HEADINGS*** 

WRITE(1,6666) 

WRITE(1,6667) MAXREP 
WRITE(1,6668) K 
WRITE(1,6669) 

IF(K.EQ. 5) THEN 
WRITEC1,6699) 

ELSE 
END IF 

WRITEC1,6670) MSTRQ 
WRITEC1,6671) 

WRITEC1,C1C15) 

WRITEC1,3334) AI 
WRITEC1,0007) 

WRITEC1,C1C15) 

WRITEC1,3334) QI 
WRITEC1,0005) 

WRITEC1,C1C15) 

WRITEC1,3335) N 
WRITEC1,6674) 

***COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. CSERIES SYSTEM)*** 
*** AND FOR THE 5-COMPONENT BRIDGE STRUCTURE*** 
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CALL RSRS(QI,K,RS) 

WRITECl , ’ ( " ' ',///-/' 'THE TRUE SERIES SYSTEM ' ' , 

+ " RELIABILITy VALUE IS: ” ,T51,F8. 5) ’ ) RS 
CALL RBRIDG(QI,K,RSBRDG) 

IE(K.EQ.5) THEN 

WRITECl,'( " ",////''THE TRUE BRIDGE STRUCTURE '', 

+ "RELIABILITY VALUE IS:' ',T51,F8. 5) ' ) RSBRDG 
ELSE 
END IF 

WRITEC1,6675) 

***COMPUTE THE DIFFERENCE 'DELTAR' BTWN. RS AND RSHAT OF THE THEO*** 
***RETICAL quantile given by ALFACMUST USE SORITD RSHAT ARRAY)*** 

IFCFLAG. EQ. 1) THEN 
WRITEC1,5755) 

ELSE 
END IF 

DO 450 ALF=1, MAXALF 

QUANTLCALF) = MAXREP * Cl - ALFACALF)) 

DELTARCALF) = RS - RSHATCALF,QUANTLCALF)) 

DELBRGCALF) = RSBRDG - RSHTBRCALF,QUANTLCALF)) 

IFCFLAG. EQ. 1) THEN 

DELSTRCALF) = RS - RHTSTRCALF,QUANTLCALF)) 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5656) RHTSTRCALF,QUANTLCALF)) 

WRITEC1,5657) DELSTRCALF) 

ELSE 
END IF 

IFCK.EQ.5) THEN 

DELBRGCALF) = RSBRDG - RSHTBRCALF,QUANTLCALF)) 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5666) RSHTBRCALF,QUANTLCALF)) 

WRITEC1,5667) DELBRGCALF) 

ELSE 
END IF 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 

WRITEC1,5556) RSHATCALF,QUANTLCALF)) 

WRITEC1,5557) DELTARCALF) 

450 CONTINUE 

PRINT*, 'QUANTLC 1) IS:', QUANTLC 1) 

PRINT*, 'QUANTLC2) IS:', QUANTLC2) 

***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

**>v ******** RSHAT ******* *** 

WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - xTSHATCALF,REPS) 

500 CONTINTJE 

DO 600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 


37 




WRITE(1,',( ' ' " ,/' 'TRUE CONFIDENCE LIMIT IS: " , 

+ E8.4)') 

+ (TRUQNTCALF) / REAL(MAXREP)) * 100. 

GO TO 620 

ELSEIF(DIFF(REPS).LT.0.) THEN 
TRUQNTCALF) = REPS 
GO TO 610 

ELSE 
END IF 

600 CONTINUE 

610 IF(TRUQNTCALF). EQ. 0. ) THEN 

WRITEC1,4443) ALFA(ALF) 

WRITE(1,'(" ",/"THE SMALLEST’' , 

+ '' DIFFERENCE BETWEEN RS AND RSHAT IS: " ,F10. 5)' ) DIFF( 

+ MAXREP) 

ELSEIF(TRUQNT(ALF).EQ. 1. ) THEN 
WRITEC 1,4442:) ALFACALF) 

WRITEC1,'C" ",/"all RSHAT" , 

+ ' ' ARE GREATER THAN RS" ) ' ) 

ELSEIFCABSCDIFFCTRUQNTCALF'')).LE.ABSCDIFFCTRUQNTCALF) - 1))) 
+ THEN 

WRITEC1,4444) ALFACALF), 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)) 

WRITEC 1,4446;) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)-!) 

WRITEC1,4447) 

620 END IF 

400 CONTINUE 


***FiND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** *5 V**a*5V5V RSHTBR C BRIDGE) ******* *** 


IFCK.EQ.5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 501 REPS=1, MAXREP 

DIFFCREPS) = RSBRDG - RSHTBRCALF,REPS) 

501 CONTINUE 

DO 601 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C" ",/"TRUE CONFIDENCE LIMIT IS:", 
+ F8.4)') 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 621 

ELSEIFCDIFFCREPS).LT. 0. ) THEN 
TRUQNTCALF) = REPS 
GO TO 611 

ELSE 
END IF 

601 CONTINTJE 
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611 IF(TRUQNT(ALF).EQ.0. ) THEN 

WRITEC1,4443) ALFA(ALF) 

WRIT£(" ,/''THE SMALLEST'', 

+ '' DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS:'', 

+ FIO.50’) DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ.l. ) THEN 
TOITE(1,4442) ALFA(ALF) 

WRITEC 1,' ( " ", /'' ALL RSHTBR ’ ' , 

+ "ARE GREATER THAN RSBRDG'’)') 

ELSEIF(ABS(DIFFtTRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITEC 1,4444) ALFACALF);, 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF.TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

621 END IF 

401 CONTINUE 
ELSE 
END IF 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

*** ***5V5V*Vr* RHTSTR ****5V**;V 

IFCFLAG, EQ. 1). THEN 
DO 4400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 5500 REPS=1, MAXREP 

DIFFCREPS) = RS - RHTSTRCALF,REPS) 

5500 CONTINUE 

DO 6600 REPS=1, MAX'REP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC 1, ’ C " " ,/"TRUE CONFIDENCE LIMIT IS: " , 
+ F8.4)') 

+ CTRUQ.NTCALF) / REALCMAXREP)) * 100. 

GO TO 6620 

ELSEIFCDIFFCREPS).LT. 0, ) THEN 
TRUQNTCALF) = REPS 
GO TO 6610 

ELSE 
END IF 

6600 CONTINUE 

6610 IFCTRUQNTCALF).EQ. 0.) THEN 
TOITEC1,4443) ALFACALF) 

WRITEC1,'C" ",/"THE SMALLEST", 

+ ’' DIFFERENCE BETWEEN RS AND RHTSTR IS: " , 

+ F9.5)') DIFFCMAXREP) 

ELSEIFCTRUONTCALF). EQ. 1.) THEN 
WRITEC1,4442) ALFACALF) 

WRITEClj'C" ALL RHTSTR" , 
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+ "ARE GREATER THAN RS" ) ' ) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE. ABSCDIFF(TRUQNT(ALF) - 1) ) ) 
+ THEN 

TOITEC1,44^4) ALFA(ALF), 

+ (TRUQNT(ALF) / REAL(MAaREP)) * 100. 

WRITE(T,4448) RHTSTR(AEF,TRUQNT(ALF)) 

WRITE(1:,4446) 

ELSE 

WRITECl,4444) ALFA(ALF), 

+ ((TRUQNT(ALF)-l) / REAL(MAXREP)) * 100. 

WRITE(1,4448) RHTSTR(ALF,TRUQNTC ALF)-1) 

WRITE(1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION*** 

IF(PRNT. EQ. 1) THEN 
1 = 1 

185 WRITEC1,REPSHD) ALFA(SELCTA), ALFA(SELCTA), 

+ALFA(SELCTB),ALFA(SELCTB),ALFA(SELCTA),ALFA(SELCTA),ALFACSELCTB)., 
+ALFA(SELCTB) 

175 IFCI.GE. (MAXREP + 1)) THEN 
GOTO 180 

ELSE 

IF( (I. EQ. 71). OR. (I. EQ. 211). OR. (I. EQ. 351). OR. (I. EQ. 491). OR. 

+ (I.EQ.631).0R. (I.EQ.771).OR.(I.EQ.9ll).OR.(I.EQ. 1051) ) THEN 

I = I + 70 
WRITECl,'(" + ")') 

GOTO 185 

ELSE 

WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,I), QHTUPR(1,I), 

+ CHISQ(2,I), QHTUPR(2,I) 

END IF 

IF(I + 70. LE. MAXREP) THEN 

WRITEC1,3337) 1+70,INT(DEGFR(1+70)),CHISQ(1,1+70), 

+ QHTUPRC1,1+70),CHISQC 2,1+70),QHTUPR(2,1+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENTIIF 

9999 WRITECl,'("THE TOTAL NO OF REPS WAS:", 18)') TOTREP 

WRITECl,'("THE TOTAL NO OF EFFECTIVE REPS WAS: ",I8)') LOOP 
WRITECl,'("THE TOTAL NO OF NO FAILURE RUNS WAS; ",I8)') ZFAILS 
WRITECl,'("AVERAGE NO. OF COMPONENTS PER REPLICATION WITH ", 

+''NO FAILURES:'',F5.2)') ZFPREP / MAXREP 
WRITECl,’("THE TOTAL NO OF RUNS WITH FAILURES WAS: ",I8)') FAILS 
0008 FOR.MAT (/ 3X,'C l',5X,'C 2', 

+5X,'C 3',5X,’C 4',5X,'C 5',5X,'C 6’,5X,'C 7',5X, 

+'C 8',5X,'C 9',5X,'C 10',4X,'C ll',4X, 
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+'C 12',4X,'C 13',4X,'C 14'.4X,'C 15’) 

0009 F0RHAT(/1X,’REP NO’,2X.’DF‘,IX ’CHISQR(’,F4. 3,’)’,IX, 

+ 'QKTUPR(’,F4. 3,’)',IX,'CKISQR(‘,F4. 3,')’,IX,'QKTUPR(',F4.3,’)', 
+2X,’REP NO' ,2X,'DF',1X,’CHISQR(',F4. 3,')\1X, 

+’QHTUPR(’,F4. 3,')',IX,'CHISQR(’,F4. 3,’)’,IX,'QHTUPR(’,F4. 3,’)‘/) 
0001 FORMAT (///’UNIFORM RANDOM DEVIATES ARE:') 

0002 FORMAT (///’BERNOULLI TRIALS ARE:’) 

0003 FORMAT (///’TOTAL NO. OF FAILURES FOR EACH COMPONENT: ’) 

0004 FORMAT (///’ESTIMATED UNRELIABILITY FOR EACH COMPONENT:') 

0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS: ’) 

0006 FORMAT (///’ESTIMATED WEIGHTS FOR EACH COMPON'ENT: ') 

0007 FORMAT (///’Q I FOR EACH COMPONENT: ’) 

1111 FORMAT (15F8.5) 


2222 FORMAT (/IX,15(14,4X)) 

3333 FORMAT (/IX,15(I4,4X)) 

3334 FORMAT (/15F8.5) 

3335 FORMAT (/IX,15(14,4X)) 

3336 FORMAT (T3,I4,T9,I3,T13,F11. 5,T27,F8. 5,T39,F11.5,T53,F8. 5) 

3337 FORMAT (’ + ',T67,I4,T73,I3,T77,F11. 5,T91,F8. 5,T103,F11. 5,T117,F8.5) 

4442 FORMAT (’ ',///'THE RESULTING (1 - ',F4.3,’) CONFIDENCE ’, 

+'LIMIT IS:'.T50,' 00.000 ’) 

4443 FORMAT (' ’,///'THE RESULTING (1 - ',F4.3,’) CONFIDENCE ’, 

+’LIMIT IS: ',T50,’100. 0000’) 

4444 FORMAT (’ ',///'tKE RESULTING (1 - ’,F4.3,’) CONTIDENCE ', 

+’LIMIT IS:’,T50,F8.4) 

4445 FORMAT (’ ’,/'THE RSHAT VALUE CLOSEST TO RS IS: ’,T51,F8.5) 

4446 FORMAT (’ ’,/’(FIRST NEGATIVE DIFFERENCE)’) 

4447 FORMAT (’ ',/'(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)’) 

4448 FORMAT (' ’,/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4449 FORMAT (' ',/'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T51,F8.5) 

5555 FORM.AT (’ ',///'THE ' , 14, ’ (1-',F4. 3,' ) QUANTILE IS: ’ ,T49 ,F8. 3) 

5556 FOR.MAT (' ',/’THE VALUE OF RSHAT FOR THAT QUANTILE IS: ’,T51,F8.5) 

5557 FORMAT (’ ’,/'THE DIFFERENCE(RS - RSHAT) IS:’,T51,F8.5) 

5656 FORMAT (’ ’,/'THE VALUE OF RHTSTR FOR THAT QUANTILE IS: ',T51,F8.5) 

5666 FORMAT (’ ',/’THE VALUE OF RSHTBR FOR THAT QUANTILE IS:’,T51,F8.5) 

5657 FORM.AT (’ ',/'THE DIFFERENCE(RS - RHTSTR) IS: ’ ,T51 ,F8. 5) 

5667 FORMAT (’ ’,/’THE DIFFERENCE(RS - RSHTBR) IS:’,T51,F8. 5) 

5755 FORMAT (' ',///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR’, 

+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY ’, 

+'’’RHTSTR'' IS COMPUTED’) 

6666 FORMAT ('+’,'* 

+' **?V******** pujv} INPUX SETTINGS ■>WrftTr*A***********A**********5VA*' ^ 

6667 FORMAT (' ’,//’NUMBER OF REPLICATIONS: ’,T50,14) 

6668 FORMAT (' ’,//’N'UMBER OF COMPONENTS: ' ,T50,14) 

6669 FORMAT (' './/'SYSTEM RELIABILITY FUNCTION: ’,T50,’SERIES’) 

6699 FORMAT (’ ’.//’SYSTEM RELIABILITY FUNCTIO.N: ' ,T50,’BRIDGE') 

6670 FORMAT (' './/'MASTER UNRELIABILITY USED: ',T50,F8.5) 

6671 FORMAT (' './/'INPUT WEIGHTS(A SUB I"S): ’) 

6674 FORMAT (' + ' ///’ ****VbV*A****;V********Vn;rfr*A*****iV******5Wt***rir***' ^ 

+’**R UN r’eSULT S*****’'^*********'^***************— 

+' A-*************' J 

+ ' ESTIMATE ERRORS i^*5r**A5V*A;VjV*A»V*******T:jV**Tr*Tr**T:**«?V*A***iV***’ ^ 


4448 FORMAT (' ’ 

4449 FORMAT (' ' 
5555 FORM.AT (’ ' 
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6676 FORMAT ('■!■' J / / / ' ********:V*****iV5VA********5V*iVi'r*******T:***;V****5V* ' 
+' TRUE CONFIDENCE LIMITS **********’'**************************** 

+' * * 5'r5VA*******5V*' ^ 


? 


END 
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APPENDIX D. FORTRAN CODE FOR THE PREFERRED LOWER 
CONFIDENCE LIMIT FOR SYSTEM RELIABILITY (SERIES PARALLEL 

SYSTEM) 


PROGRAM ZFYSCN 


A ^ ^ ^ j * y * 

* 


* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE 

* ZERO FAILURES ALLOWED; NO SCALING 

* AUTHOR: E. F. BELLINI, LT, USN 

* MODIFIED BY: LT VALERIE A. COVINGTON,USN (MAR 90) 

* DATE: NOV 89 


THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 
RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 
OF THEIR COMPONENTS * 


* IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 

* TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT * 

* SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING ’El EXEC'.* 

* THE REXX EXEC PROGRAM * 

* 'Bl' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT * 

* FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITI.NG THE 

* LNDEX COUNTERS I, J, K OF THE 'Bl' EXEC ONE CAN RUN A.NY USER- * 

* SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. 

it it 

* VARIABLES USED * 


it 

AHATI 

it 

AI 

it 

ALFA 

it 

BIGF 

it 

CHISQ 

it 

C1C15 

it 

DEGFR 

it 

DELBRG 

it 

DELSTR 

it 

DELTAR 

it 

DIFF 

it 

EPS 

it 

ERROR 

it 

FAILS 

it 

FI 

it 

FLAG 

it 

INC 

* 

KEYl 

it 

KEY2 

it 

KEYS 

it 

KEY4 

it 

KfC 

it 

LOOP 

it 

MAXALF 


NFIGHT ESTIMATES FOR EACH COMPONF.NT * 
INPUT WEIGHTS FOR EACH COMPONENT * 
LEVELS OF SIGNIFICANCE * 
TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHI-SQUARE RANDOM VARIABLE VALUE * 
FORMAT LABEL * 
DEGREES OF FREEDOM * 
DIFFERENCE FOR BRIDGE SYSTEM * 
DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DIFFERENCE FOR SERIES SYSTEM * 
DIFFERENCE (TRUE REL. - ESTIMATED REL. ) * 
SMALL QUANTITY(CONSTANT) * 
PARAMETER FOR IMSL ROUTINE * 
COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
NO. OF FAILU-RES FOR EACH COMPON'ENT(ALL MISSION TST)* 
1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INCREMENT STEP SIZE FOR ROUTINF USMNILX * 
ARRAY OF INDE^,ES FOR ROUTI.NE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONDNTS* 
COUNTS NO. OF REPLICATION PERFOR.MED * 
MAX NO. OF SIGNIFICANT LEVELS DESIRED(ARRAy SIZING)* 






* 

MAXREP 


MAX NO. OF REPLICATIONS 

•k 

* 

MAXRUN 


MAX NO. OF PROGRAM ITERATIONS ALLOWED 

k 

* 

MSTRQ 


MASTER UNRELIABILITY: USED WITH Al 'S TO CALC. Ql'S) 

k 

* 

MULT 


MULTIPLIER FOR RANDOM NO. GENERATOR SR.ND 

* 

* 

N 


NO. OF mission; test FOR EACH COMPONENT * 

* 

NIMAX 


MAX NO. OF MISSION TESTS 

* 

* 

NIMIM 


MIN NO. OF MISSION TESTS 

* 

* 

NINDX 


INDEX NO. OF MAX NO. OF MISSION TESTS 

* 

* 

NIREAL 


NO. OF MISSION TESTS TRANSFORMED TO REAL 

k 

* 

NMAX- 


MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 

k 

* 

NPRNT 


FLAG FOR DETAILED REPORT OUTPUT 

k 

* 

PRNT 


SAME AS ABOVECPARAMETER) 

k 

* 

QHATI 


UNRELIABILITY ESTIMATES FOR EACH COMPONENT 

* 

* 

QHTMAX 


LARGEST QHATI 

k 

* 

QHTUPR 


UPPER LIMIT ON SYSTEM UNRELIABILITY 

k 

* 

QI 


INPUT UNRELIABILIY FOR EACH COMPONENT 

k 

* 

QINDX 


INDEX 

k 

* 

QUANTL 


QUANTILE 

k 

* 

REPSHD 


REPLICATIONS HEADING FORMAT NUMBER 

k 

* 

RHTSTR 


SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 

k 

* 

RS 


TRUE SERIES SYSTEM RELIABILITY 

k 

* 

RSBRDG 


TRUE BRIDGE SYSTEM RELIABILITY 

k 

* 

RSHAT 


SERIES SYSTEM RELIABILITY ESTIMATE 

k 

* 

RSHTBR 


BRIDGE SYSTEM RELIABILITY ESTIMATE 

k 

is 

SEED 


PARAMETER 

k 

* 

SELCTA 


SIGNIFICANCE LEVEL SELECTION 

k 

* 

SELCTB 


SIGNIFICANCE LEVEL SELECTION 

k 

is 

SORT 


PARAMETER FOR ROUTINE SRND 

k 

* 

SUMNAI 


SUM OF THE PRODUCT OF NT'S ANT) Al'S 

k 

* 

TEMP 


TEMPORARY ARRAY 

k 

* 

TOTREP 


TOTAL NUMBER OF PROGRAM ITERATIONS 

k 

is 

TRANBR 


TEMPORARY ARRAY 

k 

* 

TRANSQ 


TEMPORARY ARRAY 

k 

* 

TRANSR 


TEMPORARY ARRAY 

k 

is 

TRIALS 


BERNOULLI TRIALS ARRAY ,(2-DIM) 

k 

is 

TRNSTR 


TEMPORARY ARRAY 

k 

is 

TRUQNT 


TRUE QUANTILE 

k 

is 

UNIRV 


UNIFORM RANDOM DEVIATES (2-DIM) 

k 


ZFAILS 


TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES 

k 

is 

ZFPREP 


NO. OF COMPNTS. WITH ZERO FAILUTIES PER REPLICATION 

k 


* 




PARAMETER (paK=10,MAXALF=2,NPRNT=0) 

PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=. 000001) 

REAL*4 UNIRVC15,1000),TEMP(1000),QI(KK),AI(KK),AHATI(KK) 

REAL*4 QHATI(KK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFA(MAXALF) 
REAL*4 DF(5),AALFA(5),SUMNAI,RSHAT(MAXALF,MAXREP),RS 
REAL*4 KEYl(MAXREP),KEY2(MAXREP),KEY3(MAXREP),TRNSTR(MAXREP) 

REAL*4 DEGFRC MAXREP), QHTUPR(MAXALF,MAXTIEP),CHISQ(MAXALF,MAXREP) 
REAL*4 QUPAKMAXREP), QUPA2(MAXREP),RHTSTR(MAXALF,MAXREP) 

REAL*4 DELTARCMAXALF), TRANSQ(MAXREP),TRANSRCMAXREP),DIFF(MAXREP) 
REAL*4 DELSTRCMAXALF),NIMIM,NIMAX,NIREAL(KK) 

REAL*4 RSHTBRCMAXALF,MAXREP),DELBRG(MAXALF),KEY4(MAXREP) 

REAL*4 TRANBRCMAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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REAL*4 AVGN,SUC,STUD,PTEMP(10),FDES1,FDEG2,P3,S3,S4 

INTEGER SEED, MULT, SORT, TRTALS(15,1000), BIGF, FI(KK), N(KK) 
INTEGER NTKDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTREP 
INTEGER C1C15, REPSHD, SELCTB, ALF, FLAG, LOOP,PRNT 
INTEGER QUANTLC HAXALF), TRUQNT(MAXALF),ZFAILS,FAILS,INC 
INTEGER NTEST,FCT,HFI 

CHARACTER*8 LOOPSOCMAXREP) 

DATA SEED/123457/, MULT/1/, INC/1/ 

DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 

DATA ALFA/. 20,. 050/ 

DATA SORT/0/ 

ASSIG.N 8 TO C1C15 
ASSIGN 9 TO REPSKD 

CALL COMPRS 
PR.NT = NPR.NT 

DO 12 1=1,KK 

Aid) = 9999. 

N(l) = 99999999 
12 CONTINUE 

PvEAD(03,*)K,MSTRQ 

DO 11 1=1,K 
READ(03,*) AI(I),N(I) 

11 CONTINUE 


IF(K.NE.5) THEN 

WRITEC1,'(''WARNING: BRIDGE STRUCFuRE '', 

+"ONLY USES THE FIRST 5 COMPONENTS")') 

ELSE 
END IF 

***INTTIALIZE THE QHTUPR ARRAY OF UNRELIABILITl’ REPLICATIONS, *** 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 

* AND RhTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 

***ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NUMBER *** 

DO 172 ALF=1,MAXALF 

DO 173 REPS=1,MAX'REP 

QHTU?R(ALF,REPS) = 0. 

RSHATCALF,REPS) = 0. 

RHTSTRCALF,REPS) = 0. 

RSHTBRCALF,REPS) = 0. 

LOOPSO(REPS)=' **iHhv* • 

173 CONTINUE 

172 CONTINUE 

***SET FLAG TO 1 IF ALL COMPONENTS HA\E SAME NO. OF MISSION TESTS**** 
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FLAG=1 

DO 50 1=1,K -1 

IF((N(I) - N(I+l)).NE.O) THEN 
FLAG=0 

ELSE 
END IF 
50 CONTINUE 

PRINT *, 'FLAG IS:', FLAG 

***MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)*** 

ZFPREP = 0. 

2FAILS = 0 
FAILS = 0 
TOTRCP = 0 
LOOP = 0 

10 IF(LOOP. LT.MAXREP) THEN 
LOOP = LOOP + 1 
IF(TOTREP.LT.MAXRUN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT’^** 
***OF THE SUBROUTINE SHSORT *** 

DO 95 REPS=1, MAXREP 
KEY1(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
95 CONTINUE 

***CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE*** 
***CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO.*** 

CALL IMAX(N,K,NMAX,NINDX) 


***CALCULATE THE Ql'S FROM THE GIVEN MASTER Q AND THE Al'S*** 

DO 115 1=1, K 

QI(I) = MSTRQ * AI(I) 

115 CONTINUE 
C S3=QI(3) 

S4=QI(1)**.50 
DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINUE 

*** DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS*** 
DO 130 1=1, K 
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CALL SRNDCSEED, TEMP, N(I), MULT, SORT) 
DO 135 J=l, NCI) 

UNIRV(I,J) = TEMP(J) 

IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 
TRIALS(I,J) = 0 

VT 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
130 CONTINUE 


***CALCULATE THE NO. OF FAILURES FOR EACH COMPONENT*** 

DO 150 1=1, K 
FI(I) = 0 
150 CONTINUE 
lONECT = 0 

***CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES*** 
BIGF = 0 
DO 155 1=1, K 

DO 160 J=l, NCI) 

FICI) = FICI) + TRIALSCI,J) 

160 CONTINUE 

IF CFICI) .EQ. 0) THEN 
ZFPREP=ZFPREP+1 
ELSE 
ENDlf 

BIGF = BIGF + FICI) 

QHATICI) = REALCFICI)) / NCI) 

155 CONTINUE 

***CHANGES FOR SERIES PARALLEL SYSTEM CCOMPONENT 2)*** 

IF CFIC2) .EQ. 0) THEN 
ZFPREP=ZFPREP - 1 
ENDIF 

BIGF = BIGF - FIC2) 

FICI) = 0 
DO 161 1=1,NC2) 

HFI=0 

DO 162 J=l,2 

CALL SRNDCSEED, PTEMPCJ), NC2), MULT, SORT) 

IF CPTEMPCJ) .GT. 1-S4) THEN 
HFI = HFI + 1 
E.NDIF 
162 CONTINUE 

IF CHFI .EQ. 2) THEN 
FIC2) = FIC2) + 1 
ENDIF 

161 CONTINUE 

BIGF = BIGF + FIC2) 

IF CFIC2) .EQ. 0) THEN 
ZFPREP = ZFPREP + 1 
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MULT, SORT) 


ELSE 
EKDIF 

C IF (FI(3) .EQ. 0) THEN 

-C ZFPREP=ZFPREP - 1 

C ENDIF 

C BIGF = BIGF “ FI(3) 

-C FI(3) = 0 

C DO 163 I=1,N(3) 

C HFI=0 

C °° CALL'^sk’DCSEED, PTEMP(J). N(3) 

C IF (PTEMP(J) .GT. 1-S3) THEN 

G HFI = HFI + 1 

C ENDIF 

C164 CONTINUE 
C IF (HFI .GE. 2) THEN 

C FI(3) = FI(3) + 1 

C ENDIF 

C163 CONTINUE 
C BIGF = BIGF + Fl(3) 

C IF (FI(3) .EQ. 0) THEN 

C ZFPREP = ZFPREP + 1 

C ELSE 

C ENDIF 

DO 19 1=1,K 

QI(I) = MSTRQ * Aid) , ^ 

QHATICD = REALCFKD) / N(I) 

19 CONTINUE 
C QI(4) = QI(4)**2 

r OT('3lVi^^((3.*P3**2*QI(3))+(P3**3)) 

L.COUCTS NWIBER OF COMPONENTS THAT HAVE FAILED**" 


IF^CFia) -NE. 0) IONECT=IONECT+I 
156 CONTINUE 

***CASE WHERE NO COMPONENTS HAVE ANY FAILURES*’^* 


200 


* 

205 


IF(BIGF.EQ.O) THEN 

LOOPSO(LOOP)=’ *ZERO” 
ZFAILS = ZFAILS + 1 


1 


AVGN=0.0 
DO 200 1=1,K 

AVGN=AVGN+REAL(N(I)) 


CONTINUE 

AVGN=AVGN/REAL(K) 

DO 205 ALF=1, MAXALF 
RSHAT(ALF,LOOP)= 
IFCFLAG.EQ. 1) THEN 

RHTSTR(ALF,LOOP) 


ALFA(ALF)**(1. /AVGN) 
=ALFA(ALF)**C1./N(l)) 


ELSE 
END IF 

PRINT *,LOOP 
CONTINUE 


, LOOP,’RSHAT’,RSHATC ALF,LOOP) 
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DEGFR(LOOP) = 2. 

GO TO 10 

FT FF 

FAILS = FAILS + 1 
EKD IF 

***COUNTS NUMBER OF COMPONENTS THAT FAIL RECORDS NO, COMPT TESTS*** 
FCT=0 

DO 202 1=1,K 

IF (FI(I) .NE. 0) THEN 
FCT=FCT+1 
NTEST=N(I) 

ENTIIF 

202 CONTINUE 

***FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES*** 

CALL RMAXCQHATI, K, QHTMAX, QINDX) 

IF (LOOP .EQ. 1) THEN 
ENT)IF 

***CALCULATE THE AHAT SUB I’S (KTIGHT ESTIMATES;)*** 

***IF COMPONENT HAS NO FAILURES AHAT SUB I IS ZERO*** 

SUMNAI = 0. 

DO 165 1=1, K 

AHATI(I) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + N(I) * AHATI(I) 

165 CONTINUE 

***1 COMPONENT FAILURE SERIES SYSTEM*** 

IF (FCT .EQ. 1) THEN 

LOOPSO(LOOP)=' *ONECF* ' 

DO 305 ALF=1,MAXALF 
SUC=REAL(NTEST-BIGF) 

FDEG1=2.*(REAL(BIGF)+1, ) 

FDEG2=2. *SUC 

STUD=FIN(1. -ALFA(ALF),FDEG1,FDEG2) 
RSHAT(ALF,LOOP)=SUC/(SUC+(REAL(BIGF)+l. )*STUD) 

* PRINT *,'SUC=',SUC,'FAIL=’,BIGF,’NTEST=’,NTEST 

* PRINT *,'FIN=’.STUD 

* PRINT *,'ALFA=',ALFA(ALF),'RSHAT=',RSHAT(ALF,LOOP) 

IF (FLAG .EQ.1) THEN 

RHTSTR(ALF,LOOP)=RSHAT(ALF,LOOP) 

ELSE 

ENDIF 

305 CONTINUE 
ENDIF 

***CALCULATE 1 REPLICATION OF UPPR ALFA C. L. ON SYSTEM RELIABILITY*** 
DEGFR(LOOP) = 2 * (1 + BIGF) 

DO 170 ALF=1, MAXALF 

CALL MDCHKl - ALFA(ALF),DEGFR(LOOP) ,CHISQ(ALF,LOOP), ERROR) 
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QHTUPR(ALF,LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAI) 

IFCTLAG. EQ. 1) THEN 

RHTSTR(ALF,LOOP) = 1 -(CHISQ(ALF,LOOP) / REAL(2*N(1))) 

ELSE 
END IF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES*** 
IF (FCT .NE. 1) THEN 

CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF,LOOP)) 
ENDIF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE*** 

IF (lONECT .NE. 1) THEN 

CALL RHTBRGC QHTUPR(ALF, LOOP)., AHATI, K,RSHTBR( ALF, LOOP)) 
ENDIF 

170 CONTINUE 

***EXACTLY 1 COMPONENT FAILS AND REDUNDANT COMPONENT*** 

IF ((lONECT .EQ. 1) .AND. (K .EQ.5)) THEN 
DO 207 1=1, K 

NIREAL(I) = REALCNCI)) 

207 CONTINUE 

CALL USMNMX(NIREAL,K,INC,NIMIN,NIMAX) 

DO 206 ALF=1,MAXALF 

RSHTBR(ALF,LOOP)=ALFA(ALF)**(1. /NIMIN) 

206 CONTINUE 
ENDIF 

***THIS ELSE AND ENT)IF ARE FOR THE TEST AGAINST MAXRUN*** 

ELSE 

WRITE(1,'( " " ,/ " PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+" ALLOWED 0F:",I6)‘) TOTREP 
GOTO 9999 

END IF 
GOTO 10 
END IF 

C WRITE(2,'( "UNSORTED RSHAT 1 IS: " ,/10(F8. 5)) ' ) 

C +(RSHAT(1,L00P), L00P=1, MAXREP) 

C WRITE(2,T"UNSORTED RSHAT 2 IS: " ,/10(F8. 5))' ) 

C +(RSHAT(2,L00P), L00P=1, MAXREP) 

C IFCFLAG.EQ. 1) THEN 

C WRITE(2,'( "UNSORTED RHTSTR 1 IS: " ,/10(F8. 5))') 

C +(RHTSTR(l,LOOP), L00P=1, MAXREP) 

C WRITE(2,'(''UNSORTED RHTSTR 2 IS: " ,/10(F8.5))') 

C +(RHTSTR(2,L00P), LOOP=l, MAXREP) 

C ELSE 

C END IF 

C IFCK.EQ.5) THEN 

C WRITE(2,'( "UNSORTED RSHTBR 1 IS: " ,/10(F8. 5))' ) 

C +(RSHTBR(l,LOOP), LOOP=l, MAXREP) 
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WRITE(2,'( "UNSORTED RSHTBR 2 IS: " ,/10(F8. 5)) ') 
+(RSHTBR(2,LOOP), LOOP=1, MAXREP) 

ELSE 
END IF 

WRITE (2,'( " ZERO AND ONE FAILURE REPS:'’,/10(A8))') 

+ (LOOPSOCLOOP),LOOP=l,MAXREP) 

***SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL)*** 

DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPR(ALF,REPS) 

TRANSR(REPS) = RSHAT(ALF,REPS) 

TRNSTR(REPS) = RHTSTR(ALF,REPS) 

TRANBR(REPS) = RSHTBR(ALF,REPS) 

800 CONTINUE 

CALL SHS0RT(TRANSQ,KEY1,MAXREP) 

CALL SHSORT(TRANSR,KEY2,MAXREP) 

CALL SHSORT(TRNSTR,KEYS,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 

RSHAT(ALF,REPS) = TRANSR(REPS) 

RHTSTR(ALF,REPS) = TRNSTR(REPS) 

RSHTBR(ALF,REPS) = TRANBR(R£PS) 

900 CONTINUE 

700 CONTINUE 

***PRINT OUTPUT REPORT HEADINGS*** 

WRITEC1,6666) 

WRITEC1,6667) MAXREP 
WRITEC1,6668) K 
raiTEC1,6669) 

IF(K. EQ. 5) THEN 
WRITEC1,6699) 

ELSE 
END IF 

WRITEC1,6670) MSTRQ 
WRITEC1,6671) 

WRITEC1,C1C15) 

WRITEC1,3334) AI 
WRITEC1,0007) 

WRITEC1,C1C15) 

WRITEC1,3334) QI 
WRITEC1,0005) 

WRITEC1,C1C15) 

WRITEC1,3335) N 
WRITEC1,6674) 

***COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL, FNCTN.CSERIES SYSTEM)*** 
*** AND FOR THE 5-COMPONENT BRIDGE STRUCTURE*** 


CALL RSRSCQI,K,RS) 
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WRITE(1,'(" ",////"THE TRUE SERIES SYSTEM ", 

+ ''RELIABILITY VALUE IS; '',T51,F8. 5)') RS 
CALL RBRIDG(QI,K,RSBRDG) 

IF(K.EQ.5) THEN 

WRITECl,'( " " ,//// ''THE TRUE BRIDGE STRUCTURE '', 

+’'RELIABILITY VALUE IS:'',T51,F8. 5)') RSBRDG 
ELSE 
END IF 

WRITE(1,6675) 

***COMPUTE THE DIFFERENCE 'DELTAR' BTWN. RS AND RSHAT OF THE TKEO*** 
**.vretICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY)*** 

IFCFLAG.EQ. 1) THEN 
WRITEC1,5755) 

ELSE 
END IF 

DO 450 ALF=1, MAXALF 

QUANTL(ALF) = MAXREP * (1 - ALFA(ALF)) 

DELTARCALF) = RS - RSHAT(ALF,QUANTL(ALF)) 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF)) 

IFCFLAG.EQ.1) THEN 

DELSTR(ALF) = RS - RHTSTR(ALF,QUANTL(ALF)) 

WRITECl,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5656) RHTSTRCALF,QUANTLC ALF)) 

WRITEC1,5657) DELSTRCALF) 

ELSE 
END IF 

IFCK.EQ.5) THEN 

DELBRGCALF) = RSBRDG - RSHTBRCALF,QUANTLCALF)) 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5 666) RSHTBRC ALF,QUANTLC ALF)) 

WRITEC1,5667) DELBRGCALF) 

ELSE 
END IF 

WRITEC 1,5555) MAXREP, ALFACALF), REALCQUA.NTLCALF)) 

WRITE C1,5 5 5 6) RSHATC ALF,QUANTLC ALF)) 

WRITEC1,5557) DELTARCALF) 

450 CONTINUE 

PRINT*, 'QUANTLC1) IS:', QUANTLC1) 

PRINT*, 'QUANTLC2) IS:', QUANTLC2) 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** ******** RSHAT ******* *** 


WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - RSHATCALF,REPS) 

500 CONTINUE 

DO 600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C " " ,/ " TRUE CONTIDENCE LIMIT IS: " , 
+ F8.4)') 
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^ GO^TO^'620^^ ^ REAL(MAXREP)) * 100. 

ELSEIF(DIFF(REPS).LT.0.) THEN 
TRUQNT(ALF) = REPS 
GO TO 610 

ELSE 

600 CONTINUE 

610 IFCTRUQNT(ALF). EQ. 0-) THEN 

WRITE(1,4443) ALFA(ALF) 

+ f'' >/' 'the smallest' ’, 

+ MAXrE) “ IS: " ,no. 5)') DIFF( 

ELSEIF(TRUQNT(ALF). EQ. 1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,'C’' ",/"ALLRSHAT", 

+ are GREATER THAN RS" I'l 

+ - 1 ))) 

WRITEC1,4444) ALFACALF), 

! REAL(MAXREP)) * jqO 

WIKU Itw) W( AIF) ) 

ELSE 

WRITE(1,4444) ALFA(ALF), 

> REAL(MAXREP)) * 100 
SlTEa’44471 ^^^■^^ALF.TRUQNTCALF)-!) 

620 END IF ^ 

400 CONTINUE 

----- 

IF(K.EQ.5) then 
DO 401 ALF=1,MAXALF 
TRUQNT(ALF) = 0 

DO 501 REPS=1, MAXREP 

501 CO.NTINUr^^^^^^ "" ^ShBBG - RSHTBR(ALF,REPS) 

DO 601 REPS=1, MAXREP 

IF(ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT(ALF) = REPS 

+ P^i^4)')’'^" ^O^IDENCE LIMIT IS: ”, 

■*■ (TRUQNT(ALF) / REAL(MAXREP)) * 100 

GO TO 621 

ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQNT(ALF) = REPS 
GO TO 611 

ELSE 
END IF 

601 CONTINUE 

611 IE(TRUQNTCALF). EQ. 0. ) THEN 

WRITEC1,4443) ALFACALF) 
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raiTEC1,'( " " ,/ " THE SMALLEST'', 

+ '' DIFFERENCE BETVEEN RSBRDG AND RSHTBR IS: " , 

+ FIO.5)') DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,'(' ' " ,/’' ALL RSHTBR'', 

+ "ARE GREATER THAN RSBRDG'')’) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNI’(ALF) - 1))) 
+ THEN 

WRITEC1,4444) ALFACALF), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITECl,4447) 

621 END IF 

401 CONTINUE 
ELSE 
END IF 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

•kii-k ***5^**** RHTSTR ******** *** 

IFCFLAG.EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNTTCALF) = 0 

DO 5500 REPS=1, MAXREP 

DIFFCREPS) = RS - RHTSTRCALF,REPS) 

5500 CONTINUE 

DO 6600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C'' '',/''TRUE CONFIDENCE LIMIT IS:'', 

+ F6.4)') 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 6620 

ELSEIFCDIFFCREPS).LT.O. ) THEN 
TRUQNTCALF) = REPS 
GO TO 6610 

ELSE 
END IF 

6600 CONTINUE 

6610 IFCTRUQNTCALF).EQ.0.) THEN 

WRITEC1,4443) ALFACALF) 

WRITEC1,'C" ",/"THE SMALLEST" , 

+ '' DIFFERENCE BETWEEN RS AND RHTSTR IS: " , 

+ F9.5)') DIFFCMAXREP) 

ELSEIFCTRUQNTCALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFACALF) 

WRITEC1,'C" ",/"ALL RHTSTR" , 

+ "ARE GREATER THAN RS" )') 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE. ABSCDIFFCTRUQNTCALF) - 1))) 
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+ THEN- 

WRITE(1/4444) ALFA(ALF), 

+ (TRUQNT(ALF) / REALCMAXREP)) * 100. 

WRITEC1,4448) RHTSTR(ALF,TRUQNT(^ALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,^444) ALFA(ALF), 

+ ((TRUQNT(ALF)-l) / REAL(KAXREP)) * 100. 

WRITEC1,4448) RHTSTRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION*** 

IFCPRKT.EQ. 1) THEN 
I = 1 

185 WRITEC1,REPSHD) ALFACSELCTA), ALFACSELCTA), 

+ALFACSELCTB),ALFACSELCTB).ALFACSELCTA).ALFACSELCTA).ALFACSELCTB), 
+ALFACSELCTB) 

175 IFCI.GE. CMAXREP + 1) ) THEN 
GOTO 180 

ELSE 

IFC CI. EQ. 71). OR. Cl. EQ. 211). OR. CI. EQ. 351). OR. CI. EQ. 491). OR. 

+ CI.EQ.631).OR.CI.EQ.771).OR.CI.EQ.911).OR.CI.EQ.1051) ) THEN 

I = I + 70 
WRITEC1,'C" + ")') 

GOTO 185 

ELSE 

WRITEC1,3336) I. INTCDEGFRCI)), CHISQCl.I), QHTUPRCl,!), 

+ CHISQC2.I), QHTUPRC2,I) 

END IF 

IFC I + 70.LE.HAXREP) THEN 

WRITEC1,3337) 1+70,INTCDEGFRC1+70)).CHISQC1,1+70), 

+ QHTUPRC1,1+70).CHISQC 2.1+70).QHTUPRC 2.1+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 

9999 WRITEC 1,’("THE TOTAL NO OF REPS WAS:",18)') TOTREP 

WRITEC 1,'C "THE TOTAL NO OF EFFECTIVE REPS WAS: ",18)') LOOP 
WRITECl.'C'THE TOTAL NO OF NO FAILURE RUNS WAS: ' ,18)') ZFAILS 
WRITEC 1,’("AVERAGE NO. OF COMPONENTS PER REPLICATION WITH ", 
+"N0 FAILURES: " ,F5. 2) ’ ) ZFPREP / MAXREP 
WRITECl.'C'THE TOTAL NO OF RUNS WITH FAILURES WAS: ",I8)') FAILS 
0008 FORMAT (/ 3X,'C l’,5X,'C 2', 

+5X,'C 3',5X,'C 4',5X,'C 5',5X,'C 6',5X,'C 7',5X, 

+'C 8',5X,'C 9',5X,'C 10',4X,'C ll',4X, 

+'C 12',4X,'C 13',4X,'C 14',4X,'C 15') 

0009 FORMATC/lX.'REP NO',2X,'DF',IX,'CHISQRC',F4.3,')',IX, 
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+ 'QHTUPR(’ ,F4.3,')' IX.'CHISQRC, F4. 3,')', 1X,'QHTUPR(' ,F4. 5,’)' , 
+2X, 'REP NO' ,2X, 'DF^1X 'CHISQR( \F4. 3,' ) ', IX, 

+'QHTUPR(',F4.3,')',1X,'CHISQR(',F4. 3, )',IX/QHTUPR(',F4. 3,')'/) 
0001 FORMAT (///'UNIFORM RANDOM DEVIATES ARE: ) 

0002 FORMAT (///'BERNOULLI TRIALS ARE:') 

0003 FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT:') 

0004 FORMAT (///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT:') 

0005 FORMAT (///'TOTAL NTJMBER OF MISSION TESTS: ') 

0006 FOR.MAT (///'ESTIMATED WEIGHTS FOR EACH COMPONENT: ') 

0007 FORMAT (///'Q I FOR EACH COMPONENT:') 

1111 FORMAT (15F8,5) 


',F4.3,') CONFIDENCE ', 
',F4.3,') CONTIDENCE ', 


2222 FORMAT (/IX,15(14,4X)) 

3333 FORMAT (/IX,15(14,4X)) 

3334 FORMAT (/15F8.5) 

3335 FORMAT (/1X,15(I4,4X)) 

3336 FORMAT (T3,14,T9,13,T13,F11.5,T27,F8.5,T39,F11.5,T53,F8. 5) 

3337 FORMAT ( '+\T67,14,T73,13,T77,F11. 5 ,T91 ,F8. 5 ,T103,F11. 5,T117,F8. 5) 

4442 FORMAT (' ',///'THE RESULTING (1 - ‘,F4.3,') CONTIDENCE ', 

+'LIMIT IS:',T50,' 00.000 ') 

4443 FORMAT (' ',///‘tHE RESULTING (1 - ',F4.3,') CONFIDENCE ', 

+'LIMIT IS:',T50,'l00.0000') 

4444 FORMAT (' ',///'THE RESULTING (1 - ',F4.3,') CONTIDENCE ', 

+'LIMIT IS:',T50,F8.4) 

4445 FORMAT (' ',/'THE RSHAT VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4446 FORMAT ('; ',/'(FIRST NEGATIVE DIFFERENCE)') 

4447 FORMAT (' ',/'(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)') 

4448 FORMAT (' ',/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4449 FORMAT (' ',/'THE RSHTBR VALUE CLOSEST TO RSBPvDG IS: ',T51,F8.5) 

5555 FORMAT (' ',///'THE ',14,'(1-',F4.3,') QUANTILE IS:',T49,F8. 3) 

5556 FORMAT (' ',/'THE VALUE OF RSHAT FOR THAT QUANTILE IS: ’,T51,F8.5) 

5557 FORMAT (' ',/'THE DIFFERENCE(RS - RSHAT) IS:',T51,F8.5) 

5656 FORMAT (' ',/'THE VALUE OF RHTSTR FOR THAT QUANTILE IS:',T51,F8. 5) 

5666 FORMAT (' ',/'THE VALUE OF RSHTBR FOR THAT QUANTILE IS:',T51,F8.5) 

5657 FORMAT (' ',/'THE DIFFERENCE(RS - RHTSTR) IS:',T51,F8. 5) 

5667 FORMAT (' ',/'THE DIFFERENCE(RS - RSHTBR) IS:',T51,F8. 5) 

5755 FORMAT (' ',///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR', 

+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY ', 

+'''RHTSTR'' IS COMPUTED') 

6666 F0RMA7 C 

+' AwiV****!^*** RUN INPUT SETTINGS ****T:************jr*******»V**5T5r*' ^ 

6667 FORMAT (' './/'NUMBER OF REPLICATIONS: ',T50,14) 

6668 FORMAT (' './/'NUMBER OF COMPONENTS:',T50,14) 

6669 FORMAT (' './/'SYSTEM RELIABILITY FUNCTION:',T50,'SERIES') 

6699 FORMAT (' './/'SYSTEM RELIABILITY FUNCTION: ',T50,'BRIDGE') 

6670 FORMAT (' './/'MASTER UNRELIABILITY USED:',T50,F8.5) 

6671 FORMAT (' './/'INPUT WEIGHTS(A SUB I"S):') 

6674 FORMAT ('+',///'***********************************************', 

+' UN RESULT S***************************’'^'***********' , 

+ ' *** 51 ?********** ' ) 

6675 FORMAT ('+' , ///' **********#r5?r***********5ff*****fryf*****5*****5*******' ^ 
+ ' ESTIMATE ERRORS ****A********5V***********5Jr*******************' ^ 

6676 FORMAT ('+' ,///' ****5V*5V*****5rr****5V**5V****5t****5V**5V****5V*****5V** ' ^ 

+ ' TRUE CONFIDENCE LIMITS ***^-****''f***********5’f******************' , 
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+' **************' ■) 

end ^ 
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APPENDIX E. FORTRAN CODE FOR THE PREFERRED LOWER 
CONFIDENCE LIMIT FOR SYSTEM RELIABILITY (SERIES-PARALLEL 
SYSTEM WITH A 2/3 COMPONENT 


PROGRAM ZEYSCN 


******************** 

* 


TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE * 

ZERO FAILURES ALLOWED; NO SCALING * 

AUTHOR: E. F. BELLINI, LT, USN * 

MODIFIED BY: LT VALERIE A. COVINGTON,USN (MAR 90) 

DATE: NOV 89 * 

THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 
RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 
OF THEIR COMPONENTS * 

* 

IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 

TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT * 
SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC'.* 
THE REXX EXEC PROGRAM * 

'Bl' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT 
FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE * 

INDEX COUNTERS I, J, K OF THE ‘Bl’ EXEC ONE CAN RUN ANY USER- * 

SPECIFIC RUN FROM JUST ONE RU-N TO ALL 12. * 

it 

VARIABLES USED * 


* 

AHATI 

* 

AI 

it 

ALFA 

it 

BIGF 

it 

CHISQ 

it 

C1C15 

it 

DEGFR 

* 

DELBRG 

* 

DELSTR 

it 

DELTAR 

* 

DIFF 

* 

EPS 

* 

ERROR 

* 

FAILS 

* 

FI 

* 

FLAG 

* 

INC 

it 

KEYl 

* 

KEY2 

* 

KEY3 

* 

KEY4 

it 

KK 

it 

LOOP 

it 

MAXALF 


WEIGHT ESTIMATES FOR EACH COMPONENT * 
INPUT WEIGHTS FOR EACH COMPONENT * 
LEVELS OF SIGNIFICANCE * 
TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHI-SQUARE RANDOM VARIABLE VALUE * 
FORMAT LABEL * 
DEGREES OF FREEDOM * 
DIFFERENCE FOR BRIDGE SYSTEM * 
DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DIFFERENCE FOR SERIES SYSTEM * 
DIFFERENCE (TRUE REL. - ESTIMATED REL.) * 
SMALL QUANTITYCCONSTANT) * 
PARAMETER FOR IMSL ROUTINE * 
COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
1 IF ALL COMP. HAVE SA.ME NO. OF MISSION TESTS * 
INCREMENT STEP SIZE FOR ROUTINE USMNMX * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INTECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 


ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
COUNTS NO. OF REPLICATION PERFORMED * 
MAX NO. OF SIGNIFICANT LEVELS DESIRED(ARRAY SIZING)* 
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MAXREP : 


MAXkUN : 

•it 

MSTRQ ; 

it 

MULT : 

it 

N : 

it 

NIMAX : 

it 

NIMIM : 

* 

NINDX : 

* 

NIREAL : 

* 

NMAX : 

* 

NPRNT : 

* 

PRNT : 

* 

QHATI : 

it 

QHTMAX : 

it 

QHTUPR : 

i: 

QI : 

it 

QINDX : 

it 

QUA.NTL : 

it 

REPSHD : 


RKTSTR : 

it 

RS 

it 

RSBRDG : 

* 

RSHAT : 

it 

RSHTBR : 

it 

SEED : 

it 

SELCTA : 

it 

SELCTB : 

it 

SORT : 

it 

SUMNAI : 

it 

TEMP ; 

it 

TOTREP : 

* 

TRANBR : 

it 

TRANSQ : 

it 

TRANSR : 

it 

TRIALS : 

it 

TRNSTR : 

it 

TRUQNT : 

it 

UNIRV : 

it 

ZFAILS : 

it 

ZFPREP : 


MAX NO. OF REPLICATIONS ■ * 

MAX NO. OF PROGRAM ITER.4TIONS ALLOWED * 

MASTER U.\’RELIABILm*(USED WITH Al’S TO CALC. Ql’S) * 
MULTIPLIER FOR RANDOM NO. GENERATOR SRND * 

NO. OF MISSION TEST FOR EACH CCMPO.NENT * 

MAX NO. OF MISSION TESTS * 

MIN NO. OF MISSIO.N TESTS * 

INDEX NO. OF MAX NO. OF MISSION TESTS * 

NO. OF MISSION TESTS TRANSFORMED TO REAL * 

MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 
FLAG FOR DETAILED REPORT OUTPUT * 

SAME AS ABOXTCPARAMETER) * 

UNRELIABILITY ESTIM.ATES FOR EACH COMPONENT * 

LARGEST QHATI * 

UPPER LIMIT ON SYSTEM UNRELIABILITY * 

INPUT U.NRELIABILiy FOR EACH COMPONENT * 

INDEX * 

QUA.NTILE * 

REPLICATIONS HEADING FORMAT NUMBER 
SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) * 
TRUE SERIES SYSTEM RELIABILITY 

TRUE BRIDGE SYSTEM RELIABILITY * 

SERIES SYSTEM RELIABILITY ESTIMATE * 

BRIDGE SYSTEM RELIABILITY ESTIhLATE * 

PARAMETER * 

SIGNIFICANCE LEVEL SELECTION * 

SIGNIFICANCE LEIDL SELECTION * 

PA.RAMETER FOR ROUTINE SR.ND * 

SUTJ OF THE PRODUCT OF NI’S A.ND AI'S * 

TEMPORARY ARRAY * 

TOTAL NUMBER OF PROGRAM ITERATIONS * 

TEMPORARY ARRAY * 

TEMPORARY ARRAY * 

TEMPOPvARY ARRAY * 

BER.NOULLI TRIALS ARRAY (2-DIM) * 

TEMPORARY ARRAY * 

TRUE QUANTILE * 

UNIFORM RANDOM DEVIATES (2-DIM) * 

TOTAL NUMBER OF REPLICATIONS WITO ZERO FAILURES * 
iNO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION * 


rrATfrArrfr********; 


PARAMETER (KK=10,MAXALF=2,NPRNT=0) 

PARA.METER CMAXREP=1000, MAXRUN=2000, EPS=.000001) 

REAL*4 UNIRV( 15,1000) ,TE.MP( 1000) ,QI(KK),AI(KK) ,AHATI(KK) 

REAL*4 QHATI(KK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFA(KAXALF) 
REAL*4 DF( 5 ), AALFA( 5 ), SUMNAI, RSH.AT( MAXAL?, MAXREP), RS 
REAL*4 KEY1(M.AXREP) ,KEY2(MAXREP) ,KEY3( MAXREP I ,TRNSTR(MAXREP) 
REAL*4 DEGFR(MAXREP), QHrUPR(MAXALF,MAXREP),CHISQ(HAX.ALF,MAXREP) 
REAL*4 QUPA1(MAXP£P), QUPA2(M.AXREP),R.HTSTR(MAXALF,MAXREP) 

REAL*4 DELTAR(MAXALF), TRANSQ(MAXREP),TRANSR(MAXREP),DIFF(MAXREP) 
REAL*4 DELSTR(M.AXALF) ,NIMIN,NIMAX,NIREAL(KK) 

REAL*4 RSKTBR(MAXALF,MAX'REP) ,DELBRG(HA.XALF) ,KEY4(MAXREP) 

REAL*4 TRA.NBR(MAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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REAL*4 AVGN,SUC,STUD,PTEMP(10),FDEG1,FDEG2,P3,S3,S4 

INTEGER SEED, MULT, SORT, TRIALSC15,1000), BIGF, FI(KK), N(KK) 
INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTREP 
INTLwuR RF.PSwn, SELCTB, ALF, FLAG, LOOP,PRNT 

INTEGER QUANTLCMAXALF), TRUQNT(MAXALF),ZFAILS,FAILS,INC 
INTEGER NTEST,FCT,HFI 

CiiARACTER*8 LOOPSO(MAXREP) 

DATA SEED/123457/, MULT/1/, INC/1/ 

DATA AALFA/.01,.05,.9,.95,.99/, DF/1,5,10,30,40/ 

DATA ALFA/. 20,. 050/ 

DATA SORT/0/ 

ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 

* CALL COMPRS 
PRNT = NPRNT 

DO 12 1=1,KK 

AI(I) = 9999. 

N(I) = 99999999 
12 CONTINUE 

READ(03,*)K,MSTRQ 

DO 11 1=1,K 
REAP(03,*) AI(I),N(I) 

11 CONTINUE 


IF(K.NE.5) THEN 

WRITE(1,'(’'WARNING: BRIDGE STRUCTURE '', 
+"ONLY USES THE FIRST 5 COMPONENTS")') 

ELSE 
END IF 


***INITIALIZE TIE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, *** 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 

* AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 

^’W-ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NUMBER *** 


DO 172 ALF=1,MAXALF 

DO 173 REPS=1,MAXREP 

QHTUPR(ALF,REPS) = 0. 

RSHAT(ALF,REPS) = 0. 

RHTSTRCALF,REPS) = 0. 

RSHTBR(ALI',REFS) = 0. 

LOOFSO(REPS)=' ****-.v* ' 

173 CONTINUE 

172 CONTINUE 

***SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS**** 
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FLAG=1 

DO 50 1=1,K -1 

IF((N(I) - N(I+l)).NE.O) THEN 
FLAG=0 

ELSE 
END IF 
50 CONTINUE 

PRINT 'FLAG IS:'. FLAG 

***MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)*** 

ZFPREP = 0, 

ZFAILS = 0 
FAILS = 0 
TOTREP = 0 
LOOP = 0 

10 IF(LOOP.LT.MAXREP) THEN 
LOOP = LOOP + 1 
IFCTOTREP. LT. MAXRUN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT’'^** 
***0F THE SUBROUTINE SHSORT *** 

DO 95 REPS=1, MAXREP 
KEYl(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEYACREPS) = REPS 
95 CONTINUE 

***CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE*** 
***CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO.*** 

CALL IHAX(N,K,NMAX,NINDX) 


***CALCULATE THE QI'S FROM THE GIVEN MASTER Q AND THE Al'S*** 

DO 115 1=1, K 

QI(I) = MSTRQ * AI(I) 

115 CONTINUE 

S3=QI(3) 

S4=QI(4) 

DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINUE 

*** DRAW UNIFORM (0,1) RV’S AND CONVERT TO BERNOULLI TRIALS*** 
DO 130 1=1, K 
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CALL SRNDCSEED, TEMP, N(I), MULT, SORT) 
DO 135 J=l, N(I) 

UNIRV(I,J) = TEMP(J) 

IF (UN1RV(I,J).LE. 1 - OKI)) THEN 
TRiALS(I,J) = 0 

ELSE 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
130 CONTINUE 


***CALCULATE THE NO. OF FAILURES FOR EACH C0MP0NT;NT*** 


DO 150 1=1, K 
FKI) = 0 
150 CONTINUE 
lONECT = 0 

***CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES*** 
BIGF = 0 
DO 155 1=1, K 

DO 160 J=l, N(I) 

FI(I) = FI(I) + TRIALS(I,J) 

160 CONTINUE 

IF (FI(I) .EQ. 0) THEN 
ZFPREP=ZFPREP+1 
ELSE 
ENDIF 

BIGF = BIGF + FI(I) 

QHATI(I) = REAL(FI(I)) / N(I) 

155 CONTINUE 

***Change for Series-Parallel System with 2 out of 3 component*** 


IF (FI(4) .EQ. 0) THEN 
ZFPREP=ZFPREP - 1 
ENDIF 

BIGF = BIGF - FI(4) 

FI(4) = 0 
DO 161 1=1,N(4) 

HFI=0 

DO 162 J=l,2 

CALL SRNDCSEED, PTEMP(J), N(4), MULT, SORT) 
IF (PTEMP(J) .GT. 1-S4) THEN 
HFI = HFI + 1 
ENDIF 
162 CONTINUE 

IF (HFI .EQ. 2) THEN 
FI(4) = FI(4) + 1 
ENDIF 
161 CONTINUE 

BIGF = BIGF + FI(4) 

IF (FI(4) .EQ. 0) THEN 
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ENDIF 

IF CFI(3) .EQ. 0) THEN 

nwVr ■ 

DO 163 r=l,N(3) 

HFI=0 

DO 164 J=l,3 

164 CONTINUE 

IF (HFI .GE. 2) THEN 

endif'^' = * 1 

163 CONTINUE 

BIGF = BIGF + FI(3) 

__.2IPREP = ZFPREP + 1 
£*jLot 

ENDIF 

DO 19 1=1,K 

^QI(I) = MSTRQ * Aidi 

QHATI(I) = REAL(FI(I)j / N(i) 

19 CONTINUE 

QI(4) = QI(4)**2 

P3 = 1.-QI(3) 

DO 156 1=1,K 

156 CONTINU^^^^ IO^''ECT^IO^’ECT+l 

*’VoVCASE WHERE NO COMPONENTS HAVE ANY FAILURES>v*vv 

IF(BIGF.EQ, 0) THEN 

LOOPSO(I,OOP)=' *ZERO* ' 

ZFAILS = ZFAILS + 1 
AVGN=0. 0 
DO 200 1=1,K 

200 

AVGN=AVGN/REAL(K) 

DO 205 ALF=1, MAXALF 

ALFA(/-«*(i, /wGN) 

EISE ’®'®’^<*^-K>OP>=AiFA' ■J)**(l./N(l)) 
^ END IF 

PEINT MODE’ .LOOP, 'RS,(AT' .REHATCALF.LOOP) 



205 CONTINUE 

DEGFR(LOOP) = 2. 

GO TO 10 

ELSE 

FAILS = FAILS + 1 
END IF 

***COUNTS NUMBER OF COMPONENTS THAT FAIL RECORDS NO. COMPT TESTS*** 
FCT=0 

DO 202 1=1,K 

IF (FI(I) .NE. 0) THEN 
FCT=FCT+1 
NTEST=N(I) 

ENDIF 

202 CONTINUE 

***FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES*** 

CALL RMAXCQHATI, K, QHTMAX, QINDX) 

***CALCULATE THE AHAT SUB I’S (WEIGHT ESTIMATES)*** . 

***IF COMPONENT HAS NO FAILURES AHAT SUB I IS ZERO*** 


SUMNAI = 0. 

DO 165 1=1, K 

AHATI(I) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + N(I) * AHATI(I) 

165 CONTINUE 

***1 COMPONENT FAILURE SERIES SYSTEM*** 

IF (FCT .EQ. 1) THEN 

LOOPSO(LOOP)=' *ONECF* ' 

DO 305 ALF=1,MAXALF 
SUC=REAL(NTEST-BIGF) 

FDEG1=2.*(REAL(BIGF)+1. ) 

FDEG2=2. *SUC 

STUD=FIN(1. -ALFA(ALF),FDEG1,FDEG2) 

RSHATC ALF,LOOP)=SUC/(SUC+(REAL(BIGF)+1. )*STUD) 

* PRINT *,'SUC=',SUC,'FAIL=’,BIGF,'NTEST=',NTEST 

* PRINT *,'FIN=' STUD 

* PRINT *,'ALFA=',ALFA(ALF),'RSHAT=’,RSHAT(ALF,LOOP) 

IF (FLAG .EQ.1) THEN 

RHTSTR(ALF,LOOP)=RSHAT(ALF,LOOP) 

ELSE 

ENDIF 

305 CONTINUE 
ENDIF 

***CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY*** 
DEGFR(LOOP) = 2 * (1 + BIGF) 

DO 170 ALF=1, MAXALF 

CALL MDCHKl - ALFA(ALF) ,DEGFR(LOOP) ,CHISQ(ALF,LOOP), ERROR) 
QHTUPR(ALF,LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAI) 
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IF(FLAG.EQ.1) THEN 

RHTSTR(ALF,LOOP) = 1 -(CHISQ(ALF,LOOP) / REAL(2*N(1))) 

ELSE 
END IF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES*** 
IF (FCT .NE. 1) THEN 

CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF,LOOP)) 
ENDIF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE*** 

IF (lONECT .NE. 1) THEN 

CALL RHTBRGC QHTUPR(ALF,LOOP),AHATI,K,RSHTBRC ALF,LOOP)) 
ENDIF 

170 CONTINUE 

***EXACTLY 1 COMPONENT FAILS AND REDUNDANT COMPONE.NT*** 

IF ((lONECT .EQ. 1) .AND. (K .EQ. 5)) THEN 
DO 207 1=1, K 

NIREALCI) = REAL(N(I)) 

207 CONTINUE 

CALL USMN.MX( NIREAL, K, INC,NIMIN, NIMAX) 

DO 206 ALF=1,MAXALF 

RSHTBR(ALF,LOOP)=ALFA(ALF)**(1. /NIMIN) 

206 CONTINUE 
ENDIF 

***THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN*** 

ELSE 

WRITE(1,'( " " ,/ " PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+'' ALLOWED OF: " ,16)') TOTREP 
GOTO 9999 

END IF 
GOTO 10 
END IF 

C WRITE(2,'( "UNSORTED RSHAT 1 IS: " ,/10(F8. 5))' ) 

C +(RSHAT(1,L00P), L00P=1, MAXREP) 

C WRITE(2,'( "UNSORTED RSHAT 2 IS: " ,/10(F8. 5)) ' ) 

C +(RSHAT(2,LOOP), L00P=1, MAXREP) 

C IFCFLAG.EQ.1) THEN 

C WRITE(2,'( " UNSORTED RHTSTR 1 

C +(RHTSTR(1,L00P), L00P=1, MAXREP) 

C WRITE(2,'(''UNSORTED RHTSTR 2 

C +(RHTSTR(2,L00P), L00P=1, MAXREP) 

C ELSE 

C END IF 

C IF(K.EQ.5) THEN 

C WRITE(2,’( " UNSORTED RSHTBR 1 

C +(RSHTBR(l,LOOP)j L00P=1, MAXREP) 

C \sRITE( 2, ’ ( ' ’ UNSORTED RSHTBR 2 


IS: '•,/10(F8.5))') 
IS: ’',/10(F8.5))’) 

IS: ’',/10(F8.S))') 
IS: " ,/10(F8.5))') 
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C +(RSHTBR(2,L00P), L00P=1, MAXREP) 

C ELSE 

C END IF 

C WRITE (2,'( " ZERO AND ONE FAILURE REPS:’',/10(AS))') 

C + (LOOPSO(LOOP),LOOP=l,NAXREP) 

***SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL)*** 


DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPR(ALF,REPS) 
TRANSR(REPS) = RSHAT(ALF,REPS) 
TRNSTR(REPS) = RHTSTR(ALF,REPS) 
TRANBR(REPS) = RSHTBR(ALF,REPS) 
800 CONTINUE 

CALL SHS0RT(TRANSQ,KEY1,MAXREP) 

CALL SHSORT(TRANSR,KEY2,MAXREP) 

CALL SHSORT(TRNSTR,KEY3,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 
RSHAT(ALF,REPS) = TRANSR(REPS) 
RHTSTR(ALF,REPS) = TRNSTR(REPS) 
RSHTBR(ALF,RF.PS) = TRANBR(REPS) 
900 CONTINUE 

700 CONTINUE 

***PRINT OUTPUT REPORT HEADINGS*** 


WRITEC1,6666) 

WRITEC1,6667) MAXREP 
WRITEC1,6668) K 
WRITEC1,6669) 

IFCK.EQ. 5) THEN 
WRITEC1,6699) 

ELSE 
END IF 

V.'RITEC 1,6670) MSTRQ 
WRITEC1,6671) 

WRITEC1,C1C15) 

WRITEC1,3334) AI 
WRITEC1,0007) 

WRITEC1,C1C15) 

WRITEC1,3334) QI 
WRITEC1,0005) 

W'RITEC1,C1C15) 

WRITEC1,3335) N 
WRITEC1,6674) 

***COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN.CSERIES SYSTEM)*** 
*** AND FOR THE 5-COMPONENT BRIDGE STRUCTURE*** 


CALL RSRSCQI,K.RS) 

WRITECl.’C" ,////"THE TRUE SERIES SYSTEM ", 
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+’'RELIABILITY VALUE IS: " ,T51,F8. 5)') RS 
CALL RBRIDG(QI,K,RSBRDG) 

IF(:K.EQ. 5) THEN 

WRITE(1,'('' '’ , 1111 ''THE TRUE BRIDGE STRUCTURE '', 
+"RELIABILITY VALUE IS; " ,T51,F8. 5) ' ) RSBRDG 
ELSE 
END IF 

WRITEC1,6675) 

***COMPUTE THE DIFFERENCE 'DELTAR' BTWN. RS AND RSHAT OF THE THEO*** 
***RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY)*** 

IFCFLAG.EQ. 1) THEN 
TOITEC1,5755) 

ELSE 
END IF 

DO 450 ALF=1, MAXALF 

QUANTL(ALF) = MAXREP * (1 - ALFA(ALF)) 

DELTAR(ALF) = RS - RSHAT(ALF,QUANTL(ALF)) 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF)) 

IF(FLAG. EQ. 1) THEN 

DELSTR(ALF) = RS > RHTSTR(ALF,QUANTL(ALF)) 

WRITEC1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
WRITEC1,5656) RHTSTRCALF,QUANTLCALF)) 

WRITEC1,5657) DELSTRCALF) 

ELSE 
END IF 

IFCK.EQ. 5) THEN 

DELBRGCALF) = RSBRDG - RSHTBRCALF,QUANTLCALF)) 
WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5666) RSHTBRCALF,QUANTLCALF)) 

■ WRITEC1,5667) DELBRGCALF) 

ELSE 
END IF 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 

WRITEC1,5556) RSHATC ALF,QUANTLC ALF)) 

WRITEC1,5557) DELTARCALF) 

450 CONTINUE 

PRINT*, 'QUANTLCl) IS:', QUANTLCl) 

PRINT*, 'QUANTLC2) IS:', QUANTLC2) 

***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

*** ******** RSHAT ******* *** 

WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - RSHATCALF,REPS) 

500 CONTINUE 

DO 600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C" ",/"TRUE CONFIDENCE LIMIT IS:", 
+ F8.4)') 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 
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GO TO 620 

ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQNT(ALF) = REPS 
GO TO 610 

ELSE 
END IF 
CONTINUE 

IF(TRUQNT(ALF) .EQ. 0. ) THEN 
WRITE(1,4443) ALFA(ALF) 


V?RITE(1,'(’ 


'THE SMALLEST' 


h ” DIFFERENCE BETWEEN RS AND RSHAT IS:'',F10.5)') DIFF( 

h MAXREP) 

ELSEIF(TRUQNT(ALF).EQ.1.) THEN 
TOITEC1,4442) ALFA(ALF) 

WRITEC1,’( " " ,/ " ALL RSHAT’', 

I- "ARE GREATER THAN RS" ) ’ ) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.ABSCDIFF(TRUQNT(ALF) - 1))) 
h THEN 

WRITEC1,4444) ALFA(ALF), 

»■ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITE(1,4445) RSHAT(ALF,TRUQNT(ALF)) 

WRITEC1,4446) 

ELSE 

WDITEC1,4444) ALFACALF), 

»■ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTC ALF)-1) 

WRITEC1,4447) 

END IF 
CONTINUE 


***FIND the true CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** **A5V**** rsHTBR C BRIDGE) ******* *** 


IFCK.EQ. 5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 501 REPS=1, MAXREP 

DIFFCREPS) = RSBRDG - RSHTBRCALF,REPS) 

CONTINUE 

DO 601 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,’C" "./"TRUE CONFIDENCE LIMIT IS:", 
f F8.4)') 

f CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 621 

ELSEIFCDIFFCREPS). LT. 0.) THEN 
TRUQNTCALF) = REPS 
GO TO 611 

ELSE 
END IF 
CONTINUE 

IFCTRUQNTCALF).EQ.0.) THEN 
WRITEC1,4443) ALFACALF) 


WRITEC1,'C' 


'THE SMALLEST' 
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+ ’' DIFFERENCE BETCEEN RSBRDG AND RSHTBR IS:'’, 

+ F10.5)’) DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,’(” ALL RSHTBR” , 

+ '' ARE GREATER THAN RSBRDG'’)') 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITE(1,4444) ALFA(ALF), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4449) RSHTBR(ALF,TRUQNT(ALF)) 
y^T?ITE( 1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

621 END IF 

401 CONTINUE 
ELSE 
END IF 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*7r5V ***iV**5vyf RHTSTR *******5V *5V* 


IFCFLAG.EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 5500 REPS=1, MAXREP 

DIFFCREPS) = RS - RHTSTRCALF,REPS) 

5500 CONTINUE 

DO 6600 REPS=1, MAXREP 

IF(ABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C " " ,/ ” TRUE CONFIDENCE LIMIT IS: ” , 

+ F8.4)') 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 6620 

ELSEIFCDIFFCREPS). LT. 0. ) THEN 
TRUQNTCALF) = REPS 
GO TO 6610 

ELSE 
END IF 

6600 CONTINUE 

6610 IFCTRUQNTCALF).EQ.O.) THEN 

WRITEC1,4443) ALFACALF) 

WRITECl,'C" THE SMALLEST”, 

+ *' DIFFERENCE BETWEEN RS AND RHTSTR IS: ” , 

+ F9.5)') DIFFCMAXREP) 

ELSEIFCTRUQNTCALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFACALF) 

WRITEC1,'C" ”,/”ALL RHTSTR”, 

+ "ARE GREATER THAN RS” )' ) 

ELSEIFCABSCDIFFCTRUQNT(ALF))).LE.ABSCDIFFCTRUQNTCALF) - 1))) 
+ THEN 
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WRITE(1,A444) ALFA(AEF), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4448) RHTSTRCALF.TRUQNTCALF)) 

WRITE(1,4446) 

ELSE 

WRITE(1,4444) ALFA(ALF), 

+ ((TRUQNT(ALF)-l) / REAL(MAXREP)) * 100. 

WRITEC1,4448) RHTSTRCALF,TRUQNT(ALF)-1) 

WRITE(1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION*** 

IFCPRNT.EQ. 1) THEN 
1 = 1 

185 WRITEC1,REPSHD) ALFA(SELCTA), ALFA(SELCTA), 

+ALFA(SELCTB),ALFA(SELCTB),ALFA(SELCTA),ALFA(SELCTA),ALFA(SELCTB), 
+ALFA(SELCTB) 

175 IF(I.GE. (MAXREP + 1)) THEN 
GOTO 180 

ELSE 

IF( (I.EQ. 71).0R. (I.EQ.211).0R. (I. EQ. 351). OR. Cl. EQ. 491). OR. 

+ (I.EQ.631).OR.(I.EQ.771).OR.(I.EQ.911).OR.(I.EQ. 1051) ) THEN 

I = I + 70 
U’RITECl,'(" + ")') 

GOTO 185 

ELSE 

WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,I), QHTUPR(1,I), 

+ CHISQ(2,I), QHTUPR(2,I) 

END IF 

IF(I + 70. LE.MAXREP) THEN 

WRITEC1,3337) 1+70,INT(DEGFR(1+70)),CHISQC1,1+70), 

+ QHTUPRC1,1+70),CHISQC 2,1+70),QHTUPR(2,1+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 

9999 WRITEC1,'("THE TOTAL NO OF REPS WAS: ",I8)’) TOTREP 

WRITEC 1,'C "THE TOTAL NO OF EFFECTIVE REPS WAS: ".18)') LOOP 
WRITEC1,'("THE TOTAL NO OF NO FAILURE RUNS WAS: '^I8)') ZFAILS 
WRITE(1,'( "AVERAGE NO. OF COMPONENTS PER REPLICATION WITH ", 

+''NO FAILURES:'’,F5.2)') ZFPREP / MAXREP 

WRITEC 1,'( "THE TOTAL NO OF RQ-NS WITH FAILURES WAS:", 18)') FAILS 
0008 FOR.MAT (/ 3X,'C l',5X,'C 2', 

+5X,'C 3',5X,'C 4',5X,'C 5',5X,'C 6',5X,'C 7',5X, 

+'C 8',5X,'C 9',5X,'C 10',4X,'C ll',4X, 

+'C 12',4X,'C 13',4X,'C 14’',4X,'C 15') 

0009 F0RMAT(/1X,'REP NO',2X,'DF',1X,'CHISQRC’,F4.3,')',1X, 

+'QHTUPRC',F4.3,')',IX,'CHISQRC’,F4.3,')',IX,'QHTUPRC',F4. 3,' ) ' , 
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+2X,'REP NO’,2X,'DF',IX 'CHISQRC',F4. 3.’)’,IX, 

+’QHTUPR(',F4.3,')’,1X, CHISQR(’.F4.3, )',IX,'QHTUPR(’,F4. 3,')’/) 
0001 FORMAT (///'UNIFORM RANDOM DEVIATES ARE;') 

0002 FORMAT (///'BERNOULLI TRIALS ARE:') 

-0003 FORMAT (///’TOTAL NO. OF FAILURES FOR EACH COMPONENT:’) 

0004 FORMAT (///’ESTIMATED UNRELIABILITY FOR EACH COMPONENT:’) 

0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS; ') 

0006 FORMAT (///’ESTIMATED WEIGHTS FOR EACH COMPONENT: ’) 

0007 FORMAT (///’Q I FOR EACH COMPONENT: ’) 

1111 FORMAT (15F8.5) 


2222 FORMAT (/IX,15(14,4X)) 

3333 FORMAT (/1X,15(I4,4X)) 

3334 FORMAT (/15F8.5) 

3335 FORMAT (/IX,15(I4,4X)) 

3336 FORMAT (T3,I4,T9,I3,T13,F.ll. 5,T27,F8. 5,T39,F11. 5,T53,F8. 5) 

3337 FORMAT (' + ',T67,I4,T73,I3,T77,FIl. 5,T91,F8. 5,T103,F11. 5,T117,F8. 5) 

4442 FORMAT (’ ',///‘THE RESULTING (1 - ',F4.3,’) CONFIDENCE ', 

+'LIMIT IS:’,T50,’ 00.000 ’) 

4443 FORMAT-(' ',///‘THE RESULTING (1 - ’,F4.3,’) CONFIDENCE ’, 


+'LIMIT IS: 
4444 FORMAT (' 


IS:',T50,’100.0000') 

(' ’,///‘TKE RESULTING- (1 - 


+'LIMIT IS: ',T50,F8.4) 


',F4. 3 ,’) CONFIDENCE ', 
’,F4. 3,’) CONFIDENCE ’, 


4445 FORMAT (' ',/'rriE RSHAT VALUE CLOSEST,TO RS IS: ',T51,F8.5) 

4446 FORMAT (' ',/’(FIRST NEGATIVE DIFFERENCE)'' ) 

4447 FORMAT (' ',/'(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)’) 

4448 FORMAT (' ',/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4449 FORMAT C ',/'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T5l,F8.-5) 

5555 FORMAT (’ ’,///'THE ’,14,’(!-’,F4. 3,’) QUANTILE IS:’,T49,F8. 3) 

5556 FORMAT (’ ’,/'THE VALUE OF RSHAT FOR THAT QUANTILE IS: ’,T51,F8.5) 

5557 FORMAT (' ',/'THE DIFFERENCECRS - RSHAT) IS; ',T51,F8.5) 

5656 FORMAT (’ ’,/’THE VALUE OF RHTSTR FOR THAT QUANTILE IS:’,T51,F8. 5) 

5666 FORMAT (' ’,/'THE VALUE OF RSHTBR FOR THAT QUANTILE IS: ’,T51,F8.5) 

5657 FORMAT (' ',/'THE DIFFERENCE(RS - RHTSTR) IS: ’,T51,F8.5) 

5667 FORMAT (' ’,/'TKE DIFFERENCE(RS - RSHTBR) IS:’,T51,F8.5) 

5755 FORMAT (’ ’,///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR’, 

+’ ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY’ ’, 

+' ''RHTSTR'' IS COMPUTED') 

6666 FORMAT (’+’,'’t******?^***************************************' ^ 
+'**;■;******** INPUT SETTINGS ******************************** , 

+ ' *lV!V!V********5V*;V*5V5V!V* ' J 

6667 FORMAT (’ ',//'NUMBER OF REPLICATIONS: ’,T50,14) 

6668 FORMAT (' ’,//'NUMBER OF COMPONENTS: ',T50,14) 

6669 FORMAT (' './/'SYSTEM RELIABILITY FUNCTION:',T50,'SERIES') 

6699 FORMAT (’ ’.//’SYSTEM RELIABILITY FUNCTION:’,T50,’BRIDGE’) 

6670 FORMAT (’ //'MASTER UNRELIABILITY USED: ’,T50,F8.5) 

6671 FORMAT (’ './/'INPIT WEIGHTS(A SUB l"S); ') 

6674 FORMAT ('+’,///'*’Wf***T'f**********************^f*****************’, 

+’**R result S****************************************’ 

+' ' ) 

6675 format ('+',///'************************************************', 
+’ ESTIMATE ERRORS **********************^f******y:***************' ^ 

+ ***?V***5V*5V*5V***' ) 

6676 FORMAT ('+',///'*****************************************vr*****'^ 
+' TRUE CONFIDENCE LIMITS ***************************************', 
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APPENDIX F. FORTRAN CODE FOR THE PREFERRED LOWER 
CONFIDENCE LIMIT FOR SYSTEM RELIABILITY (PARALLEL 

SYSTEM) 


PROGRAM ZFYSCN . 

* 

* TITLE; BINOMIAL INTERVAL ESTIMATION PROCEDURE 

* ZERO FAILURES ALLOWED; NO SCALING 

* AUTHOR: E. F. BELLINI, LT, USN 

* MODIFIED BY: LT VALERIE A. COVINGTON,USN (MAR 90) 

* DATE: NOV 89 


* THIS PROGRAM COMPUTES THE TRUE CONTJDENCE LE\TL FOR THE ESTIMATE * 

* RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 

* OF THEIR COMPONENTS * 


* IN ITS PRESENT CONTIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 

* TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT * 

* SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC'.* 

* THE REXX EXEC PROGRAM * 

* 'Bl' CALLS THE I^ ' FILES TO BE READ AND NAMES THE 12 OUTPUT * 

* FILES RESULTING TiE 12 CONSECUTIVE RUNS. BY EDITING THE * 

* INDEX COU.NTERS 1, X OF THE 'Bl' EXEC ONE CAN RUN ANY USER- * 

* SPECIFIC RUN FRO.M JUST ONE RUN TO ALL 12. * 

'i: 

* VARIABLES USED * 


i: 

is 



"it 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


AHATI : WEIGHT ESTIMATES FOR EACH COMPONENT * 
AI : INPUT WEIGHTS FOR EACH COMPONENT * 
ALFA ; LEVELS OF SIGNIFICANCE * 
BIGF : TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHISQ : CHI-SQUARE RANDOM VARIABLE VALUE * 
C1C15 ; FORMAT L.ABEL * 
DEGFR ; DEGREES OF FREEDOM 

DELBRG : DIFFERENCE FOR BRIDGE SYSTEM * 
DELSTR : DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DELTAR ; DIFFERENCE FOR SERIES SYSTEM * 
DIFF ; DIFFERENCE (TRUE REL. - ESTIMATED REL.) * 
EPS : SMALL QUANTITY(CONSTANT) * 
ERROR : PARAMETER FOR IHSL ROUTINE * 
FAILS : COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
FI : NO. OF FAILURES FOR EACH COMPONTNT(ALL MISSION TST)* 


FLAG : 1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INC : INCREMENT STEP SIZE FOR ROUTINE USMNm * 
KEYl : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEY2 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEYS : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEY4 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KK : ARRAY SIZING PARAMETER FOR THE MAX NO OF COHPONDNTS* 
LOOP : COUNTS NO. OF REPLICATION PERFORMED * 


MAXALF : MAX NO. OF SIGNTFICA.NT LEVELS DESIRED(ARR/iY SIZING)* 
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* 

MAXREP 


MAX NO. OF REPLICATIONS 

it 

* 

MAXRUN 


MAX NO. OF PROGRAM ITERATIONS ALLOWED 

* 

* 

MSTRQ 


MASTER UNRELIABILITYCUSED WITH AI’S TO CALC. Ql'S) 

* 

* 

MULT 


MULTIPLIER FOR RANDOM NO. GENERATOR SRND 

* 

* 

N 


NO. OF MISSION TEST FOR EACH COMPONENT * 

* 

NIMAX 


MAX NO. OF MISSION TESTS 

* 

* 

NIMIM 


MIN NO. OF MISSION TESTS 

* 

* 

NINDX 


INDEX NO. OF MAX NO. OF MISSION TESTS 

* 

* 

NIREAL 


NO. OF MISSION TESTS TRANSFORMED TO REAL 

* 

* 

NMAX 


MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 

* 

* 

NPRNT 


FLAG FOR DETAILED REPORT OUTPUT 

* 

* 

PRNT 


SAME AS ABOVECPARAMETER) 

* 

* 

QKATI 


UNRELIABILITY ESTIMATES FOR EACH COMPONENT 

* 

* 

QHTMAX 


LARGEST QHATI 

* 

* 

QHTUPR 


UPPER LIMIT ON SYSTEM UNRELIABILITY 

* 

* 

QI 


INPUT UNRELIABILIY FOR EACH COMPONENT 

* 

* 

QINDX 


INDEX 

* 

* 

QUANTL 


QUANTILE 

* 

* 

REPSHD 


REPLICATIONS HEADING FORMAT NUMBER 

* 

* 

RHTSTR 


SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 

* 

* 

RS 


TRUE SERIES SYSTEM RELIABILITY 

* 

* 

RSBRDG 


TRUE BRIDGE SYSTEM RELIABILITY 

* 

* 

RSHAT 


SERIES SYSTEM RELIABILITY ESTIMATE 

* 

* 

RSHTBR 


BRIDGE SYSTEM RELIABILITY ESTIMATE 

* 

* 

SEED 


PARAMETER 

* 

iV 

SELCTA 


SIGNIFICANCE LEVEL SELECTION 

* 

* 

SELCTB 


SIGNIFICANCE LEVEL SELECTION 

* 

* 

SORT 


PARAMETER FOR ROUTINE SRND 

* 

Vf 

SUMNAI 


SUM OF m PRODUCT OF Nl'S AND Al'S 

* 

* 

TEMP 


TEMPORARY ARRAY 

* 

* 

TOTREP 


TOTAL NUMBER OF PROGRAM ITERATIONS 

* 

* 

TRANBR 


TEMPORARY ARRAY 

* 

* 

TRANSQ 


TEMPORARY ARRAY 

* 

* 

TRANSR 


TEMPORARY ARRAY 

* 

* 

TRIALS 


BERNOULLI TRIALS ARRAY (2-DIM) 

* 

* 

TRNSTR 


TEMPORARY ARRAY 

* 


TRUQNT 


TRUE QUANTILE 

* 

* 

UNIRV 


UNIFORM RANDOM DEVIATES (2-DIM) 

* 

* 

ZFAILS 


TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES 

* 

* 

* 

ZFPREP 


NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION 

* 

* 


************************************************************************ 


PARAMETER (KK=10,MAXALF=2,NPRNT=0) 

PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=.000001) 

REAL*4 UNIRV(15,1000),TEMP(1000),QI(KK),AI(KK),AHATI(KK) 

REAL*4 QHATI(KK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFA(MAXALF) 
P>EAL*4 DF(5) AALFA(5) ,SUMNAI,RSHAT(MAXALF,MAXREP) ,RS 
REAL*4 KEYl(MAXREP),KEY2(MAXREP),KEY3(MAXREP),TRNSTR(MAXREP) 
REAL*4 DEGFR(MAXREP), QHTUPR(MAXALF,MAXREP),CHISQ(MAXALF,MAXREP) 
REAL*4 QUPAl(MAXREP), QUPA2(MAXREP),RHTSTR(MAXALF,MAXREP) 

REAL*4 DELTARCMAXAIJ-), TRANSQ(MAXREP),TRANSR(MAXREP),DIFF(MAXREP) 
REAL*4 DELSTR(MAXALF),NIM1N,NIMAX,NIREAL(KK). 

REAL*4 RSHTBR(MAXALF,MAXREP)5DELBRG(MAXALF),KEY4(MAXREP) 

REAL*4 TRANBR(MAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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R£AL*4 AVON, SUC,STUD,PTEMP(10),FDEG1,FDEG2,P3,S3,S4 

INTEGER SEED, MULT, SORT, TRIALS(15,1000), BIGF, FI(KK), N(KK) 
INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTREP 
INTEGER G1C15, REPSHD, SELCTB, ALF, FLAG, LOOP,PRNT 
INTEGER QUANTL(MAXALF), TRUQNT(MAXALF),ZFAILS,FAILS,INC 
INTEGER NTEST,FCT,HFI 

CHARACTER*8 LOOPSO(MAXREP) 

DATA SEED/123457/, MULT/1/, INC/1/ 

DATA AALFA/.01,.05,.9,,95,.99/, DF/1,5,10,30,40/ 

DATA ALFA/. 20,. 050/ 

DATA SORT/0/ 

ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 

CALL COMPRS 
PRNT = NPRNT 


DO 12 1=1,KK 

AI(I) = 9999. 
N(I) = 99999999 
12 CONTINUE 

READ(03,*)K,HSTRQ 

DO 11 1=1,K 
READ(03,*) AI(I),N(I) 
11 CONTINUE 


IF(K.NE.5) THEN 

WRITE(1,’( " WARNING: BRIDGE STRUCTURE ’', 
+''ONLY USES THE FIRST 5 COMPONENTS'')') 

ELSE 
END IF 


***INITIALIZE THE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, *** 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 

* AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 

***ALL the COMPONENT MISSION TESTS ARE EQUAL IN NUMBER *** 


DO 172 ALF=1,MAXALF 

DO 173 REPS=1,MAXREP 

QHTUPRCALF,REPS) = 0. 

RSHATCALF,REPS) = 0. 

RHTSTRCALF,REPS) = 0. 

RSHTBR(ALF,REPS) = 0. 

LOOPSO(REPS)=' ****** ' 

173 CONTINUE 

172 CONTINUE 

***SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS**** 
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FLAG=1 

DO 50 1=1,K -1 

IFCCNCD - N(I+1)).NE, 0) THEN 
FMG=0 

ELSE 
END IF 
50 CONTINUE 

PRINT *, 'FLAG IS:', FLAG 

***MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)*** 

ZFPREP = 0. 

ZFAILS = 0 
FAILS = 0 
TOTREP = 0 
LOOP = 0 

10 IF(LOOP. LT. MAXREP) THEN 
LOOP = LOOP + 1 
IFCTOTREP. LT.MAXRUN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT*** 
***OF THE subroutine; SHSORT *** 

DO 95 REPS=1, MAXREP 
KEYl(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
95 CONTINUE 


***CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE*** 
***CALCUMTE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO.*** 

CALL IMAX(N,K,NMAX,NINDX) 


***CALCULATE THE Ql'S FROM THE GIVEN MASTER Q AND THE Al'S*** 

DO 115 1=1, K 

QI(I) = MSTRQ * AI(I) 

115 CONTINUE 
C S3=QI(3) 

S4=QI(1)**.20 
DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINTJE 

*** DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS*** 
DO 130 1=1, K 
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CALL SRNDCSEED, TEMP, N(I), MULT, SORT) 
DO 135 J=l, N(I) 

UNIRV(I,J) = TEMP(J) 

IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 
TRIALS(I,J) = 0 

ELSE 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
130 CONTINUE 


***CALCULATE THE NO. OF FAILURES FOR EACH COMPONENT*** 

DO 150 1=1, K 
FI(I) = 0 
150 CONTIN’UE 
lONECT = 0 

**^'CALCULATE THE F SUB I'S AND THE GRAND TOTAL NO. OF FAILURES 
BIGF = 0 
DO 155 1=1, K 

DO 160 J=l, N(I) 

FI(I) = FI(I) + TRIALS(I,J) 

160 CONTINUE 

IF (FI(I). .EQ. 0) THEN 
ZFPREP=ZFPREP+1 
ELSE 
ENDIF 

BIGF = BIGF + FI(I) 

QHATI(I) = REALCFKI)) / N(I) 

155 CONTINUE 

***Changes for parallel system*** 


IF (FI(1) .EQ. 0) THEN 
ZFPREP=ZFPREP - 1 
ENDIF 

BIGF = BIGF - FI(1) 

FI(1) = 0 
DO 161 1=1,N(l) 

HFI=0 

DO 162 J=l,5 

CALL SRNDCSEED, PTEMP(J), N(l), MULT, SORT) 
IF (PTEMP(J) .GT. 1-S4) THEN 
HFI = HFI + 1 
ENDIF 
162 CONTINUE 

IF (HFI .EQ, 5) THEN 
FI(1) = FI(l) + 1 
ENTDIF 

161 CONTINUE 

BIGF = BIGF + FI(1) 

IF (FI(1) .EQ. 0) THEN 
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ZFPREP = ZFPREP + 1 

ELSE 

ENDIF 

C IF (FI(3) .EQ. 0) THEN 
C ZFPREP=ZFPREP - 1 

C ENDIF 

C BIGF = BIGF - FI(3) 

C FI(3) = 0 

C DO 163 1=1,N(3) 

C HFI=0 

C DO 164 J=l,3 

C CALL SRND(SEED, PTEMP(J), N(3), MULT, SORT) 

C IF (PTEMP(J) .GT. 1-S3) THEN 

C HFI = HFI + 1 

C ENDIF 

C164 CONTINUE 

C IF (HFI .GE. 2) THEN 

C FI(3) = FI(3) + 1 

C ENDIF 

C163 CONTINUE 

C BIGF = BIGF + FI(3) 

C IF (FI(3) .EQ. 0) THEN 
C ZFPREP = ZFPREP + 1 

C ELSE 
C ENDIF 

DO 19 1=1,K 
01(1) = MSTRQ * AI(I) 

QHATI(I) = REAL(FI(I)) / N(I) 

19 CONTINUE 
C QI(4) = QI(4)**2 
C P3 = 1.-QI(3) 

C QI(3) = 1.-((3.*P3**2*QI(3))+(P3**3)) 

***COUNTS NUMBER OF COMPONENTS THAT HAVE FAILED*** 

DO 156 1=1,K 

IF (FI(I) .NE. 0) I0NECT=I0N’ECT+1 
156 CONTINUE 

***CASE Vf'HERE NO COMPONENTS HAVE ANY FAILURES*** 

IF(BIGF.EQ. 0) THEN 

LOOPSO(LOOP)=' *ZERO* ’ 

ZFAILS = ZFAILS + 1 
AVGN=0.0 
DO 200 1=1,K 

AVGN=AVGN+REAL(N(I)) 

200 CONTINUE 

AVGN=AVGN/REAL(K) 

DO 205 ALF=1, MAXALF 

RSHAT(ALF,LOOP)= ALFA(ALF)**(1. /AVGN) 
IF(FLAG. EQ. 1) THEN 

RHTSTR(ALF,L00P)=ALFA(ALF)**(1. /N(l)) 

ELSE 
ENT) IF 

* PRIN’^ *,LOOP' ,LOOP, 'RSHAT' ,RSHAT(ALF,L00P) 
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205 


CONTINUE 

DEGFR(LOOP) = 2. 

GO TO 10 

ELSE 

FAILS = FAILS + 1 
END IF 

***COUNTS NUMBER OF COMPONENTS THAT FAIL RECORDS NO, COMPT TESTS*** 
FCT=0 

DO 202 1=1,K 

IF (FI(I) .NE. 0) THEN 
FCT=FCT+1 
NTEST=N(I) 

ENDIF 

202 CONTINUE 

***FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES*** 

CALL RMAXCQHATI, K, QHTMAX, QINDX) 

IF (LOOP .EQ. 1) THEN 
ENDIF 

***CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES)*** 

***IF COMPONENT HAS NO FAILURES AHAT SUB I IS ZERO*** 

SUMNAI = 0. 

DO 165 1=1, K 

AHATI(I) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + N(I) * AKATI(I) 

165 CONTINUE 

***1 COMPONENT FAILU-RE SERIES SYSTEM*** 

IF (FCT .EQ. 1) THEN 

LOOPSO(LOOP)=’ *ONECF* ’ 

DO 305 ALF=1,MAXALF 
SUC=REAL(NTEST-BIGF) 

FDEG1=2. *(REAL(BIGF)+1. ) 

FDEG2=2.*SUC 

STUD=FIN(1.-ALFA(ALF),FDEG1,FDEG2) 

RSHAT(ALF,LOOP)=SUC/(SUC+(REAL(BIGF)+1. )*STUD) 

* PRINT *,'SUC=’,SUC,'FAIL=',BIGF,'NTEST=',NTEST 

* PRINT *,'FIN=’ STUD 

* PRINT *,'ALFA=',ALFA(ALF),'RSHAT=',RSHAT(ALF,LOOP) 

IF (FLAG ,EQ.1) THEN 

RHTSTR(ALF,LOOP)=RSHAT(ALF,LOOP) 

ELSE 

ENDIF 

305 CONTINUE 
ENDIF 

***CALCULATE 1 REPLICATION OF UPPR ALFA C.L. ON SYSTEM RELIABILITY*** 
DEGFR(LOOP) = 2 * (1 + BIGF) 

DO 170 ALF=1, MAXALF 
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CALL MDCHI(1 - ALFA(ALF),DEGFR(LOOP),CHISQ(ALF,LOOP), ERROR) 
QHTUPR(ALF,LOOP') = CHISQ(ALF,LOOP) / (2 * SUMNAI) 

IFCFLAG.EQ. 1) THEN 

RHTSTR(ALF,LOOP) = 1 -(CH1SQ(ALF,LOOP) / REAL(2*N(1))) 

ELSE 
END IF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN SERIES*** 
IF (FCT .NE. 1) THEN 

CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHAT(ALF,LOOP)) 
ENDIF 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE*** 

IF (lONECT .NE. 1) THEN 

CALL RHTBRGC QHTUPR(ALF,LOOP),AHATI,K,RSHTBR(ALF,LOOP)) 
ENDIF 

170 CONTINUE 

***EXACTLY 1 COMPONENT FAILS ANT) REDUNDANT COMPONTiNT*** 

IF ((lONECT .EQ. 1) .AND. (K .EQ.5)) THEN 
DO 207 1=1, K 

NIREAL(I) = REAL(N(I)) 

207 CONTINUE 

CALL USMNMXC NIREAL,K,INC,NIMIN,NIMAX) 

DO 206 ALF=1,MAXALF 

RSHTBRCALF,LOOP)=ALFA(ALF)**(1. /NIMIN) 

206 CONTINUE 
ENDIF 

***THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN*** 

ELSE 

WRITE(1,'( " " ,/ " PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 

+ '' ALLOWED OF: " ,16)') TOTREP 
GOTO 9999 

END IF 
GOTO 10 
END IF 

C WRITE(2,'( "UNSORTED RSHAT 1 IS: " ,/10(F8. 5))' ) 

C +(RSKAT(1,L00P), LOOP=l, MAXREP) 

C WRITE(2,'( "UNSORTED RSHAT 2 IS: " ,/10(F8. 5))' ) 

C +(RSHAT(2,L00P), LOOP=l, MAXREP) 

C IFCFLAG. EQ. 1) THEN 

C WRITE(2,'( "UNSORTED RHTSTR 1 IS: " ,/10(F8. 5))' ) 

C +CRHTSTR(l,LOOP), LOOP=l, MAXREP) 

C VRITE(2,'(''UNSORTED RHTSTR 2 IS: " ,/10(F8.5))') 

C +(RHTSTR(2,LOOP), LOOP=l, MAXREP) 

C ELSE 

C ENT) IF 

C IFCK.EQ.S) THEN 

C WRITE(2,'( "UNSORTED RSHT3R 1 IS: " ,/10(F8. 5))') 
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+(RSHTBR(1,L00P), L00P=1, MAXREP) 

WRITE(2,'(''UNSORTED RSHTBR 2 IS: " ,/10(F8.5))') 
+(RSHTBR(2,LOOP), LOOP=l, MAXREP) 

ELSE 
END IF 

WRITE (2,'("ZERO AND ONE FAILURE REPS: " ,/10(A8))' ) 

+ (LOOPSO(LOOP),LOOP=1,MAXREP) 

***SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL)*** 

DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPR(ALF,REPS) 

TRANSR(REPS) = RSHAT(ALF,REPS) 

TRNSTR(REPS) = RHTSTR(ALF,REPS) 

TRANBR(REPS) = RSHTBR(ALF,REPS) 

800 CONTINUE 

CALL SHSORT(TRANSQ,KEyi,MAXREP) 

CALL SHS0RT(TRANSR,KEY2,MAXREP) 

CALL SKSORT(TRNSTR,KEYS,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 

RSHAT(ALF,REPS) = TRANSR(REPS) 

RHTSTR(ALF,REPS) = TRNSTR(REPS) 

RSHTBR(ALF,REPS) = TRANBR(R£PS) 

900 CONTINUE 

700 CONTINUE 

***PRINT OUTPUT REPORT HEADINGS*** 

WRITE(1,6666) 

WRITE(1,6667) MAXREP 
WT1ITE( 1,6668) K 
WRITEC1,6669) 

IF(K. EQ. 5) THEN 
WRITE(1,6699) 

ELSE 
END IF 

WRITE!1,6670) MSTRQ 
WRITE!1,6671) 

WRITE(1,C1C15) 

WRITE!1,3334) AI 
WRITEC1,0007) 

WRITE!1,C1C15) 

WRITE!1,3334) QI 
WRITEC1,0005) 

WRITE!1,C1C15) 

WRITE!1,3335) N 
WRITE!1,6674) 

***COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN.!SERIES SYSTEM)*** 
*** AND FOR THE 5-COMPONENT BRIDGE STRUCTURE*** 
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CALL RSRS(QI,K.RS) 

WRITEC 1,' ( " ' ‘ .////"THE TRUE SERIES SYSTEM ' ' , 

+ "RELIABILITY VALUE IS; " ,T51,F8. 5) ’ ) RS 
CALL RBRIDG(QI,K,RSBRDG) 

IF(K.EQ.5) THEN 

WRITEC 1,' ( " " ,////"THE TRUE BRIDGE STRUCTURE ’ ’ , 

+ "RELIABILITY VALUE IS: " ,T51,F8. 5)' ) RSBRDG 
ELSE 
END IF 

WRITEC1.6675) 

***COMPUTE THE DIFFERENCE 'DELTAR* BTVN. RS AND RSHAT OF THE THEO*** 
***RETICAL QUANTILE GIVEN BY ALFACMUST USE SORTED RSHAT ARRAY)*** 

IFCFLAG.EQ. 1) THEN 
WRITEC1,5755) 

ELSE 
END IF 

DO A50 ALF=1, MAXALF 

QUANTLCALF) = MAXREP * Cl “ ALFACALF)) 

DELTARCALF) = RS - RSHATCALF,QUANTLCALF)) 

DELBRGCALF) = RSBRDG - RSHTBRCALF.QUANTLCALF)) 

IFCFLAG.EQ. 1) THEN 

DELSTRCALF) = RS - RHTSTRCALF,QUANTLCALF)) 
WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5656) RHTSTRCALF,QUANTLCALF)) 

WRITEC1,5657) DELSTRCALF) 

ELSE 
END IF 

IFCK.EQ. 5) THEN 

DELBRGCALF) = RSBRDG - RSHTBRCALF.QUANTLCALF)) 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 
WRITEC1,5666) RSHTBRC ALF,QUANTLC ALF)) 

WRITEC1,5667) DELBRGCALF) 

ELSE 
END IF 

WRITEC1,5555) MAXREP, ALFACALF), REALCQUANTLCALF)) 

WRITE C1,5 5 5 6) RSHATCALF,QUANTLC ALF)) 

WRITEC1,5557) DELTARCALF) 

450 CONTINUE 

PRINT*, 'QUANTLC 1) IS:', QUANTLC 1) 

PRINT*, 'QUANTLC2) IS:', QUANTLC2) 

***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

■Mrit ******** RSHAT *^?***** *** 

WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - RSHATCALF,REPS) 

500 CONTINTJE 

DO 600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC 1,' C " " TRUE CONTIDENCE LIMIT IS: " , 
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F8.4)') 

(TRUQNT(ALF) / REAL(MAXREP)) * 100. 

GO TO 620 

ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQNT(ALF) = REPS 
GO TO 610 

ELSE 
END IF 
CONTINUE 

IF(TRUQNT(ALF).EQ. 0. ) THEN 
WRITEC1,4443) ALFA(ALF) 

WRITE(1,'C' '',/"THE SMALLEST”, 
h '' DIFFERENCE BETVEEN RS AND RSHAT IS:‘',F10.5)') DIFFC 

(• MAXREP) 

ELSEIF(TRUQNT(ALF).EQ.1.) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITEC 1,'(” ”,/”ALL RSHAT”, 
t- " ARE GREATER THAN RS” ) ’ ) 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.ABSCDIFFCTRUQNTCALF) - 1))) 
h THEN 

WRITEC1,4444) ALFACALF), 

I- CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

h CCTHUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

END IF 
CONTINUE 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** vv******,v rsHTBR CBRIDGE) * * ** * ** *** 


IFCK.EQ. 5) THEN 
DO 401 ALF=l,MAXArxF 
TRUQNTCALF) = 0 

DO 501 REPS=1, MAXREP 

DIFFCREPS) = RSBRDG - RSHTBRCALF,REPS) 

CONTINUE 

DO 601 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C” ”,/"TRUE CONFIDENCE LIMIT IS: 
f F8.4)' ) 

f CTRUQNTCALF) / REALCMAXP^P)) * 100. 

GO TO 621 

ELSEIFCDIFFCREPS). LT.0.) THEN 
TRUQNTCALF) = REPS 
GO TO 611 

ELSE 
END IF 
CONTINUE 

I?(TRUQNTCALF).EQ.O.) THEN 
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WRITEC 1,4443) ALFA(ALF); 

WRITE(1,'C" ”,/"THE SMALLEST", 

+ '’ DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS: " , 

+ F10.5)') DIFF(MAXREP) 

ELSEIFCTRUQNT(ALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITEC 1,’(" ",/"ALL RSHTBR", 

+ " ARE GREATER THAN RSBRDG*')’) 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.ABSCDIFFCTRUQNTCALF) “ 1))) 
+ THEN 

WRITEC1,4444) ALFACALF), 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

621 END IF 

401 CONTINUE 
ELSE 
END IF 


***FIND THE TRUE CONFIDENCE LE'.TEL OF THE SYSTEM REL. ESTIMATE*** 
*** ******** RHTSTR ******** *** 


IFCFLAG. EQ. 1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 5500 REPS=1, MAXREP 

DIFFCREPS) = RS - RHTSTRCALF,REPS) 

5500 CONTINTJE 

DO 6600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUONTCALF) = REPS 

WRITEC1,'C’' " ,/’‘TRUE CONFIDENCE LIMIT IS: " , 
+ F8.4)') 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 6620 

ELSEIFCDIFFCREPS).LT.0.) THEN 
TRUQNl'CALF) = REPS 
GO TO 6610 

ELSE 
• END IF 

6600 CONTINUE 

6610 IFCTRUQNTCALF).EQ.O. ) THEN 

WRITEC1,4443) ALFACALF) 

WRITEC1,*C" ",/"THE SMALLEST", 

+ ’' DIFFERENCE BETWEEN RS AND RHTSTR IS; " , 

+ F9.5)') DIFFCMAXREP) 

ELSEIFCTRUQNTCALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFACALF) 

WRITEC1,’C" “,/"ALL RHTSTR", 

+ "ARE GREATER THAN RS" ) ’ ) 
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ELSEIF(ABS(DIFFCTRUQNT(ALF))).LE.ABSCDIFF(TRUQNT(ALF) - 1))) 
+ IHEN 

WRITE(1,4444) ALFA(ALF), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4448) RHTSTR(ALF,TRUQNT(ALF)) 

WRITE(1,4446) 

ELSE 

WRITE(1,4444) ALFA(ALF),- 

+ ((TRUQNT(ALF)-l) / REAL(MAXREP)) * 100. 

WRITE(1,4448) RHTSTR(ALF,TRUQNT(ALF)-1) 

WRITE(1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION*** 

IF(PRNT.EQ. 1) THEN 
1 = 1 

185 WRITE(1,REPSHD) ALFA(SELCTA), ALFA(SELCTA), 

+ALFA(SELCTB),ALFA(SELCTB),ALFA(SELCTA),ALFA(SELCTA),ALFA(SELCTB), 
+ALFA(SELCTB) 

175 IFd.GE. (MAXREP + 1)) THEN 
GOTO 180 

ELSE 

IF( (I.EQ. 71). OR. (I.EQ. 211). OR. (I. EQ. 351). OR. (I.EQ. 491). OR. 

+ (I.EQ.631).OR. (I.EQ.771).OR. (I.EQ.911).OR.(I.EQ. 1051) ) THEN 

I = I + 70 
WRITEd,'(”+”)’) 

GOTO 185 

ELSE 

WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,I), QHTUPR(1,I), 

+ CHISQ(2,I), QHTUPR(2,I) 

END IF 

IF(I + 70. LE. MAXP^P) THEN 

WRITE(1,3337) 1+70,INT(DEGFR(1+70)),CHISQ(1,1+70), 

+ QHTUPR(1,1+70),CHISQ(2,1+70),QHTUPR(2,1+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 

9999 WRITE(1,'("THE TOTAL NO OF REPS WAS; ",I8)') TOTREP 

WRITE(1,'("THE TOTAL NO OF EFFECTIVE REPS WAS: ".18)') LOOP 
WRITE(1,'("THE TOTAL NO OF NO FAILURE RUNS WAS:'‘,18)') ZFAILS 
WRITE(1,'("AVERAGE NO. OF COMPONENTS PER REPLICATION WITH ", 

+''NO FAILURES: '',F5. 2)*) ZFPREP / MAXREP 

WRITE(1,'("THE TOTAL NO OF RUNS WITH FAILURES WAS; ",I8)') FAILS 
0008 FORMAT (/ 3X,’C l',5X,'C 2’, 

+5X,'C 3’,5X,’C 4',5X,'C 5’.5X,'C 6',5X,'C 7’,5X, 

+'C 8'jSX,'C 9',5X,'C 10’,4X,'C 11',4X, 

+'C 12‘,4X,'C 13',4X,'C 14',4X,'C 15') 
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0009 F0RMAT(/1X, ’REP NO',2X 'DF',1X 'CHISQRC',F4. 3,')',1X, 

+'QHTUPRC',F4. 3,')',IX,'CHISQRC ,F4.3,')'.IX,'QHTUPRC',F4. 3,')', 
+2X, 'REP NO' ,2X, 'DF^IX 'CHISQR( ' ,F4. 3) ^1X, 

+'QHTUPRC' ,F4.3,')' ,1X,'CHISQR(' ,F4.3,')’,1X,'QHTUPR(' ,F4.3,')'/) 
0001 FORMAT (///’UNIFORM RANDOM DEVIATES ARE:') 

0002 FORMAT (///’BERNOULLI TRIALS ARE:') 

0003 FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT: ') 

0004 FORMAT (///’ESTIMATED UNRELIABILITY FOR EACH COMPONENT: ') 

0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS:') 

0006 FORMAT (///’ESTIMATED VEIGHTS FOR EACH COMPONENT: ') 

0007 FORMAT (///’Q I FOR EACH COMPONENT: ') 

1111 FORMAT (15F8.5) 


2222 FORMAT (/IX,15(14,4X)) 

3333 FORMAT (/iX,15(14,4X)) 

3334 FORMAT (/15F8.5) 

3335 FORMAT (/1X,15(I4,4X)) 

3336 FORMAT (T3.I4,T9,I3,T13,F11.5,T27,F8. 5,T39,F11. 5,T53,F8. 5) 

3337 FORMAT ('+^T67,I4,T73,I3,T77,F11. 5.T91,F8. 5,T103,F11. 5,T117,F8. 5) 

4442 FORMAT (' ’,///'THE RESULTING (1 - ^F4.3,') CONFIDENCE ’, 

+’LIMIT IS:',T50,' 00.000 ’) 

4443 FORMAT (' ',///*THE RESULTING (1 - ’,F4.3,') CONFIDENCE ', 

+'LIMIT IS: ’,T50,'100.0000’) 

4444 FORMAT (' ',///‘THE RESULTING (1 - ',F4.3,') CONFIDENCE ', 

+’LIMIT IS: ',T50,F8.4) 

4445 FORMAT (’ ',/'THE RSHAT VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4446 FORMAT (' ’./’(FIRST NEGATIVE DIFFERENCE)') 

4447 FORMAT (’ ’./’(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)’) 

4448 FORMAT (’ ’./’THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 

4449 FORMAT (' './’THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T51.F8.5) 

5555 FORMAT (’ ’.///’THE ’,14,’(1-',F4.3,') QUANTILE IS:',T49,F8.3) 

5556 FORMAT (' './'THE VALUE OF RSHAT FOR THAT QUANTILE IS: ',T51,F8.5) 

5557 FORMAT (' './’THE DIFFERENCE(R3 - RSHAT) IS: ',T51,F8. 5) 

5656 FORMAT (' ’./'THE VALUE OF RHTSTR FOR THAT QUANTILE IS:’,T51,F8. 5) 

5666 FORMAT (' './’THE VALUE OF RSHTBR FOR THAT QUANTILE IS:'.T51,F8.5) 

5657 FORMAT (' ’./’THE DIFFERENCE(RS - RHTSTR) IS:',T51,F8.5) 

5667 FORMAT (' ’./'THE DIFFERENCE(RS - RSHTBR) IS:',T51,F8.5) 

5755 FORMAT (' './//'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR’. 

+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY '. 

+"'RHTSTR'I I S COMPl^p'). 

+ 'ir!rkMr/Hrlr>z*:* RUN INPUT SETTINGS , 


4 -«» ) 

6667 FORMAT (' ’.//'NUMBER OF REPLICATIONS: ',T50,14) 

’.//’NUMBER OF COMPONENTS: '.T50,I4) 

'.//'SYSTEM RELIABILITY FUNCTION:'.T50,'SERIES') 
'.//’SYSTEM RELIABILITY FUNCTION: ’,T50.'BRIDGE') 
'.//’MASTER UNRELIABILITY USED: ',T50.F8.5) 


6668 FORMAT (' 

6669 FORMAT (' 
6699 FORMAT (' 

6670 FORMAT (' 

6671 FORMAT (' 


'.//'INPUT WEIGHTS(A SUB l"S): ') 


6674 FORMAT ('+*,///' ^ 

6675 FORMAT ('+’ /// * ^ 

+1 ESTIMA TE ERRO RS 

6676 FORMAT ('+’,///'-^ 
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+1 TRUE CONF^ENCE LIMITS 
END ^ 






APPENDIX G. FORTRAN CODE FOR ALTERNATE PROCEDURE A 
FOR ESTIMATING THE LOWER CONFIDENCE LIMIT FOR SYSTEM 

RELIABILITY 

PROGRAM ZFYSCN 


* * 

* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE * 

* ZERO FAILURES ALLOWED; NO SCALING * 

* AUTHOR: E. F. BELLINI, LT, USN * 

* DATE: NOV 89 * 

* * 

* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 

* RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 

* OF THEIR COMPONENTS * 

* * 


* IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 

* TIMES- EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT * 

* SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC'.* 


* THE REXX EXEC PROGRAM * 

* 'Bl' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT * 

* FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE * 

* INDEX COUNTERS I, J, K OF THE 'Bl' EXEC ONE CAN RUN ANY USER- * 

* SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. * 

it it 


* VARIABLES USED 

it 


it 

* 


* 

* 

* 

* 

it 

* 

* 

* 

* 

* 

•* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

it 

it 

* 

* 


AHATI : WEIGHT ESTIMATES FOR EACH COMPONENT * 
AI : INPUT WEIGHTS FOR EACH COMPONENT * 
ALFA : LEVELS OF SIGNIFICANCE * 
BIGF ; TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHISQ ; CHI-SQUARE RANDOM VARIABLE VALUE * 
C1C15 : FORMAT LABEL * 
DEGFR : DEGREES OF FREEDOM * 
DELBRG : DIFFERENCE FOR BRIDGE SYSTEM * 
DELSTR : DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DELTAR : DIFFERENCE FOR SERIES SYSTEM * 
DIFF : DIFFERENCE (TRUE REL. - ESTIMATED REL. ) * 
EPS : SMALL QUANTITYCCONSTANT) * 
ERROR : PARAMETER FOR IMSL ROUTINE * 
FAILS : COUNTS NO. OF REPLICATIONS WITH AT LST. ? FAILURE * 
FI : NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
FLAG : 1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INC ; INCREMENT STEP SIZE FOR ROUIINE USMNMX * 
KEYl : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEY2 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEYS : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KEY4 : ARRAY OF INDECES FOR ROUTINE SHSORT * 
KK : ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
LOOP ; COUNTS NO. OF REPLICATION PERFORMED * 


MAXALF : MAX NO. OF SIGNIFICANT LEVELS DESIRED(ARRAY SIZING)* 
MAXREP : MAX NO. OF REPLICATIONS * 
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* 

MAXRUN 


MAX NO. OF PROGRAM ITERATIONS ALLOWED 

k 

* 

MSTRQ 


-MASTER: UNRELIABILITY(USED WITH AI'S TO CALC. QI'S) 

k 

* 

MULT 

• 

MULTIPLIER FOR RANDOM NO. GENERATOR SRND 

k 

* 

N 

* 

* 

NO. OF MISSION TEST FOR EACH COMPONENT * 

* 

NIMAX 


MAX NO. OF MISSION TESTS 

k 

* 

NIMIM 


MIN NO; OF MISSION TESTS 

* 

* 

NINDX 


INDEX NO. OF MAX NO. OF MISSION TESTS 

k 

* 

NIREAL 


NO. OF MISSION TESTS TRANSFORMED TO REAL 

k 


NMAX 


MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 

k 

* 

NPRNT 

♦ 

FLAG FOR DETAILED REPORT OUTPUT 

k 

* 

PRNT 

• 

• 

SAME AS ABOVECPARAMETER) 

k 

* 

QHATI 


UNRELIABILITY ESTIMATES FOR EACH COMPONENT 

k 

•k 

QHTMAX 

• 

LARGEST QHATI 

k 

* 

QHTUPR 

• 

• 

UPPER LIMIT ON SYSTEM UNRELIABILITY 

k 

■k 

QI 


INPUT UNRELIABILIY FOR EACH COMPONENT 

k 

★ 

QINDX 


INDEX 

k 

* 

QUANTL 


QUANTILE 

k 

* 

REPSHD 


REPLICATIONS HEADING FORMAT NUMBER 

k 

* 

RHTSTR 


SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 

k 

■k 

RS 


TRUE SERIES SYSTEM RELIABILITY 

k 

* 

RSBRDG 

• 

TRUE BRIDGE SYSTEM RELIABIIITY 

k 

■k 

RSHAT 

• 

* 

SERIES SYSTEM RELIABILITY ESTIMATE 

k 

* 

RSHTBR 

* 

• 

BRIDGE SYSTEM RELIABILITY ESTIMATE 

k 

* 

SEED 

• 

• 

PARAMETER 

k 

* 

SELCTA 

* 

• 

SIGNIFICANCE LEVEL SELECTION 

k 

* 

SELCTB 

• 

• 

SIGNIFICANCE LEVEL SELECTION 

k 

* 

SORT 

• 

« 

PAPJ^TER FOR ROUTINE SRND 

k 

* 

SUMNAI 

• 

* 

SUM OF THE PRODUCT OF Nl'S AND Al'S 

k 

k 

TEMP 

« 

TEMPORARY ARRAY 

k 

k 

TOTREP 

• 

* 

TOTAL NUMBER OF PROGRAM ITERATIONS 

k 

k 

TRANBR 

• 

• 

TEMPORARY ARRAY 

k 

k 

TRANSQ 

• 

• 

TEMPORARY ARRAY 

k 

k 

TRANSR 

• 

• 

TEMPORARY ARRAY 

k 

k 

TRIALS 

• 

BERNOULLI TRIALS ARRAY (2-DIM) 

k 

k 

TPvNSTR 

• 

• 

TEMPORARY ARRAY 

k 

k 

TRUQNT 

• 

• 

TRUE QUANTILE 

k 

k 

UNIRV 

* 

* 

UNIFORM RANDOM DEVIATES (2-DIM) 

k 

k 

ZFAILS 

• 

TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES 

k 

k 

k 

ZFPREP 


NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION 

k 

k 


PARAMETER (KK=10,MAXALF=2,NPRNT=0) 

PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=. 000001) 

REAL*4 UNIRVC15,1000),TEMP(1000),QI(KK),AI(KK),AHATI(KK) 

REAL*4 QHATICKK), NMAX, NNMAX, QHTMAX, CHISQR(5,5), ALFA(MAXALF) 
REAL*4 DF(5),AALFA(5),SUMNAI,RSHAT(MAXALF,MAXREP),RS 
REALM KEYKMAXREP) ,KEy2(MAXREP) ,KEY3(MAXREP) ,TRNSTR(MAXREP) 
REALM DEGFRCMAXREP), QHTUPR(MAXALF,MAXREP),CHISQ(MAXALF,MAXREP) 
REALM QUPAl(MAXREP), QUPA2(MAXREP),RHTSTR(MAXALF,MAXREP) 

REALM DELTAR(MAXALF), TRANSQ(MA)®EP),TRANSR(MAXREP),DIFF(MAXREP) 
REALM DELSTR(MAXALF),NIMIN,NIMAX,NIREAL(KK) 

REALM RSHTBR(MAXALF,MAXREP),DELBRG(MAXALF),KEY4(MAXREP) 

REALM TRANBR(MAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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INTEGER SEED, MULT, SORT, TRIAES(15,1000), BIGF, FICKK), N(KK)= 
INTEGER NIOTX, QINDX, ERROR, REPS, SELCTA,= SELCTK, TOTREP 
INTEGER C1C15, REPSHD, SELCTB, ALF, FLAG,. LOOP,PRNT 
INTEGER QU^LC MAXALF) , TRUQNTCMAXALF), ZFAILS,FAILS,INC 

DATA SEED/123A57/, MULT/1/, INC/1/ 

DATA AALFA/. 01,. 05,. 9,. 95,. 99/, DF/1,5,10',30,40/ 

DATA ALFA/. 20,. 050/ 

DATA SORT/0/ 

ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 

* CALL COMPRS 
PRNT = NPRNT 

DO 12 I=1,KK 

AI(I) = 9999. 

N(I) = 99999999 
12 CONTINUE 

READ(03,*)K,MSTRQ 

DO 11 1=1,K 
READ(03,*) AI(I),N(I) 

11 CONTINUE 


IF(K.NE.5) THEN 

WRITE01,'( "WARNING: BRIDGE STRUCTURE ", 

+"ONLY USES THE FIRST 5 COMPONENTS")') 

ELSE 
END IF 

//*** 


//VWf* 

DO 172 ALF=1,MAXALF 

DO 173 REPS=1,MAXREP 

QHTUPRCALF,REPS) = 0. 

RSHATCALF,REPS) = 0. 

RHTSTRCALF,REPS) = 0. 

RSHTBRCALF,REPS) = 0. 

173 CONTINUE 

172 CONTINUE 


***// INITIALIZE THE QHTUPR ARRAY OF UNRELIABILITY REPLICATIONS, 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS 

* AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN 
***// all the component MISSION TESTS ARE EQUAL IN NUMBER 


***// SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS**** 
FLAG=1 

DO 50 1=1,K -1 . 

IF((N(I) - N(I+l)).NE.O) THEN 
FLAG=0 

ELSE 
END IF 
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50 CONTINUE 

PRINT *, 'FLAG IS:', FLAG 

*^// MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)//*** 

ZFPREP = 0. 

ZFAILS== 0 
FAILS = 0 
TOTREP = 0 
LOOP = 0 

10 IF(LOOP.LT.MAXREP) THEN 
LOOP = LOOP + 1 
IFCTOTREP. LT. MA}®UN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***// FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT 

***// OF THE- SUBROUTINE SHSORT 
DO 95 REPS=1, MAXREP 
KEY1(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
95 CONTINUE 

***// calculate: nmAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE 

***// CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NO. //*** 

CALL IMAX(N,K,NMAX,NINDX) 


***// CALCULATE THE Ql'S FROM THE GIVEN MASTER Q AND THE AI'S 

DO 115 1=1, K 

QI(I) = MSTRQ * AI(I) 

115 CONTINUE 

DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINUE 

***// DRAW UNIFORM (0,1) RV'S AND CONVERT TO BERNOULLI TRIALS //*** 
DO 130 1=1, K 

CALL SRNDCSEED, TEMP, NCI), MULT, SORT) 

DO 135 J=l, N(I) 

UNIRV(I,J) = TEMPCJ) 

IF (UNIRV(I,J).LE. 1 - QI(I)) THEN 
TRIALS(I,J) = 0 

ELSE 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
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130 CONTINUE 


***/;/ CALCULATE THE NO. OF FAILURES FOR: EACH COMPONENT //*** 

DO 150 1=1, K 
FI(I) = 0 
150 CONTINUE 


***// CALCULAIE THE F SUB I'S AND THE GiiAND TOTAL NO. OF FAILURES 
BIGF =0 
DO 155 1=1, K 

DO 160 J=l, N(I) 

FI(I) = FI{I) + TRIALS(I,J) 

160 CONTINUE 

IF(FI(I).EQ.O) THEN 

ZFPREP = ZFPREP + 1 

ELSE 
END; IF 

CALCULATE THE QHAT SUB I’.S: F SUB I’S DIVIDED BY N SUB I'S 
QHATK I) = REALCFKD ) 7 N(I) 

BIGF = BIGF + FI(I> 

155 CONTINUE 

***// CASE WHERE NO^ COMPONENTS: HAVE ANY FAILURES //*** 

IF(BIGE.EO.O) THEN 

ZFAILS = ZFAILS + 1 
DO 200 1=1, K 

NIREAL(I) = REALCNd)) 

200 COi.TINUE 

CALL USMNMX(NIREAL,K,INC,NIMIN,NIMAX) 

DO 205 ALF=1, MAXALF 

CALL MDCHICl - ALFA(ALF) ,2. ,CHISQ(ALF,LOOP) ,ERROR-) 
RSHATCALF,LOOP)= 1 - (CHISQ(ALF,LOOP) / REAL(2 * NIMIN)) 
IFCFLAG.EQ. 1) THEN 

RHTSTRCALF,L00P)=1-(CKISQ(ALF,L00P) / REAL(2 * N(l))) 

ELSE 
END IF 

205 CONTINUE 

IFCPRNT.EQ. 1) THEN 
WRITEC1,0007) 

WRITE(1,C1C15) 

WRITEC1,3334) QI 
WRITEC1,0001) 

WRITEC1,C1C15) 

DO 141 J=1,NMAX 

WRITECljllll) CUNIRVCI,J), 1=1, K) 

141 CONTINUE 

WRITEC1,0002) 

WRITEC1,C1C15) 

DO 146 J=1,NMAX 

WRITEC1,2222) CTRIALSCI,J), 1=1, K) 

146 CONTINUE 

WRITEC1,0003) 

WRITEC1,C1C15) 
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WRITEC1,3333) FI 
WRITE(1,0005) 

WRITEC 1,C1C 15): 

WRITEC1,3335) N 
WRITEC1,0004) 

WRITEC1,C1C15) 

WRITEC1,3334) QHATI 

WRITEC1,’C/ THE MAXIMUM Q HAT SUB I IS:*', T40, F8.5)’) QHTMAX 
WRITEC1,’C/’’THE MAXI Q HAT SUB I IS ELMNT NO.: '’,T40,I5)’) QINDX 
WRITEC1,’C/’’THE GRAND TOTAL NO. OF FAILURES IS:’',T40, 15)') BIGF 
ELSE 
ENDIF 

DEGFRCLOOP) = 2. 

GO TO 10 

ELSE 

FAILS = FAILS + 1 
END IF 

***// FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES 
CALL RMAXCQHATI, K, QHTMAX, QINDX) 

***// PRINT THE RESULT OF THE MAIN OPERATING ELEMENTS OF THE PROGRAM 

IFCPRNT.EQ.1) THEN 
WRITEC1,0007) 

WRITEC1,C1C15) 

WRITEC1,3334) QI 
WRIIEC1,0001) 

WRITEC1,C1C15) 

DO 140 J=1,NMAX 

WRITEC1,1111) CUNIRVCI.J), 1=1, K) 

140 CONl'INUE 

WRITEC1,0002) 

WRITEC1,C1C15) 

DO 145 J=1,NMAX 

WRITEC1,2222) CTRIALSCI,J), 1=1, K) 

145 CONTINUE 

WRITEC1,0003) 

WRITEC1,C1C15) 

WRITEC1,3333) FI 
WRITEC1,0005) 

WRITEC1,C1C15) 

WRITEC1,3335) N 
WRITEC1,0004) 

WRITEC1,C1C15) 

WRITEC1,3334) QHATI 

WRITEC 1,'C/' THE MAXIMUM Q HAT SUB I IS:", T40, F8.5)') QHTMAX 
WRITEC1,'C/"THE MAXI Q HAT SUB I IS ELMNT NO.: ^' ,T40,I5)') QINDX 
WRITEC 1,’C/”THE GRAND TOTAL NO. OF FAILURES IS:",T40, 15)') BIGF 
ELSE 
ENDIF 

***// CALCULATE THE AHAT SUB I’S CWEIGHT ESTIMATES) 

SUMNAI = 0. 

DO 165 1=1, K 
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AHATI(I) = QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + NCI) * AHATI(I) 

165 CONTINUE 

IFCPRNT.EQ.1) THEN 
WRITEC1,0006) 

WRITE(1,C1C15) 

WRITECli3334) AHATI 

ELSE 
END IF 

***// CALCULATE 1 REPLICATION OF UPPR ALFA C. L. ON SYSTEM RELIABILITY 
DEGFR(LOOP) = 2 * (i + BIGF) 

DO 170 ALF=1, MAXALF 

CALL MDCHKl - ALFA(ALF) ,DEGFR(LOOP) ,CHISQ(ALF,LOOP), ERROR) 
QHTUPR(ALF,LOOP) = CHISQ(ALF,LOOP) / (2 * SUMNAI) 

IFCFLAG.EQ.1) THEN 

RHTSTR(ALF,LOOP) = 1 -(CHISQ(ALF,LOOP) / REAL(2*N(1))) 

ELSE 
END IF 

* + (ALF.LOOP), ALFA(ALF) 

***// CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR COMPNTS. IN. SERIES 

CALL RHTSRS(QHTUPR(ALF,LOOP), AHATI,K, RSHATCALF,LOOP)) 

* +T40,F8.5)') RSHATCALF.LOOP) 

CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE *** 

CALL RHTBRGC QHTUPRCALF,LOOP),AHATI,K,RSHTBRC ALF,LOOP)) 

170 CONTINUE 

***// this else and ENDIF are for the test against MAXRUN ***vt******* 
ELSE 

WRITE(1,'(" 'program exceeded THE max no. OF RUNS", 

+'' ALLOWED OF: " ,16)‘) TOTREP 
GOTO 9999 


END IF 
GOTO 10 
END IF 

WRITE(2,'( "UNSORTED RSHAT 1 IS: " ,/10(F8.5))' ) 
+(RSHAT(l,LOOP), LOOP=l, MAXREP) 

WRITE(2,'( "UNSORTED RSHAT 2 IS; " ,/10(F8. 5)>' ) 
+(RSHAT(2,LOOP), LOOP=l, MAXREP) 

IFCFLAG.EQ. 1) THEN 

WRITEC2,'C "UNSORTED RHTSTR 1 IS: " ,/10CF8. 5)) ’ ) 
+CRHTSTRCl,LOOP), LOOP=l, MAXREP) 

WRITEC2,'C UNSORTED RHTSTR 2 IS: " ,/10CF8.5))’) 
+CRHTSTRC2,L00i'), LOOP=l, MAXREP) 

ELSE 
END IF 

IFCK.EQ.5) THEN 

WRITEC2.'C "UNSORTED RSHTBR 1 IS: " ,/10CF8.5))') 
+CRSHTBRC1,LOOP), LOOP=l, MAXREP) 

WRITEC2,'C UNSORTED RSHTBR 2 IS; " ,/10CF8.5))') 
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+(RSHTBR(2,L00P), L00P=1, MAXREP) 

ELSE 
END IF 

'V**// SORT THE ARRAYS OF SYSTEM UNRELIABILITIESCl FOR EACH CONF. LEVEL) 


DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPRCALF,REPS) 
TRANSR(REPS) = RSHAT(ALF,REPS) 
TRNSTR(REPS) = RHTSTR(ALF,REPS) 
TRANBR(REPS) = 'RSHTBR(ALF,REPS) 
800 CONTINUE 

CALL SHS0RT(TRANSQ,KEY1,MAXREP) 

CALL SHS0RT(TRANSR,KEY2,MAXREP) 

CALL SHSORTCTPJs’STR,KEYS,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 
RSHAT(ALF,REPS) = TRANSR(REPS) 
RHTSTR(ALF,REPS) = TRNSTR(REPS) 
RSHTBR(ALF,REPS) = TRANBR(REPS) 
900 CONTINUE 

700 CONTINUE 


***// PRINT OUTPUT ^PORT HEADINGS 

WRITEC1,6666) 

WRITEC1,6667) MAXREP 
WRITEC1,6668) K 
WRITEC1,6669) 

IFCK.EQ.5) THEN 
WRITEC1,6699) 

ELSE 
ENT) IF 

WRITEC1,6670) MSTRQ 
WRITEC1,6671) 

WRITEC1,C1C15) 

WRITEC1,3334) AI 
WRITEC1,0007) 

WRITEC1,C1C15) 

WRITEC1,3334) QI 
WRITEC1,0005) 

WRITEC1,C1C15) 

WRITEC1,3335) N 
WRITEC1,6674) 

WRITEC2,6666) 

WRITEC2,6667) MAXREP 
WRITEC2,6668) K 
WRITEC2,6669) 

IFCK.EQ.5) THEN 
WRITEC1,6699) 

ELSE 
END IF 

WRITEC2,6670) MSTRQ 
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C WRITE(2,6671) 

C WRITE(2,C1C15) 

C WRITE(2,3334) AI 

C WRITE(2,0007) 

C WRITE(2,C1C15) 

C WRITE(2,3334) QI 

C WRITEC2,0005) 

C WRITE(2,C1C15) 

C WRITE(2,3335) N 

C WRITE(2,6674) 

C 

C WRITE(2,’("SORTED RSHAT 1 IS: " ./10(F8.5))’) 

C +(RSHAT(1,REPS), REPS=1, MAXREP) 

C WRITE(2,'( "SORTED RSHAT 2 IS: " ,/10(F8. 5)) ' ) 

C +(RSHAT(2,REPS), REPS=1, MAXREP) 

C IF(FLAG.EQ. l) THEN 

C WRITE(2,’C "SORTED RHTSTR 1 IS: " ,/10(F8. 5))’ ) 

C +(RHTSTR(1,REPS), REPS=I, MAXREP) 

C WRITE(2, ‘(‘'S0RTED RHTSTR 2 IS: " ,/10(F8. 5))‘) 

C +(RHTSTR(2,REPS), REPS-1, MAXREP) 

C ELSE 

C END IF 

C IF(K.EQ.5) THEN 

C WRITE(2,'( "SORTED RSHTBR 1 IS: " ,/10(F8.5))') 

C +(RSHTBR(1,REPS), REPS=1, MAXREP) 

C WRITE(2,’(‘'SORTED RSHTBR 2 IS: " ,/10(F8.5))') 

C +(RSHTBR(2,REPS), REPS=1, MAXREP) 

C ELSE 

C END IF 

***// COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. (SERIES SYSTEM) 
***// AND FOR THE 5-COMPONENT BRIDGE STRUCTURE 
CALL RSRS(QI,K,RS) 

WRITEC 1,'C" ",////"THE TRUE SERIES SYSTEM ' ', 

+"RELIABILITY VALUE IS: " ,T51,F8. 5)' ) RS 
CALL RBRIDG(QI,K,RSBRDG) 

IF(K.EQ.5) THEN 

WRITE(1,'(" ",////"THE TRUE BRIDGE STRUCTURE ", 

+"RELIABILITY VALUE IS: " ,T51,F8. 5)') RSBRDG 
ELSE 
END IF 

WRITEC1,6675) 

***// COMPUTE THE DIFFERENCE 'DELTAR' BIVN. RS AND RSHAT OF THE THEO 
***// RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY) 

IFCFLAG.EQ. 1) THEN 
WRITEC1,5755) 

ELSE 
END IF 

DO 450 ALF=1, MAXALF 

QUANTL(ALF) = MAXREP * (1 “ ALFA(ALF)) 

DELTAR(ALF) = RS - RSHAT(ALF,QUAI^L(ALF)) 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF)) 

IFCFLAG.EQ. 1) THEN 

DELSTR(ALI’) = RS - RHTSTR(ALf ,QUANTL(ALF)) 







WRITE(1,S555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
WRITEC1,5656) RHTSTR(ALF,QUANTL(ALF)) 

WRITEC1,5657) DELSTR(ALF^ 

ELSE 
END IF 

IF(K.EQ.5) THEN 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTLCALF^) 
WRITE(1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 
WRITE(1,5666) RSHTBR(ALF,QUANTL(ALF)) 

WRITEC1,5667) DELBRG(ALE) 

ELSE 
END IF 

WRITE(1,5555) MAXREP, ALFACALF), REAL(QUANTL(ALF)) 

WRITEC1,5556) RSHATCALF,QUANTLCALF))= 

WRITEC1,5557) DELTARCALF) 

450 CONTINUE 

PRINT*, 'QUANTLCl) IS:’, QUANTLCl) 

PRINT*, 'QUANTLC2) IS:', QUANTLC2) 

***// FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*** 
***// ******** RSHAT ******* //*** 

WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - RSHATCALF,REPS) 

CONTINUE 

DO 600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITECl.'C" ",/"TRUE CONFIDENCE LIMIT IS:’’, 
F8.4)') 

CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO to 620 

ELSEIFCDIFFCREPS). LT. 0.) THEN 
TRUQNTCALF) = REPS 
GO TO 610 

ELSE 
ENT) IF 
CONTINU'E 

IFCTRUQNTCALF).EQ.0. ) THEN 
WRITEC1,4443) ALFACALF) 

WRITEC1,'C" ”,/'’THE SMALLEST", 

" DIFFERENCE BETWEEN RS AND RSHAT IS: ",F10.5)') DIFFC 
MAXREP) 

ELSEIFCTRUQNTCALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFACALF) 

WRITEC 1,'C" ",/"ALL RSHAT", 

'' AR_ GREATER THAN RS" )' ) 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.ABSCDIFFCTRUQNTCALF) - 1))) 
THEN 

WRITEC1,4444) ALFACALF), 

CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)) 

WRITEC 1,4446) 

ELSE 


500 


+ 

+ 


600 

610 


+ 

+ 

+ 

+ 

+ 
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WRITEC1,4444) ALFA(ALF), 

+ ((TRUQNT(ALF)-l) / REAL(MAXREP)) * 100. 

WRITEC1,4445) RSHAT(ALF,TRUQNT(ALF)-1) 
WRITEC1,4447) 

620 END IF 

400 CONTINUE 


***// FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //*** 
***// irf ck -k ** * * rshTBR CBRIDGE) ******* //*** 


501 


601 

611 


621 

401 


IFCK.EQ.5) THEN 
DO 401 ALP=1,MAXALF 
TRUQNTCALF) = 0 

DO 501 REPS=1, MAXREP 

DIFFCREPS) = RSBRDG - RSHTBRCALF,REPS) 

CONTINUE 

DO 601 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C " " ,/ " TRUE CONFIDENCE LIMIT IS: " , 

+ F8.4)*) 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 621 

• ELSEIFCDIFFCREPS). LT. 0. ) THEN 
TRUQNTCALF) = REPS 
GO TO 611 

ELSE 
END IF 
CONTINUE 

IFCTRUQNTCALF). EQ. 0.) THEN 
WRITEC1,4443) ALFACALF) 

WRITEC1,'C" ",/*’THE SMALLEST”, 

+ '' DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS: ” , 

+ F10.5)’) DIFFCMAXREP) 

ELSEIFCTRUQNTCALF).EQ. 1.) THEN 
WRITEC1,4442) ALFACALF) 

WRITECli'C" ALL RSHTBR”, 

+ "ARE GREATER THAN RSBRDG'')') 

ELSEIFCABSCDIFFCTRUQNTCALF))).LE.ABSCDIFFCTRUQNTCALF) - 1))) 
+ THEN 

WRITEC1,4444) ALFACALF), 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

END IF 
CONTINUE 
ELSE 
END IF 
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***// FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE //**’> 

*******iV RHTSTR **rrtW^* //**iiV 


IFCFEAG.EQ.1) THEN 
DO 4400 ALF=1,MAXALF 
TRUQNT(ALF) =0 

DO 5500 REPS=1, MAXREP 

DIFF(REPS) = RS - PJITSTR(ALF,REPS) 

5500 CONTINUE 

DO 6600 REPS=1, MAXREP 

IF(ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT(ALF) = REPS 

WRITE(1,'(‘’ ”,/'‘TRUE CONFIDENCE LIMIT IS:'*, 

+ F8.4)') 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

GO TO 6620 

ELSEIF(DIFF(REPS). LT. 0. ) THEN 
TRUQ.NT(ALF) = REPS 
GO TO 6610 

ELSE 
END IF 

6600 CONTINUE 

6610 IF(TRUQNT(ALF).EQ.O. ) THEN 

WRITEC1,4443) ALFACALF) 

WRITE(1,'(" ’',/"THE SMALLEST", 

+ '' DIFFERENCE BETWEEN RS AND RHTSTR IS: " , 

+ F9..5)') DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ.1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,'(" ",/"ALL RHTSTR", 

+ "ARE GREATER THAN RS" )' ) 

ELS£IF(ABS(DIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITEC1,4444) ALFACALF), 

+ CTRUQNTCALF) / REALCMAXREP)) * 100. 

WRITEC1,4448) RHTSTRCALF,TRUQNTCALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) /REALCMAXREP)) * 100. 

WRITEC1,4448) RHTSTRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***// PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION ***** 

IFCPRNT.EQ.1) THEN 
1 = 1 

185 WRITECl.REPSHD) ALFACSELCTA). ALFACSELCTA), 

+ALFACSELCTB),ALFACSELCTB) ,aLFACSELCTA) ,ALFACSELCTA).ALFACSELCTB), 
+ALFACSELCTB) 

175 IFCI.GE.CMAXREP + 1)) THEN 
GOTO 180 
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ELSE 


IF( (I. EQ. 71). OR. (I.EQ. 211). OR. (I. EQ; 351). OR. (I. EQ. 491). OR. 

+ (I.EQ.631).OR.(I.EQ.771).OR.(I.EQ.911).OR.(I.EQ.1051) ) THEN 

I = I + 70 
WRITEd,'("+”)• ) 

GOTO 185 

ELSE 

WRITE(1,3336) I, INT(DEGFR(I)), CHISQ(1,T), QHTUPR(1,I), 

+ CHISQ(2,I), QHTUPR(2,I) 

END IF 

IFCI + 70. LE. MAXREP) THEN 

WRITEC1,3337) 1+70,INT(DEGFR(1+70)),CHISQ(1,1+70), 

+ QHTUPRC1,1+70),CHISQ(2,I+70),QHTUPR(2,I+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 

9999-WRITE(1,'( "THE TOTAL NO OF REPS WAS:", 18)') TOTREP 

WRITE(1,'("THE TOTAL NO OF EFFECTIVE REPS WAS: ",I8)’) LOOP 
WRITE(1,‘("THE TOTAL NO OF NO FAILURE RUNS WAS:’^I8)') ZFAILS 
WRITEC1,'(''AVERAGE NO. OF COMPONENTS PER REPLICATION WITH '', 

+''NO FAILURES:'',F5.2)') ZFPREP / MAXREP 

WRITEC 1,'( "THE TOTAL NO OF RUNS WITH FAILURES WAS:",18)') FAILS 
0008 FORMAT (/ 3X,'C l',5X 'C 2', 

+5X,'C 3',5X,'C 4',5X,'C 5',5X,'C 6',5X,'C 7',5X, 

+'C 8',5X,'C 9',5X,'C 10',4X,'C 11’,4X, 

+'C 12',4X,'C 13',4X,'C 14',4X,'C 15') 

0009 F0RMAT(/1X,'REP NO',2X,'DF',1X 'CHISQR(’,F4. 3,’)',1X, 

+'QHTUPRC’,F4.3,')’,1X,'CHISQR(‘,F4.3,')',1X,'QHTUPR(',F4. 3,’)’, 
+2X,'REP NO' ,2X,'DF^1X,'CHISQR(^F4.3 ')\1X, 

+'QHTUPRC',F4. 3,')',IX,'CHISQRC',F4.3,')’,IX,'QHTUPRC',F4.3,')'/) 
0001 FORMAT C///’UNIFORM RANDOM DEVIATES ARE:') 

0002 FORMAT C///'BERNOULLI TRIALS ARE:') 

0003 FORMAT C///'TOTAL NO. OF FAILURES FOR EACH COMPONENT:') 

0004 FORMAT C///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT:') 

0005 FORMAT C///'TOTAL NUMBER OF MISSION TESTS: ') 

0006 FORMAT C///’ESTIMATED WEIGHTS FOR EACH COMPONENT:') 

0007 FORMAT C///'Q I FOR EACH COMPONENT: ') 

1111 FORMAT C15F8.5) 


2222 FORMAT 

3333 FORMAT 

3334 FORMAT 

3335 FORMAT 

3336 FORMAT 

3337 FORMAT 

4442 FORMAT 
+’LIMIT 

4443 FORMAT 
+'LIMIT 

4444 FORMAT 
+'LIMIT 


C/1X,15CI4,4X)) 

C/1X,15CI4,4X)) 

C/15F8.5) 

C/1X,15CI4,4X)) 

CT3.I4,T9,I3,T13,F11.5,T27,F8.5,T39,F11.5,T53,F8.5) 
C’+‘,T67.I4,T73,I3,T77,F11.5.T91.F8.5,T103,F11.5,T117,F8.5) 
C' ',///'THE resulting Cl - ,F4.3,’) CONFIDENCE ’, 

IS:',T50,’ OO.OOO ') 

C' ’,///^THE RESULTING Cl " ',F4.3,’) CONFIDENCE ', 

IS:',T50,'100.0000') 

C ',///'THE RESULTING Cl - ',F4.3,’) CONFIDENCE ', 

IS: ',T50,F8. 4) 
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4445 FORMAT (' 

4446 FORMAT (' 

4447 FORMAT (' 

4448 FORMAT (' 

4449 FORMAT (’ 

5555 FORMAT (' 

5556 FORMAT (* 

5557 FORMAT (' 

5656 FORMAT (' 

5666 FORMAT (' 

5657 FORMAT C‘ 

5667 FORMAT (' 
5755 FORMAT (’ 


'RHTSTR' 
6666 FORMAT ('+ 
+’ 


,/'THE RSHAT VALUE CLOSEST TO RS IS: ',T51,F8.5) 
./'(FIRST NEGATIVE DIFFERENCE)') 

./'(ELEMENT PRECEEDING FIRST NEGATIVE DIFFERENCE)') 
./'THE RHTSTR VALUE CLOSEST TO RS IS: '.T51.F8.5) 

./'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: '.T51.F8.5) 
.///'THE '.I4.'(l-'.F4.3.') QUANTILE IS: ',T49.F8.3) 
./'THE VALUE OF RSHAT FOR THAT QUANTILE IS: '.T51.F8. 5) 
./'THE DIFFERENCE(RS - RSHAT) IS: '.T51.F8. 5) 

./'THE VALUE OF RHTSTR FOR THAT QUANTILE IS: ',T51.F8.5) 
./'THE VALUE OF RSHTBR FOR THAT QUANTILE IS: '.T51.F8. 5) 
./'THE DIFFERENCE(RS - RHTSTR) IS:'.T51.F8.5) 

./'THE DIFFERENCE(RS - RSHTBR) IS: '.T51.F8.5) 

.///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR'. 


+’ ALL COMPONENTS THE CLOSED FORM: SERIES SYSTEM RELIABILITY 


IS COMPUTED') 


* RUN INPUT SETTINGS 

***tri-**' ) 

.//'NUMBER OF REPLICATIONS:',T50.I4) 

.//’NUMBER OF COMPONENTS: '.T50.14) 

.//’SYSTEM RELIABILITY FUNCTION: ',T50.'SERIES') 
.//'SYSTEM RELIABILITY FUNCTION: '.T50.'BRIDGE') 
.//'MASTER UNRELIABILITY USED:',T50.F8.5) 


6667 FORMAT (' 

6668 FORMAT (' 

6669 FORMAT (' 

6699 FORMAT (' 

6670 FORMAT (' 

6671 FORMAT (' './'/'INPUT WEIGHTS(A SUB l”S):') 

6674 FORMAT ('+'.///' 

+'**R UNRESULT 
4.' **5V*A#r*******A' ) 

6675 FORMAT ('+'j///'************************************;^***********' 
+' ESTIMATE ERRORS **********************************************^’ 

+'**************'j 

6676 FORMAT ( * -f'. / / /' ***»v********************** * ********************' 

+' TRUE CONFIDENCE LIMITS 

+'**************') 

END 
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APPENDIX H. FORTRAN CODE FOR ALTERNATE PROCEDURE B 
FOR ESTIMATING THE LOWER CONFIDENCE LIMIT FOR SELECTED 

SYSTEMS 


PROGRAM ZFYSCN 

* * 


* TITLE: BINOMIAL INTERVAL ESTIMATION PROCEDURE * 

* ZERO FAILURES ALLOWED; NO SCALING * 

* AUTHOR: E. F. BELLINI, LT, USN * 

* MODIFIED BY: LT.VALERIE A. COVINGTON,USN (MAR 90) 

* DATE: NOV 09 * 

* * 

* THIS PROGRAM COMPUTES THE TRUE CONFIDENCE LEVEL FOR THE ESTIMATE * 

* RELIABILITY OF A SERIES AND BRIDGE SYSTEM GIVEN THE RELIABILITY * 

* OF THEIR COMPONENTS * 

* * 


* IN ITS PRESENT CONFIGURATION THIS PROGRAM IS SET UP TO RUN 12 * 

* TIMES EACH TIME PRODUCING 1000 REPLICATIONS USING A DIFFERENT * 

* SET OF INPUT DATA. RUN THE PROGRAM FROM CMS BY TYPING 'B1 EXEC'.* 

* THE REXX EXEC PROGRAM * 

* 'Bl' CALLS THE INPUT FILES TO BE READ AND NAMES THE 12 OUTPUT * 

* FILES RESULTING FROM THE 12 CONSECUTIVE RUNS. BY EDITING THE * 

* INDEX COUNTERS I, J, K OF THE 'Bl' EXEC ONE CAN RUN ANY USER- * 

* SPECIFIC RUN FROM JUST ONE RUN TO ALL 12. * 

* ■ it 

* VARIABLES USED * 

* * 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


AHATI 

AI 

ALFA 

BIGF 

CHISQ 

C1C15 

DEGFR 

DELBRG 

DELSTR 

DELTAR 

DIFF 

EPS 

ERROR 

FAILS 

FI 

FLAG 

INC 

KEYl 

KEY2 

KEYS 

KEY4 

KK 

LOOP 

MAXALF 


WEIGHT ESTIMATES FOR EACH COMPONENT * 
INPUT WEIGHTS FOR EACH COMPONENT * 
LEVELS OF SIGNIFICANCE * 
TOTAL NO. OF FAILURES FOR EACH REPLICATION * 
CHI-SQUARE RANDOM VARIABLE VALUE * 
FORMAT LABEL * 
DEGREES OF FREEDOM * 
DIFFERENCE FOR BRIDGE SYSTEM * 
DIFFERENCE FOR SERIES SYSTEM- CLOSED FORM * 
DIFFERENCE FOR SERIES SYSTEM * 
DIFFERENCE (TRUE REL. - ESTIMATED REL.) * 
SMALL QUANTITY(CONSTANT) * 
PARAMETER FOR IMSL ROUTINE * 
COUNTS NO. OF REPLICATIONS WITH AT LST. 1 FAILURE * 
NO. OF FAILURES FOR EACH COMPONENT(ALL MISSION TST)* 
1 IF ALL COMP. HAVE SAME NO. OF MISSION TESTS * 
INCREMENT STEP SIZE FOR ROUTINE USMNMX * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 
ARRAY OF INDECES FOR ROUTINE SHSORT * 


ARRAY SIZING PARAMETER FOR THE MAX NO OF COMPONENTS* 
COUNTS NO. OF REPLICATION PERFORMED * 
MAX NO. OF SIGNIFICANT LEVELS DESIRED(ARRAY SIZING)* 
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■k 

MAXREP : 

MAX NO. OF REPLICATIONS 

k 

* 

MAXRUN : 

MAX NO. OF PROGRAM ITERATIONS ALLOWED 

k 

* 

MSTRQ : 

MASTER UNRELIABILITYCUSED WITH Al'S TO CALC. Ql’S) 

k 

* 

MULT : 

MULTIPLIER FOR RANDOM NO. GENERATOR SRND 

k 

* 

N : 

NO. OF MISSION TEST FOR EACH COMPONENT * 

* 

NIMAX : 

MAX NO. OF MISSION TESTS 

k 

k 

NIMIM : 

MIN NO. OF MISSION TESTS 

k 

k 

NINDX : 

INDEX NO. OF MAX NO. OF MISSION TESTS 

k 

k 

NIREAL : 

NO. OF MISSION TESTS TRANSFORMED TO REAL 

k 

k 

NMAX : 

MAX NO. OF MISSION TESTS FOR OUTPUT CONTROL 

k 

k 

NPRNT : 

FLAG FOR DETAILED REPORT OUTPUT 

k 

k 

PRNT : 

SAME AS ABOVE(PARAMETER) 

k 

k 

QHATI : 

UNRELIABILITY ESTIMATES FOR EACH COMPONENT 

k 

k 

QHTMAX : 

LARGEST QHATI 

k 

k 

QHTUPR : 

UPPER LIMIT ON SYSTEM UNRELIABILITY 

k 

k 

QI : 

INPUT UNRELIABILIY FOR EACH COMPONENT 

k 

* 

QINDX : 

INDEX 

k 

* 

QUANTL : 

QUANTILE 

k 

* 

REPSHD : 

REPLICATIONS HEADING FORMAT NUMBER 

k 

k 

RHTSTR : 

SERIES SYSTEM RELIABILITY ESTIMATE(CLOSED FORM) 

k 

k 

RS : 

TRUE SERIES SYSTEM RELIABILITY 

k 

k 

RSBRDG : 

TRUE BRIDGE SYSTEM P.ELIABILITY 

k 

k 

RSHAT : 

SERIES SYSTEM RELIABILITY ESTIMATE 

k 

k 

RSHTBR : 

BRIDGE SYSTEM RELIABILITY ESTIMATE 

k 

* 

SEED : 

PARAMETER 

k 

k 

SELCTA : 

SIGNIFICANCE LEVEL SELECTION 

k 

k 

SELCTB : 

SIGNIFICANCE LEVEL SELECTION 

k 

k 

SORT : 

PARAMETER FOR ROUTINE SRND 

k 

k 

SUMNAI : 

SUM OF THE PRODUCT OF Nl'S AND Al'S 

k 

k 

TEMP : 

TEMPORARY ARRAY 

k 

k 

TOTREP : 

TOTAL NUMBER OF PROGRAM ITERATIONS 

k 

k 

TRANBR : 

TEMPORARY ARRAY 

k 

k 

TRANSQ : 

TEMPORARY ARRAY 

k 

k 

TRANSR : 

TEMPORARY ARRAY 

k 

k 

TRIALS : 

BERNOULLI TRIALS ARRAY (2-DIM) 

k 

k 

TRNSTR : 

TEMPORARY ARRAY 

k 

k 

TRUQNT : 

TRUE QUANTILE 

k 

k 

UNIRV : 

UNIFORM RANDOM DEVIATES (2-DIM) 

k 

k 

ZFAILS : 

TOTAL NUMBER OF REPLICATIONS WITH ZERO FAILURES 

k 

k 

k 

ZFPREP : 

NO. OF COMPNTS. WITH ZERO FAILURES PER REPLICATION 

k 

k 


PARAMETER (KK=10,MAXALF=2,NPRNT=0) 

PARAMETER (MAXREP=1000, MAXRUN=2000, EPS=. 000001) 

REAL*4 UNIRVC15,1000),TEMP(1000),QI(KK),AI(KK),AHATI(KK) 

REAL*4 QHATI(KK), UMAX, NNMAX, QHIMAX, CHISQR(S,5), ALFA(MAXALF) 
REAL*4 DF(5),AALFAC 5),SUMNAI,RSHAT(MAXALF,MAXREP),RS 
REAL*4 KEYKMAXREP) ,KEY2(MAXREP) ,KEY3(MAXREP) ,TRNSTR(MAXREP) 
REAL*4 DEGFRCMAXREP), QHTUPR(MAXALF,MAXREP),CHISQCMAXALF,MAXREP) 
REAL*4 QUPAl(MAXREP), QUPA2(MAXREP),RHTSTR(MAXALF,MAXREP) 

REAL*4 DELTARCMAXALF), TRANSQ(MAXREP),TRANSR(MAXREP),DIFF(MAXREP) 
REAL*4 DELSTR(MAXALF),NIMIN,NIMAX,NIREAL(KK) 

REAL*4 RSHTBR(MAXALF,MAXREPJ,DELBRG(MAXALF),KEY4(MAXREP) 

REAL*4 TRANBRCMAXREP), RSBRDG ,MSTRQ 
REAL*4 ZFPREP 
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REAL*4 RSHATI(KK),SI(KK)=,QMTIU(KK) 

REAL*4 MXQHAT,RI(KK),SUMRN 

REAL*4 QHATMU(MAXALF),CHIVAL(MAXALF),P(KK),PTEMP(10) 

INTEGER SEED, MULT, SORT, TRIALS(15,1000), BIGF, FI(KK), N(KK) 
INTEGER NINDX, QINDX, ERROR, REPS, SELCTA, SELCTK, TOTREP 
INTEGER C1C15, REPSHD, SELCTB, ALF, FLAG, LOOP,PRNT,HFI 
II^GER QUANTL(MAXALF) , TRUQNT(MA^LF) ,ZFAILS,FAILS,INC 

CHARACTER*8 LOOPSO(MAXREP) 

DATA SEED/123457/, MULT/1/, INC/1/ 

DATA AALFA/. 01,.0S,.9,.95,.99/, DF/1,5,10,30,40/ 

DATA ALFA/. 20,.050/ 

DATA SORT/0/ 

ASSIGN 8 TO C1C15 
ASSIGN 9 TO REPSHD 

* CALL COMPRS 
PRNT = NPRNT 


DO 12 1=1,KK 

AI(I) = 9999. 
N(I) = 99999999 
12 CONTINUE 

READ(03,*)K,MSTRQ 

DO 11 1=1,K 
READ(03,*) AI(I),N(I) 
11 CONTINUE 


IFCK.NE.5) THEN 

WRITE(1,’("WARNING: BRIDGE STRUCTURE ", 
+"ONLY USES THE FIRST 5 COMPONENTS")’) 

ELSE 
END IF 


***initIALIZE the QHTUPR array of UNRELIABILITY REPLICATIONS, *** 

* RSHAT ARRAY OF ESTIMATE SYSTEM RELIABILITY REPLICATIONS * 

* AND RHTSTR ARRAY OF EST. SYST. REL. FOR A SERIES SYST WHEN * 

***ALL THE COMPONENT MISSION TESTS ARE EQUAL IN NTJMBER *** 


DO 172 ALF=1,MAXALF 

DO 173 REPS=1,MAXREP 

QHTUPRCALF,REPS) = 0. 

RSHATCALF,REPS) = 0. 

RHTSTRCALF,REPS) = 0. 

RSHTBRCALF,REPS) = 0. 

LOOPSO(REPS)=' ****** • 

173 CONTINUE 

172 CONTINUE 

***SET FLAG TO 1 IF ALL COMPONENTS HAVE SAME NO. OF MISSION TESTS'W^^ 
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FLAG=1 

DO 50 1=1,K -1 

IF((N(I) - ^?(I+1^).NE.0) THEN 
FLAG=0 

ELSE 
END IF 
50 CONTINUE 

PRINT 'FLAG IS:', FLAG 

***MAIN PROGRAM OUTER LOOP START(EVERY LOOP IS ONE REPLICATION)*** 

ZFPREP = 0. 

ZFAILS =0 
FAILS = 0 
TOTREP = 0 
LOOP = 0 

10 IF(LOOP.LT.MAXREP) THEN 
LOOP = LOOP + 1 
IFCTOTREP. lt, MAXRUN) THEN 
TOTREP = TOTREP + 1 

SELCTA = 1 
SELCTB = 2 

***FILL ARRAY KEY(REPS) WITH INTEGERS 1 TO K TO BE USED AS OUTPUT*** 
***0F THE SUBROUTINE SHSORT *** 

DO 95 REPS=1, MAXREP 
KEYl(REPS) = REPS 
KEY2(REPS) = REPS 
KEY3(REPS) = REPS 
KEY4(REPS) = REPS 
95 CONTINUE 

***CALCULATE NMAX NOT TO PRINT LONGER THAN THE MAX SAMPLE SIZE*^* 
***CALCULATE THE MAXIMUM NO. OF TRIALS AND ITS INDEX NQ.*’'^* 

CALL IMAX(N,K,NMAX,NINDX) 


***CALCULATE THE QI'S FROM THE GIVEN MASTER Q AND THE Al'S*** 

DO 115 1=1, K 

QIC I) = MSTRQ * AI(I) 

115 CONTINUE 

DO 120 1=1,15 

DO 125 J=l,500 

UNIRV(I,J) = 999. 

TRIALS(I,J) = 99999 
125 CONTINUE 

120 CONTINUE 

*** DRAW UNIFORM (0,1) RV’S AND CONVERT TO BERNOULLI TRIALS*** 
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DO 130 1=1, K 

CALL SRND(SEED, TEMPi N(I), MULT, SORT) 
DO 135 J=l, N(;I) 

UNIRV(I,J) = TEMP(J) 

IF (UNIRV(I,J).LE. 1 - QICD) THEN 
TRIALS(I,J) = 0 

ELSE 

TRIALS(I,J) = 1 
END IF 

135 CONTINUE 
130 CONTINUE 


***CALCULATE THE NO. OF FAILURES FOR EACH COMPONENT*** 

DO 150 1=1, K 
FI(I) = 0 
150 CONTINUE 
lONECT = 0 

***CALCULATE THE F SUB I * S AND THE GRAND TOTAL NO. OF FAILURES*** 
BIGF = 0 
DO 155 1=1, K 

DO 160 J=l, N(I) 

;FI(I) = FI(I) + TRIALS(I,J) 

160 CONTINUE 
IF(FI(I).EQ.O) THEN 

ZFPREP = ZFPREP + 1 

ELSE 
END IF 

***CALCULATE THE QHAT SUB I'S: F SUB I’S DIVIDED BY N SUB I'S*** 
QHATI(I) = REAL(FI(I)) / N(I) 

BIGF = BIGF + FI(I) 

155 CONTINUE 

IF (FI(2) .EQ. 0) THEN 
ZFPREP=ZFPREP - 1 

ELSE 

ENDIF 

BIGF=BIGF-FI(2) 

FI(2) = 0 
DO 161 I=2,N(2) 

HFI=0 

DO 163 J=l,2 

CALL SRNDCSEED, PTEMPCJ), N(2), MULT, SORT) 

IF (PTEMPCJ) .GT. 1-QI(2)**.5) THEN 
HFI= HFI+1 
ENDIF 

163 CONTINUE 

IF (HFI .EQ. 2) THEN 
FIC2) = FI(2) +1 
ENDiF 

161 CONTINUE 

QHATI(2) = REAL(FI(2)) / N(2) 

BIGF= BIGF + FI(2) 

IF (FI(2) .EQ. 0) THEN 
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ZFPREP = ZFPREP + 1 
ENDIF 


***00^18 NUMBER OF COMPONENTS THAT HAV’ -AILED*** 
DO 136 1=1,K 

IF (FI(I) .NE. 0) IONECT=IOKV. 

136 CONTINUE 

***CASE WPRE NO COMPONENTS HAVE ANY FAI. ORES*** 


IF(BIGF.EQ.O) THEN 

LOOPSO(LOOP)=‘ *ZERO* ‘ 

*** SERIES ESTIMATE MODIFICATION (NO. OF FAILUBlES IRRELEVANT)*** 
DO 152 1=1,K 

SI(I) = REAL(N(I))-REAi,iFI(T)) 
RSHATI(I)=SI(I)/(SI(I)+(REAL(FI(I))+1. ) 

+ * FIN(.50,2.*(REAL(FI(I))+1. ),?*SI(I))) 

QHATIU(I)=1. -RSHATI(I) 

152 CONTINUE 

MXQHAT=0. 

DO 154 1=1,K 

IF (QHATIU(I) .GT. MXQHAT) THEN 
MXQHAT= QHATIU(I) 

ENDIF 

154 CONTINUE 

DO 156 ALF=1,MAXALF 

CALL MDCHI(ALFA(ALF),2*(1.+REAL(BIGF)),CHIVAL(ALF),ERROR) 
SUMRN=0. 

DO 157 1=1,K 

RI(I)=QHATIU(I)/MXQKAT 
SUMRN=SUMRN+SI(I)*REAL(N(I)) 

157 CONTINUE 

QHATMUC ALF)=CHIVALC ALF)/(2*SUMRN) 

156 CONTINUE 

DO 158 ALF=1,MAXALF 
RSHAT(ALF,L00P)=1. 

DO 159 1=1,K 

RSHATCALI'',LOOP)=RSHAT(ALF,LOOP)*( 1.-(RI(I)*QHATMUCALF))) 
IFCFLAG. EQ. 1) THEN 

RHTSTR(ALF,LOOP) = RSHAT(ALF,LOOP) 

ELSE 
END IF 
159 CONTINUE 

158 CONTINTJE 

***CALCULATE VALUE OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTUREt*^ 

DO 141 ALF=1, MAXALF 
DO 142 1=1,K 
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P( I')=1-RI ( I)*QHATMU( ALF) 

142 CONTINUE 

RSHTBR(ALF,LOOP)=P(l)*p(4)+p(2)*p(5)+P(1)*P(3)* 
C P(fS)+P( 2)*P( 3)*PC4) -P( 1)*P( 2)*Pf^3 )*P(4) -P( 1)* 

C PC2)*P(3)*P(5)-P(l)*P(2J*P(4)*Pt5)-P(l)*P(3)* 

C P(4)*P(5D“PC2)*P(3)*P(4)*P(5)+2*P(1)*P(2)* 

C PC3)*P(4)*F(5) 

141 CONTINUE 

ZFAILS = ZFAILS + 1 
DEGFRCLOOP'l = 2. 

GO TO 10 

els; 

FAILS = FAILS + 1 
END I? 


*^FIND THE MAX OF THE INDIVIDUAL COMPONENT UNRELIABILITIES--'** 


CALL RMAXCQHATI, K, QHTMAX, QINDX) 


•***CALCULATE THE AHAT SUB I'S (WEIGHT ESTIMATES)*** 

SUMNAI = 0. 

DO 165 1=1, K 

AHATI(I) ■■= QHATI(I) / QHTMAX 
SUMNAI = SUMNAI + NCI) * AHATI(I) 

165 CONTINUE 

***1 FAILURE ONLY SERIES SYST. 

IF (lONECT .EQ. 1) THEN 
LOOPSO(LOO?)=' *ONE* ' 

ENDIF 

***CALCULATE 1 REPLICATION OF UPPR ALFA C.L, ON SYSTEM RELIABILITY*** 
iJEGFR(LOOP) = 2 * (1 + BIGF) 

***SERI';S ESTIMATE MODIFICATION (NO. OF FAILURES IRRELEVANT)*''* 

DO 162 1=1,K 

SI(I) = REAL(N(I?)-RSA?.:CFI(I)) 

RSHATKI)=SI(I)/(Sl(D-fCREALCFKI))+l. ) 

+ * F.‘N(.50,2.*(REAL(FI(I))+1. ),2*SI(I))) 

QHATIU(I)=1. -RSHATKI) 

162 CONTINUE 

MXQHAT^O. 

DO 164 1=1,K 

IF (QHATIU(I) .GT. MXQHAT) THEN 
MXQHAT = QHATIUCI) 

ENDIF 

164 CONTINUE 
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DO 166 ALF=1,MAXALF 

CALL- MDCHI(ALFA(ALF) ,2*( 1. +REAL(BIGF) )=,CHIVAL(ALF):,ERROR) 
SUMRN=0. 

DO 167 1=1,K 

RI(I)=QHATIU(I)/MXQHAT 
SUMRN=SUMRN+RI(I)*REAL(N(I)) 

167 CONTINUE 

QHATMUCALF)=CHIVALCALF)/(2*SUHRN) 

166 CONTINUE 

DO 168 ALI’=1,MAXALF 
RSHATCALF,LOOP)=l. 

DO 169 1=1,K 

RSHAT(ALF,LOOP)=RSHAT(ALF,LOOP)*(1. -CRI(I)*QHATMUCALF))) 
IFCFLAG. EQ. 1) THEN 

RHTSTRCALF,LOOP) = RSHAT(ALF,LOOP) 

ELSE 
END IF 
169 CONTINUE 

168 CONTINUE 


***CALCULATE VALU’E OF THE SYSTEM RELIABILITY FOR BRIDGE STRUCTURE*** 

DO 170 ALF=1, MAXALF 
DO 206 1=1,K 

P(T)=1-RI(I)*QHATMU(ALF) 

206 CONTINUE 

RSHTBRCALF,LOOP)=P(1)*P(4)+P(2)*P(5)*vF(1)*P(3)* 

C P(5)+P(2)*P(3)*P(4)-P(1)*P(2)*P(3)*P(4)>P(1)* 

C PC2)*P(3)*P(5)-P(1)*P(2)*P(4)*P(5)-P(1)*P(3)* 

C P(4)*P(5)-P(2)*P(3)*P(4)*P(5)+2*P(1)*P(2)* 

C P(3)*P(4)*P(S) 

170 CONTINUE 


***THIS ELSE AND ENDIF ARE FOR THE TEST AGAINST MAXRUN*** 

ELSE 

WRITECl,'('= " 'PROGRAM EXCEEDED THE MAX NO. OF RUNS'', 
+'' ALLOWED OF: " ,16)') TOTREP 
GOTO 9999 

END IF 
GOTO 10 
END IF 

WRITE(2,'( "UNSORTED RSKAT 1 IS: " ,/10(F8. ■:))') 
+(RSHAT(1.LOOP), L00P=1, MAXREP) 

WRITE(2,‘( "UNSORTED RSHAT 2 IS: " ,/10(F8.5)) *) 
■f-CRSHAT(2,L00P), L00P=1, MAXREP) 

IFCFLAG. EQ.l) THEN 

WRITE(2,'C "UNSORTED RHTSTR 1 IS: " ,/10(F8.5))') 
+(RHTSTR(1,M0P), L00P=1, MAXREP) 

WRITE(2,’(‘'UNSORTED RHTSTR 2 IS: " ,/10(F8.5))’) 

+(RHTSTRC2,LOOP), L00P=1, MAXREP) 
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ELSE 
END IF 

lECK.EQ.S) THEN 

WRITE(2,•(''UNSORTED RSHTBR 1 IS: " ,/10(F8. 5)) ) 
+(RSHTBR(l,LOOP),_ L00P=1, MAXREP) 

WRITE(2,' ("UNSORTED RSHTBR 2 IS: " ,/10(F8. 5)) ) 
+(RSHTBR(2,L00P:), LOOP=l, MAXREP) 

ELSE 

END IF ... 

WRITE (2,’("ZERO AND ONE FAILURE REPS: " ,/10(A8)0 ) 

+ (LOOPSO(LOOP),LOOP=l,MAXREP) 

***SORT THE ARRAYS OF SYSTEM UNRELIABILITIES(1 FOR EACH CONF. LEVEL)*** 

DO 700 ALF=1, MAXALF 

DO 800 REPS=1, MAXREP 

TRANSQ(REPS) = QHTUPR(ALF,REPS) 

TRANSR(REPS) = RSHAT(ALF,REPS) 

TRNSTR(REPS) = RHTSTR(ALF,,REPS) 

TRANBR(REPS) = RSHTBRULF,REPS) 

800 CONTINUE 

CALL SHS0RT(TRANSQ,KEY1,MAXREP) 

CALL SHS0RT(TRANSR,KEY2,MAXREP) 

CALL SHSORT(TRNSTR,KEYS,MAXREP) 

CALL SHS0RT(TRANBR,KEY4,MAXREP) 

DO 900 REPS=1, MAXREP 

QHTUPR(ALF,REPS) = TRANSQ(REPS) 

RSHAT(ALF,REPS) = TRANSR(KEPS) 

RHTSTR(ALF,REPS) = TRNSTR(REPS) 

. RSHTBR(ALF,REPS) = TRANBR(REPS) 

900 CONTINUE 

700 CONTINUE 

***PRINT OUTPUT REPORT HEADINGS*** 

WRITE(1,6666) 

WRITE(1,6667) MAXREP 
WRITE(1,6668) K 
WRITE(1,6669) 

IF(K.EQ.5) THEN 
WRITE(1,6699) 

ELSE 
END IF 

WRITE(1,6670) MSTRQ 
WRITE(1,6671) 

WRITE(1,C1C15) 

WRITE(1,3334) AI 
WRITE(1,0007) 

WRITE(1,C1C15) 

WRITE(1,3334) QI 
WRITE(1,0005) 

WRITE(1,C1C15) 

WRITE(1,3335) N 
WRITE(1,6674) 
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***COMPUTE THE VALUE RS OF THE TRUE SYSTEM REL. FNCTN. (SERIES SYSTEM)*** 
*** AND FOR THE 5-COMPONENT BRIDGE STRUCTURE*** 

CALL RSRS(QI,K,RS) 

WRITE(1,'(" '^////"THE TRUE SERIES SYSTEM ”, 

+”RELIABILITY VALUE IS: ” ,T51,F8.5)') RS 
CALL RBRIDG(QI,K,RSBRDG) 

IF(K.EQ.5) THEN 

WRITE(1,'(” ”,////”THE TRUE BRIDGE STRUCTURE ”, 

+”RELIABILITY VALUE IS: ” ,T51,F8.5) ’ ) RSBRDG 
ELSE 
END IF 

WRITE(:1,6675) 

***COMPUTE THE DIFFERENCE 'DELTAR’ BTVN. RS AND RSHAT OF THE THEO*** 
***RETICAL QUANTILE GIVEN BY ALFA(MUST USE SORTED RSHAT ARRAY)*** 

IFCFLAG. EQ, T) THEN 
WRITEC1,5755) 

ELSE 
END IF 

DO 450 ALF=1, MAXALF 

QUANTL(ALF) = MAXREP * (1 - ALFA(ALF)) 

DELTAR(ALF) = RS - RSHAT(ALF,QUANTL(ALF)) 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTL(ALF)) 

IFCFLAG.EQ.1) THEN 

DELSTR(ALF) = RS - RHTSTR(ALF,QUANTL(ALF)) 

WRITEC1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 

WRITEC1,5656) RHTSTR(ALF,QUANTL(ALF)) 

WRITEC1,5657) DELSTR(ALF) 

ELSE 
END IF 

IF(K.EQ.5) THEN 

DELBRG(ALF) = RSBRDG - RSHTBR(ALF,QUANTLCALF)) 

WRITEC1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 

WRITEC1,5666) RSHTBR(ALF,QUANTL(ALF)) 

WRITEC1,5667) DELBRG(ALF) 

ELSE 
END IF 

WRITEC1,5555) MAXREP, ALFA(ALF), REAL(QUANTL(ALF)) 

WRITEC1,5556) RSHAT(ALF,QUANTL(ALF)) 

WRITEC1,5557) DELTAR(ALF) 

450 CONTINUE 

PRINT*, ’QUANTL(l) IS:’, QUANTL(l) 

PRINT*, ’QUANTLC2) IS:', QUANTL(2) 

***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 

*** ■ A - A - AA rfr A - A T fr RSHAT ******* *** 

WRITEC1,6676) 

DO 400 ALF=1,MAXALF 
TRUQNT(ALF) = 0 

DO 500 REPS=1, MAXREP 

DIFFCREPS) = RS - RSHAT(ALF,P£PS) 

500 CONTINTJE 

DO 600 REPS=i, MAXREP 
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IF(ABS(DIFF(REPS)).LE.EPS) THEN 
TRUQNT(ALF) = REPS 

WRITE(1,‘(" " 'TRUE CONFIDENCE LIMIT IS:'', 

+ F8.4)') 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

GO to 620 

ELSEIFCDIFF(REPS). LT. 0. ) THEN 
TRUQNT(ALF) = REPS 
GO TO 610 

ELSE 
END IF 

600 CONTINUE 

610 IF(TRUQNT(ALF).EQ.O.) THEN 

WRITE0;i,4443) ALFA(ALF) 

WRITE(1,'("",/’'THE SMALLEST", 

+ ’’ DIFFERENCE BEIVEEN RS AND RSHAT IS:'',F10.5)') DIFFC 

+ MAXREP) 

ELSEIF(TRUQNT(ALF).EQ. 1.) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,'(" ",/"ALL RSHAT", 

+ "ARE GREATER THAN RS*‘)’) 

ELSEIF(ABS(DIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITEC1,4444) ALFA(ALF), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTC ALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ C(TRUQNTCALF)-I) /REALCMAXREP)) * 100. 

WRITEC1,4445) RSHATCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

620 END IF 

400 CONTINUE 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** ** *A *A ;V,v rshTBR CBRIDGE) **** ** * *^'r* 


IFCK.EQ.5) THEN 
DO 401 ALF=1,MAXALF 
TRUQNTCALF) = 0 

DO 501 REPS=1, MAXREP 

DIFFCREPS) = RSBRDG - RSHTBRCALF,REPS) 

501 CONTTINUE 

DO 601 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,'C " ” 'TRUE CONFIDENCE LIMIT IS: " , 
+ F8.4)') 

+ (TRUQNTCALF) / REAL(MAXREP)) * 100. 

GO TO 621 

ELSEIFCDIFF(REPS).LT.0.) THEN 
TRUQNTCALF) = REPS 
GO TO 611 

ELSE 
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END IF 

601 CONTINUE 

611 IF(TRUQNTCALF).EQ.0.) THEN 

WRITE(1,A443) ALFA(ALF) 

WRITE(1,’(” •',/"THE SMALLEST", 

+ " DIFFERENCE BETWEEN RSBRDG AND RSHTBR IS:’', 

+ F10.5)') DIFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ. 1. ) THEN 
WRITEC1,4442) ALFA(ALF) 

WRITE(1,'(" ",/"ALL RSHTBR", 

+ "ARE GREATER THAN RSBRDG")’) 

ELSEIF(ABSCDIFF(TRUQNT(ALF))).LE.ABS(DIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITE(1,4444) ALFA(AU), 

+ (TRUQNT(ALF) / REAL(MAXREP)) * 100. 

WRITEC1,4449) RSHTBR(ALF,TRUQNT(ALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITEC1,4449) RSHTBRCALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

621 END IF 

401 CONTINIJE 
ELSE 

END IF ■ 


***FIND THE TRUE CONFIDENCE LEVEL OF THE SYSTEM REL. ESTIMATE*** 
*** * * i e*rk**ie RHTSTR **** * *** *** 


IFCFLAG.EQ. 1) THEN 
DO 4400 ALF=i,MAXALF 
TRUQNTCALF) = 0 

DO 5500 REPS=1, HAXREP 

DIFFCREPS) = RS - RHTSTRCALF,REPS) 

5500 CONTINUE 

DO 6600 REPS=1, MAXREP 

IFCABSCDIFFCREPS)).LE.EPS) THEN 
TRUQNTCALF) = REPS 

WRITEC1,’( " " 'TRUE CONFIDENCE LIMIT IS: " , 
+ F8.4)‘) 

+ (TRUQNTCALF) / REALCMAXREP)) * 100. 

GO TO 6620 

ELSEIFCDIFF(REPS). LT. 0. ) THEN 
TRUQNTCALF) = REPS 
GO TO 6610 

ELSE 


END IF 

6600 CONTINUE 

6610 IF(TRUQNT(ALF). EQ. 0. ) THEN 

WRITEC1,4443) ALFACALF) 

WRITECl,’( " " The SMALLEST'', 

+ ’' DIFFEP^NCE between RS AND RHTSTR IS: " , 

+ F9.5)') DiFF(MAXREP) 

ELSEIF(TRUQNT(ALF).EQ.1.) THEN 
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WRITE(1,4442) ALFA(ALF) 

WRITEC:!,' ( '' " ./"ALL RHTSTR'', 

+ '' ARE GREATER THAN RS" )' ) 

ELSEIF(ABS(=DIFF(TRUQNT(ALF))).LE.ABSCDIFF(TRUQNT(ALF) - 1))) 
+ THEN 

WRITEtl,4444) ALFA(ALF), 

+ (TRUQNT(ALF) / REAL(MAXREF)) * 100. 

WRITEC1,4448) RHTSTR(ALF,TRUQNT(ALF)) 

WRITEC1,4446) 

ELSE 

WRITEC1,4444) ALFACALF), 

+ CCTRUQNTCALF)-!) / REALCMAXREP)) * 100. 

WRITE C1,4448) RHTSTRC ALF,TRUQNTCALF)-1) 

WRITEC1,4447) 

6620 END IF 

4400 CONTINUE 
ELSE 
END IF 

***PRINT THE ARRAYS PERTINENT TO THE OUPUT OF EACH REPLICATION*** 


IFCPRNTiEQ. 1) THEN 
I = 1 

185 WRITEC 1,-REPSHD) ALFACSELCTA), ALFACSELCTA), 

+ALFACSELCTB).ALFACSELCTB).ALFACSELCTA).ALFACSELCTA).ALFACSELCTB), 
+ALFACSELCTB) 

175 IFCI.GE. CMAXREP + 1)) THEN 
GOTO 180 

ELSE 

IFC CI. EQ. 71). OR. Cl. EQ. 211). OR. Cl. EQ. 351). OR. Cl. EQ. 491). OR. 

+ CI.EQ.631).OR.CI.EQ.771).OR.CI.EQ.911).OR.CI.EQ.1051) ) THEN 

I = I + 70 

WRITEC1,'C"+' ')') 

GOTO 185 

ELSE 

WRITEC1,3336) I, INTCDEGFRCl)), CHISQCl.I), QHTUPRCl,!), 

+ CHISQC2,I), QHTUPRC2,I) 

END IF 

IFCI + 70. LE. MAXREP) THEN 

WRITEC1,3337) 1+70,INTCDEGFRC1+70)),CHISQC1,1+70), 

+ QHTUPRC1,1+70),CHISQC 2,1+7 0),QHTUPRC 2,1+70) 

ELSE 
END IF 
1 = 1 + 1 
GOTO 175 
180 END IF 
ELSE 
ENDIF 

9999 WRITEC1,'C"THE TOTAL NO OF REPS WAS;", 18)') TOTREP 

WRITEC1,'C'’THE TOTAL NO OF EFFECTIVE REPS WAS: ".18)') LOOP 
WRITEC 1,'C "THE TOTAL NO OF NO FAILURE RUNS WAS: '^I8)') ZFAILS 
WRITEC 1,'C"average NO. OF COMPONENTS PER REPLICATION WITH ", 
+"N0 FAILURES; ",F5. 2)') ZFPREP / MAXREP 
WRITEC1,'C"THE TOTAL NO OF RUNS WITH FAILURES WAS: ",I8)’) FAILS 
0008 FORMAT C/ 3X,'C l’,5X,'C 2', 
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+'C 12',4X,'C 13',4X,'C ^ .4X,'C 15') 

0009 F0RMAT(/1X,'REP NO^,2X,IX,'CHISQR(',F4.3,')',1X, 

+'QHTUPRC',F4. 3,')',IX,'CHISQR(',F4. 3,'j',IX,'QHTUPR(',F4.3,')', 
+2X,' REP NO', 2X,' DF', IXCHISQR(',F4. 3, ’ )', IX, 

+'QHTUPRC’,F4. 3,')',IX,'CHISQR(',F4. 3,')',IX,'qHTUPR(',F4. 3,')'/) 
0001 FORMAT (///'UNIFORM RANDOM DEVIATES ARE:') 

0002 FORMAT (///'BERNOULLI TRIALS ARE:') 

0003 FORMAT (///'TOTAL NO. OF FAILURES FOR EACH COMPONENT: ') 

0004 FORMAT (///'ESTIMATED UNRELIABILITY FOR EACH COMPONENT:’) 

0005 FORMAT (///'TOTAL NUMBER OF MISSION TESTS:') 

0006 FORMAT (///'ESTIMATED WEIGHTS FOR EACH COMPONENT: ' ) 

0007 FORMAT C///'Q I FOR EACH COMPONENT: ') 

1111 FORMAT (15F8.5) 


2222 FORMAT 

3333 FORMAT 

3334 FORMAT 

3335 FORMAT 

3336 FORMAT 

3337 FORMAT 

4442 FORMAT 
+'LIMIT 

4443 FORMAT 
+'LIMIT 

4444 FORMAT 
+'LIMIT 

4445 FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

5556 FORMAT 

5557 FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


C/1X,15CI4,4X)) 

(/IX,15(14,4X)) 

(/15F8.5) 

(A1X,15(I4,4X)) 

(T3,I4,T9,I3,T13,F11. 5,T27,F8.5,T39,F11.5,T53,F8. 5) 
('+‘,T67,I4,T73,I3,T77,F11.5,T91,F8.5,T103,F11.5,T117,F8,5) 
(' ',///'the RESULTING (1- ',F4.! 


'■,///'the resulting 

« ipcr, I 


4446 

4447 

4448 

4449 
5555 


5656 

5666 

5657 

5667 
5755 


IS: ',T50 
(' ' 

IS: ',T50 
(' ‘ 

IS: 

(’ 


00.000 ') 


100 . 0000 ') 


',///'the resulting 

' ,T50,F8..4) 

I • 


3,') CONFIDENCE 


(1 - ',F4. 3,') CONFIDENCE 
(1 - '.,F4. 3,') CONFIDENCE 


,/'THE RSHAT VALUE 
'./'(FIRST NEGATIVE 


,T51,F8.5) 


CLOSEST TO RS IS: 

DIFFERENCE)') 

',/'(El£MENT PRECEEDING FIRST NEGATIVE DIFFERENCE)') 

' ,/'THE RHTSTR VALUE CLOSEST TO RS IS: ',T51,F8.5) 
',/'THE RSHTBR VALUE CLOSEST TO RSBRDG IS: ',T51,F8.5) 


(' ',///'THE ',I4,'(1-',F4.3,') QUANTILE IS: ',T49,F8.3) 


(' 


' ,/'THE 
' ,/'THE 

I 


',/'THE 
',/'THE 


VALUE OF RSHAT FOR THAT QUANTILE IS:',T51,F8.5) 
DIFFERENCECRS - RSHAT) IS;',T51,F8.5) 

VALUE OF RHTSTR FOR THAT QUANTILE IS:',T51,F8.5) 
VALUE OF RSHTBR FOR THAT QUANTILE IS:',T51,F8.5) 
DIFFERENCECRS - RHTSTR) IS:',T51,F8.5) 


(' ',/'THE DIFFERENCECRS - RSHTBR) IS: ',T51,F8.5) 
(' ‘ 


',///'SINCE THE NO. OF MISSION TESTS IS THE SAME FOR', 

LL COMPi 

'"RHTSTR" IS COMPUTED') 


+' ALL COMPONENTS THE CLOSED FORM SERIES SYSTEM RELIABILITY 


6666 FORMAT ( *4-* , ^ 

RUN INPUT^SETTINGS 

',//'NUMBER OF REPLICATIONS; ',T50,I4) 

'.//'NUMBER OF COMPONENTS; ’,T50.I4) 

'.//'SYSTEM RELIABILITY FUNCTION;'.T50.'SERIES') 
'.//’SYSTEM RELIABILITY FUNCTION;'.T50.'BRIDGE') 
'.//'MASTER UNRELIABILITY USED:’.T50.F8.5) 

' .//'INPUT WEIGHTSCA SUB l"S);') 

('+’,///’’ 

N R E S U 


6667 

6668 

6669 
6699 

6670 

6671 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 

6674 FORMAT 
+'**R U 

+'**************') 

6675 FORMAT (’+’.///' 


(’ 

(' 

(’ 

(’ 




4 : -h 


' ESTIMATE EI^ORS *********************’V*****************aa***a'^ 
'aaaaaAaaaaaaaa') 

6676 FORMAT (’+',///'aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa' 
+' TRUE CONFIDENCE LIMITS aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa^ 

+'AAAAAAAAAAAAAA') 

END 
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APPENDIX I. SUBROUTINES 

IMSL ROUTINE NAME - USMNMX USMNOOlO 

USMN0020 

-.-.USMN0030 

USMN0040 

COMPUTER - IBM/SINGLE USMN0050 

USMN0060 

LATEST REVISION - JANUARY 1, 1978 USMN0070 

USMNdOSO 

PURPOSE - DETERMINATION OF THE MINIMUM AND MAXIMUM USMN0090 

VALUES OF A VECTOR USMNOlOO 

USMNOllO 

USAGE - CALL USMNIffi (X,N,INC,XMIN,XMAX) USMN0120 

USMN0130 

ARGUMENTS X - INPUT VECTOR OF LENGTH N FROM WHICH MINIMUM, USMN0140 

MAXIMUM VALUES ARE TO BE TAKEN. USMN0150 

N - LENGTH OF THE INPUT VECTOR X. (INPUT) USMN0160 

INC - DISPLACEMENT BETWEEN CONSECUTIVE VALUES OF X USMN0170 

TO BE CONSIDERED. USMN0180 

XMIN - OUTPUT SCALAR CONTAINING MINIMUM VALUE OF X. USMN0190 

XMAX - OUTPUT SCALAR CONTAINING MAXIMUM VALUE OF X. USMN0200 

USMN0210 

PRECISION/HARDWARE - SINGLE AND DOUBLE/H32. USMN0220 

- SINGLE/H36,H48,H60 USMN0230 

USMN0240 

REQD. IMSL ROUTINES - NONE REQUIRED USMN0250 

USMN0260 

NOTATION - INFORMATION ON SPECIAL NOTATION AND USMN0270 

CONVENTIONS IS AVAILABLE IN THE MANUAL USMN0280 

INTRODUCTION OR THROUGH IMSL ROOTINE UHELP USMN0290 

USMN0300 

COPYRIGHT - 1978 BY IMSL, INC. ALL RIGHTS RESERVED. USMN0310 

USMN0320 

WARRANTY - IMSL WARRANTS ONLY THAT IMSL TESTING HAS BEEN USMN0330 

APPLIED TO THIS CODE. NO OTHER WARRANTY, USMN0340 
EXPRESSED OR IMPLIED, IS APPLICABLE. USMN0350 


USMN0360 
USMN0370 
USMN0380 

SUBROUTINE USMNMX (X,N,INC,XMIN,XMAX) USMN0390 

C USMN0400 

DIMENSION XCN) USMN0410 

C FIRST EXECUTABLE STATEMENT USMN0420 

XMIN = X(l) USMN0430 

XMAX = X(l) USMN0440 

DO 10 1=1,N,INC USMN0450 

IF (X(I) .GE. XMIN) GO TO 5 USMN0460 

XMIN = X(I) USMN0470 

GO TO 10 USMN0480 

5 IF (X(I) .GT. XMAX) XMAX = X(I) USMN0490 

10 CONTINUE USMN0500 

RETURN USMN0510 
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END 


A. 


IDENTIFICATION: 

TITLE: NUMERICAL SORT 

ID: Ml-NPG-SHSORT (F-IV) 

PROGRAMMER: R. BRUNELL 

DATE: MARCH 1968 

MODIFIED: DEC. 1973 BY L. NOLAN; 


B. 


D. 


USMN0520 

SORTOOlO 

SORT0020 

SORT0030 

S0RT0040 

SORT0050 

SORT0060 

SORT0070 

S0RT0080 

S0RT0090 

SORTOlOO 

SORTOllO 


PURPOSE: 

TO SORT, IN ASCENDING ORDER, AN ARRAY OF SINGLE PRECISION REAL 
NUMBERS BY THE METHOD OF SHELL, AND TO PRODUCE AN ARRAY OF INDEXESS0RT0120 
SO USER CAN RE-ORDER OTHER CORRESPONDING INFORMATION ACCORDING TO S0RT0130 


ASCENDING VALUES OF "A". 
USAGE: 

1. CALLING STATEMENT: 
CALL SHSORT(A,KEY,N) 
ARGUMENTS: 


2 . 


A - ARRAY OF NUMBERS TO BE SORTED. THIS ARRAY IS SORTED 
(RE-ORDERED) BY "SHSORT". 

KEY - ARRAY, DIMENSIONED AT LEAST N IN CALLING PROGRAM, TO BE 
FILLED BY USER WITH INTEGERS FROM i TO N. AFTER EXIT 
FROM SHSORT, KEY(l) WILL CONTAIN THE ORIGINAL INDEX OF 
IHE SMALLEST ELEMENT OF "A"; KEYC2) WILL CONTAIN THE 
ORIGINAL INDEX OF THE NEXT-TO-SMALLEST ELEMENT OF "A"; 
ETC. KEY(N) WILL CONTAIN THE ORIGINAL INDEX OF THE 
LARGEST ELEMENT OF "A". 

N - NUMBER OF MEMBERS IN ARRAYS "A" AND "KEY". 

REFERENCES: 

1. "ALGORITHM 201, SHELLSORT", BOOTHROYD, J. , "COMMUNICATIONS OF 


SORT0140 

SORT0150 

SORT0160 

S0RT0170 

S0RT0180 

SORT0190 

SORT0200 

S0RT0210 

S0RT0220 

SORT0230 

SORT0240 

SORT0250 

S0RT0260 

SORT0270 

S0RT0280 

SORT0290 

SORT0300 

SORT0310 

SORT0320 


ACM", VOL 6, NO 8, AUGUST 1963, P.445. 

SORT0330 

2. "CERTIFICATION OF ALGORITHM 

201", BATTY,M.A., "COMMUNICATIONS 

S0RT0340 

OF ACM", VOL 7, NO 6, JUNE 

1964, P.349. 

SORT0350 



SORT0360 

SUBROUTINE SHSORT(A,KEY,N) 


SORT0370 

DIMENSION A(N),KEY(N) 


SORT0380 

Ml=l 


SORT0390 

6 M1=M1*2 


SORT0400 

IF (Ml .LE. N) GO TO 6 


SORT0410 

Ml=Ml/2-l 


SORT0420 

MM=MAX0(M1/2,1) 


SORT0430 

GO TO 21 


S0RT0440 

20 MM=MM/2 


SORT0450 

IF (MM .LE. 0) GO TO 100 


SORT0460 

21 K=N-MM 


S0RT0470 

22 DO 1 J=1,K 


SORT0480 

II=J 


S0RT0490 

11 IM=II+MM 


S0RT0500 

IF (A(IM) .GE. A(II)) GO TO 1 


SORT0510 

TEMP=A(II) 


SORT0520 

IT=KEY(II) 


S0RT0530 


i 
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A(IT)=A(IM) 

SORT0540 

KEY(II)=KEYCIM) 

S0RT0550 

A(IM)=TEMP 

S0RT0560 

KEY(IM)=IT 

S0RT0570 

II=II-MM 

S0RT0580 

IF (II .GT. 0) GO TO 11 

S0RT0590 

1 CONTINUE 

SORT0600 

GO TO 20 

SORT0610 

100 RETURN 

SORT0620 

END 

SORT0630 


SUBROUTINE RHTSRS(QHTUP,AAHTI,N,RRSHAT) 

***//. THIS ROUTINE CALCULATES THE VALUE OF THE SYSTEM RELIABILITY OF A 
***// SERIES SYSTEM OF'N’ NO. OF COMPONENTS WHICH HAVE UNRELIABILITY 
***// ’OHTUP’. THE FINAL SYSTEM RELIABILITY VALUE PASSED IS ’RRSHAT’ 
REAL*4 QHTUP, RRSHAT, AAHTI(N) 

INTEGER I, N 
RRSHAT = I. 

DO 100 1=1, N 

RRSHAT = RRSHAT * (1 -(AAHTI(T)* QHTUP)) 

100 CONTINUE 
END 


SUBROUTINE RSRS(QIS,N,RRS) 

***// THIS ROUTINE CALCULATES THE VALUE OF THE SYSTEM RELIABILITY OF A 
***// SERIES SYSTEM OF ’COMP’ NO. OF COMPONENTS WHICH HAVE UNRELIABILITY 
***//; 'Qis’.. THE FINAL SYSTEM RELIABILITY VALUE PASSED IS ’RRS’ 

REAL*4 QIS(N), RRS 
INTEGER I, N 
RRS = 1. 

DO 100 1=1, N 

RRS = RRS * (1 - QIS(I)) 

100 CONTINUE 
END 


***// THIS SUBROUTINE CALCULATES THE ESTIMATED RELIABILITY OF A 
***// 5-COMPONENT BRIDGE STRUCTURE. (ONLY CARRIED OUT TO THE Q-CUBED TERM 
SUBROUTINE RHTBRG(QHTUP,AHT,N,RRBRDG) 

REAL*4 QHTUP, RRBRDG, AHT(N) 

INTEGER N 

* PRINT *, 'THE VALUES FOR AHAT PASSED ARE:',AHT 
RRBRDG=1-((QHTUP**2)*(AHT(1)*AHT(2)+AHT(4)*AHT(5)))- 

CC(QHTUP**3)*(AHT(1)*AHT(3)*AHT(5)+AHT(2)*AHT(3)*AHT(4)))+ 
C(CQHTUP**4)*(AHT(1)*AHT(2)*AHT(3)*AHT(4)+AHT(1)*AHT(2)*AHT(3)* 
CAHT(5)+AHT(1)*AHT(2)*AHT(4)*AHT(5)+AHT(1)*AHTC3)*AHTC4)*AHT(5)+ 
CAHTC 2)*AHT(3)*AHT(4)*AHT(5)))- 

C2*((QHTUP**5)*(AHT(1)*AHT(2)*AHT(3)*AHT(4)*AHT(5))) 

* PRINT *, 'COMPUTED RRBRDG IS:', RRBRDG 
END 


SUBROUTINE CPARE(FI,K,BFLAG) 






INTEGER BFLAG,FI(K) 

BFLAG=0- 

IF ((FI(l) .EQ. 0) . AND. (FI(3) .EQ. 0) . AND. 

+ (FI(5) .EQ. 0)) THEN 

BFLAG=0 

ELSE IF- ((FI(1) .EQ. 0) .AOT). (FIC4) .EQ. 0)) THEN 
BFLAG =0 

ELSE IF ((FI(2) .EQ. 0) .AND. CFI(5) .EQ. 0)) THEN 
BFLAG = 0 

ELSE IF ((FI(2) .EQ. 0) .AND. (FIC3) .EQ. 0) .AND. 

+ (FI(4) .EQ. 0)) THEN 

BFLAG = 0 

ELSE IF (((FICl) .EQ. 1) .AND. (FI(2) .EQ. 1)) .AND. 

+ ((FIC2) .EQ. 0) .AND. (FI(4) .EQ. 0) .AND. 

+ (FI(5) .EQ. 0))) THEN 

BFLAG = 1 

ELSE IF (((FI(4) .EQ. 1) .AND. (FI(5) .EQ. 1)) .AND. 

+ ((FI(l) .EQ. 0) .AND. (FI(2) .EQ. 0) .AND. 

+ (FI(3) .EQ. 0))) THEN 

BFLAG = 1 

ELSE IF (C(FIG2) .EQ. 0) .AND. (FI(4) .EQ. 0)) .AND. 

+ ((FI(l) .EQ. 1) .AND. (FI(3) .EQ. 1) .AND. 

+ (FIC5) .EQ. 1))) THEN 

BFLAG = 1 

ELSE IF (((FI(1) .EQ. 0) .AND. (FI(5) .EQ. 0)) .AND. 

+■ ((FI(2) .EQ. 1) .AND. (FI(3) .EQ. 1) .AND. 

+ CFI(4) .EQ. 1))) THEN 

BFLAG = 1 

■ELSE 

BFLAG = 2 

ENDIF 

RETURN 

END 


***// THIS SUBROUTINE CALCULATES THE "TRUE" RELIABILITY OF A 5-COMPONENT 
***// BRIDGE STRUCTURE 

SUBROUTINE RBRIDG(QI,N,RRSS) 

REAL*4 QI(N), RRSS 
INTEGER N 
IFCN.NE.S) THEN 

WRITE(1,’(’'WARNING: BRIDGE STRUCTURE ONLY USES ", 

+''THE FIRST 5 COMPONENTS'')') 

ELSE 
END IF 

RRSS=( 1-QI( 1) )*(1-QI(4))+(1-QI(2))*(1-QI(5))+(l-QKi))*(1-QI(3))* 
C(1-QI(5))+(1-QIC2))*(1-QI(3))*(1-QI(4))-(1-QI(1))*(1-QIC2))* 
C(1-QI(3))*(1-QIC4))-(1-QI(1))*(1-QI(2))*(1-QI(3))*(1-QI(5))- 
C(1-QI(T))*(1-QI(2))*(1-QI(4))*(1-QI(5))-(1-QI(1))*(1-QI(3))* 
C(1-QI(4))*(1-QI(5))-(1-QI(2))*(1-QI(3))*(1-QI(4))*(1-QI(5))+ 
C2*(1-QI(1))*(1-QI(2))*(1-QI(3))*(1-QI(4))*(1-QI(5)) 

END 
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